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Abstract  (cont) 

The  seismic  arrays  operated  by  the  Norwegian  NDC  comprise  the  Norwegian  Seismic  Array 
(NOA),  the  Norwegian  Regional  Seismic  Array  (NORES),  the  Arctic  Regional  Seismic  Array 
(ARCES)  and  the  Spitsbergen  Regional  Array  (SPITS).  This  report  also  presents  statistics  for 
additional  seismic  stations  which  through  cooperative  agreements  with  institutions  in  the  host 
countries  provide  continuous  data  to  the  NORSAR  Data  Processing  Center  (NPDC).  These  sta¬ 
tions  comprise  the  Finnish  Regional  Seismic  Array  (FINES),  the  Hagfors  array  in  Sweden  and 
the  regional  seismic  array  in  Apatity,  Russia. 

The  NOA  Detection  Processing  system  has  been  operated  throughout  the  period  with  an  average 
uptime  of  100.00%.  A  total  of  2004  seismic  events  have  been  reported  in  the  NOA  monthly  seis¬ 
mic  bulletin  from  July  through  December  2001 .  On-line  detection  processing  and  data  recording 
at  the  NDC  of  NORES,  ARCES  and  FINES  data  have  been  conducted  throughout  the  period. 
Data  from  two  small-aperture  arrays  at  sites  in  Spitsbergen  and  Apatity,  Kola  Peninsula,  as  well 
as  the  Hagfors  array  in  Sweden,  have  also  been  recorded  and  processed.  Processing  statistics  for 
the  arrays  for  the  reporting  period  are  given. 

A  summary  of  the  activities  related  to  the  GSETT-3  experiment  and  experience  gained  at  the 
Norwegian  NDC  during  the  reporting  period  is  provided  in  Section  4.  Norway  is  now  contribut¬ 
ing  primary  station  data  from  two  seismic  arrays;  ARCES  and  NOA  and  one  auxiliary  array 
(SPITS).  These  data  are  being  provided  to  the  IDC  via  the  global  communications  infrastructure 
(GCI).  Continuous  data  from  all  three  arrays  are  in  addition  being  transmitted  to  the  US  NDC. 
The  performance  of  the  data  transmission  to  the  US  NDC  has  been  satisfactory  during  the 
reporting  period. 

The  PrepCom  has  encouraged  states  that  operate  IMS-designated  stations  to  continue  to  do  so 
on  a  voluntary  basis  and  in  the  framework  of  the  GSETT-3  experiment  until  the  stations  have 
been  certified  for  formal  inclusion  in  IMS.  So  far  among  the  Norwegian  stations,  the  NOA  and 
the  ARCES  array  (PS27  and  PS28  respectively)  have  been  certified.  We  envisage  continuing  the 
provision  of  data  from  these  and  other  Norwegian  IMS-designated  stations  in  accordance  with 
current  procedures. 

Summaries  of  six  scientific  and  technical  contributions  are  presented  in  Chapter  6  of  this  report. 

Section  6.1  contains  a  study  of  the  estimated  global  and  regional  detection  capability  of  the  IMS 
primary  seismic  network.  The  study  makes  use  of  the  Threshold  Monitoring  (TM)  method.  The 
TM  method  is  capable  of  using  actual  seismic  data  for  a  given  time  interval  as  the  basis  for  the 
detectability  calculations.  In  cases  when  a  seismic  station  did  not  provide  data  during  the  time 
period  under  study,  an  estimated  background  noise  level  can  be  assigned.  These  noise  estimates 
and  signal/noise  variances  can  be  based  on  results  from  earlier  studies,  or  they  can  be  taken  from 
stations  assumed  to  have  similar  noise  and  signal  characteristics. 

We  have  found  that  the  35  stations  of  the  primary  seismic  IMS  network  operational  as  of  July 
2001  have  a  three-station  detection  capability  that  is  quite  close  to  the  projected  performance  of 
the  49  stations  of  the  network.  The  largest  difference  of  about  0.5  m^  units  is  found  in  southwest 
Asia.  This  can  be  explained  by  the  fact  that  most  of  the  sensitive  array  stations  of  the  network 
are  already  in  place.  Another  uncertainty  is  the  noise  levels  assigned  to  the  planned  stations  and 
to  the  stations  that  did  not  provide  real  noise  data  during  the  period  under  investigation.  The 
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assignment  of  realistic  noise  levels  for  planned  stations  is  therefore  a  topic  requiring  additional 
studies. 

This  study  also  confirms  the  results  provided  by  Kvaema  and  Ringdal  (1999)  that  large  earth¬ 
quakes  and  the  corresponding  coda  energy  can  temporarily,  over  tens  of  minutes,  significantly 
reduce  the  detection  capability  of  the  IMS  network.  During  such  conditions  the  use  of  high-fi:e- 
quency  regional  data,  combined  with  less  stringent  event  formation  criteria  (1  or  2  station  only) 
will  be  important  for  CTBT  verification  purposes. 

Section  6.2  is  a  study  of  travel  times  and  attenuation  relations  of  regional  seismic  phases  in  the 
Barents  Sea  region.  A  database  containing  42  events  in  the  Barents  Sea  region  has  been  compiled 
and  analyzed  with  the  aim  of  evaluating  crustal  models,  travel  times  and  attenuation  relations  in 
the  context  of  performing  regional  detection  threshold  monitoring  of  this  region.  The  42  events 
are  mostly  located  around  the  circumference  of  the  study  area  due  to  the  virtually  aseismic  nature 
of  the  Barents  Sea  itself.  Regional  and  Sn  phases  were  observable  for  most  events  in  the  data¬ 
base,  while  Pg  and  Lg  phases  were  only  observable  for  events  with  ray  paths  within  continental 
crust.  This  corroborates  a  number  of  previous  observations  of  Lg- wave  blockage  across  the  Bar¬ 
ents  Sea.  Three  existing  velocity  models  were  evaluated,  with  a  model  having  slightly  lower  S- 
velocities  than  earlier  assumed  in  the  upper  mantle  giving  the  overall  best  fit  to  the  observed  arriv¬ 
als.  In  order  to  estimate  magnitudes,  short  term  average  (STA)  and  spectral  amplitude  values  were 
calculated  in  several  frequency  bands  for  aU  phase  arrivals  in  the  data  base.  There  were  no  signif¬ 
icant  differences  between  spectral  and  STA  amplitudes,  so  the  latter  were  used  as  this  parameter  is 
more  efficient  to  calculate  in  real-time  processing.  An  inversion  was  performed  in  order  to  deter¬ 
mine  a  Pn  and  Sj,  attenuation  relation  specific  for  this  region.  The  resulting  magnitudes  based  on 
Pjj  and  Sjj  phases  gave  an  internally  consistent,  reasonably  stable  set  of  values,  which  can  be  cali¬ 
brated  towards  any  existing  global  or  regional  scale.  An  attenuation  relation  was  also  determined 
for  the  Lg  phase,  but  the  low  number  of  amplitude  readings  in  this  case  renders  the  results  less 
reliable. 

Section  6.5  is  a  study  of  experimental  Site-Specific  Threshold  Monitoring  (SSTM)  applied  to  the 
Lop  Nor  test  site  in  China.  This  smdy  is  a  follow-up  on  similar  experimental  studies  applied  to 
Novaya  Zemlya  and  the  “Kursk”  accident  site  in  the  Barents  Sea.  These  two  previous  case  studies 
have  been  reported  upon  in  earlier  Semiannual  Reports,  and  have  made  use  of  seismic  stations  at 
regional  distances,  mostly  from  the  Fennoscandian  array  network,  hi  the  case  of  Lop  Nor,  we 
have  appUed  a  combination  of  3-component  stations  and  arrays,  at  both  regional  and  teleseismic 
ranges.  Our  efforts  so  far,  as  reported  in  this  contribution,  comprise  mainly  a  study  of  available 
seismic  stations,  selected  among  those  stations  which  are  most  sensitive  to  seismic  events  in  the 
Lop  Nor  general  area,  and  tuning  of  the  signal  parameters  of  these  stations  so  as  to  prepare  pro¬ 
cessing  recipes  for  the  application  of  the  threshold  monitoring  tool. 

We  have  selected  a  total  of  17  seismic  stations,  including  both  IMS  stations  (mostly  arrays  at  tele- 
seismic  distances)  and  other  available  stations  (mostly  of  the  three-component  type,  and  situated 
in  the  regional  distance  range).  It  appears  that  the  combined  threshold  monitoring  level  of  this  sta¬ 
tion  network  is  close  to  or  better  than  m^  3.0  for  the  Lop  Nor  test  site.  This  is  quite  encouraging, 
taking  into  account  that  the  IMS  network  is  not  yet  well  developed  in  central  and  eastern  Asia. 
Once  the  full  IMS  network  becomes  available,  this  capability  should  improve  considerably. 

Section  6.4  is  a  study  of  seismicity  of  the  Spitsbergen  archipelago.  Spitsbergen  and  the  adjacent 
areas  are  parts  of  a  geologically  complex  region  with  moderate  to  high  earthquake  activity.  The 
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main  seismicity  in  the  area  is  associated  with  the  North- Atlantic  Ridge,  and  especially  the 
Knipovich  Ridge.  In  addition,  some  coal  mines  are  located  in  the  area  of  Spitsbergen,  causing 
occasional  induced  seismicity.  We  have  compared  maps  of  seismic  events  in  and  near  Spitsber¬ 
gen  for  the  period  1964  - 1998  as  taken  from  ISC  bulletins  with  a  corresponding  map  for  1998- 
2001  using  the  solutions  from  the  NORSAR  Reviewed  Regional  Seismic  Bulletin.  Even  though 
the  NORSAR  regional  bulletin  covers  a  much  shorter  time  interval  than  the  ISC  bulletin,  the 
observed  seismdcity  patterns  in  the  two  plots  are  remarkably  similar.  In  particular,  a  segment  of 
the  mid-atlantic  ridge  to  the  west  can  be  observed,  and  pronounced  seismic  zones  in  Nordaust- 
landet  (northeast)  and  Heerland  (southeast)  can  also  be  identified  on  both  maps.  We  have  pur¬ 
sued  further  the  location  calibration  of  the  Spitsbergen  area,  and  have  adapted  the  Barents 
regional  velocity  model  to  this  area  by  adding  a  sedimentary  layer  of  2  km.  This  has  resulted  in 
encouraging  improvements  in  locating  events  with  known  (ground  truth)  locations. 

Section  (5.5  is  a  comparison  of  various  location  procedures  applied  to  the  Kara  Sea  seismic  event 
of  16  August  1997.  Thrs  event  has  been  the  subject  of  considerable  discussion,  and  in  particular 
it  has  been  difficult  to  obtain  a  reliable  focal  depth  estimate  from  available  seismic  recordings. 
We  have  undertaken  a  sequence  of  location  experiments  to  compare  the  results  of  a)  using  differ¬ 
ent  velocity  models  to  describe  the  travel  times  of  the  phases  and  b)  to  make  a  comparison 
between  the  use  of  a  linearized  location  algorithm  (HYPOSAT  -  Schweitzer,  2001)  and  a  fully 
non-linear  scheme  (shakeNA  —  Sambiidge  and  Kennett,  2001).  For  direct  comparisons  between 
the  two  methods  we  have  used  a  standard  least-squares  misfit  criterion,  but  we  have  also  exam¬ 
ined  the  influence  of  more  robust  choices  for  data  misfit  when  using  the  non-linear  location 
scheme. 

This  study  has  shown  both  the  importance  of  S  wave  information  in  assessing  the  depth  of 
regional  events,  and  the  need  to  get  a  good  regional  velocity  model  for  both  P  and  S  in  order  to 
place  the  strongest  constraints  on  the  location  of  the  event. 

The  conclusions  from  our  experiment  of  comparing  the  different  data  centre  solutions  with  our 
results  are  that  using  only  a  limited  data  set  but  an  adequate  travel-time  model  one  can  locate  the 
event  in  the  Kara  Sea  relatively  close  to  our  HYPOS AT-based  location.  However,  in  this  case 
there  is  very  httle  depth  resolution.  The  relative  small  error  ellipses  are  a  problem  which  arises 
when  using  only  a  limited  number  of  data.  Then  the  data  errors  do  not  usually  follow  a  normal 
distribution  but  are  biased  in  one  direction  and  suggest  an  unreasonably  high  resolution  and 
accuracy. 

The  location  estimates  for  the  whole  data  set  from  the  different  techniques  agree  quite  well,  with 
some  overlap  of  the  estimated  confidence  regions.  The  study  concludes  that  the  event  carmot  be 
shallower  than  10  km  and  is  most  likely  in  the  lower  crust  around  20-30  km  depth.  Such  deep 
crustal  events  are  often  attributed  to  the  long-term  effects  of  ice-unloading  from  the  last  glacia¬ 
tion  and  have  previously  been  observed  at  Novaya  Zemlya  (e.g.  Marshall  et  al,  1989).  The 
assertion  that  this  event  was  not  an  underwater  explosion  is  supported  by  the  occurrence  of  an 
aftershock  in  the  same  epicentral  area  about  four  hours  later  and  the  fact  that  the  observed  seis¬ 
mic  signals  do  not  show  bubble  pulse  reverberations,  typical  for  underwater  explosions. 

Section  6.6  presents  some  results  derived  from  studying  the  seismic  signals  generated  by  the 
accident  of  the  Russian  submarine  “Kursk”  on  12  August  20(X).  Two  seismic  events,  which 
occurred  about  2  minutes  and  16  seconds  apart,  were  associated  with  this  accident.  The  first  of 
these  two  events  (Kursk- 1)  was  about  two  magnitude  units  smaller  than  the  second  one  (Kursk- 
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2),  which  had  a  local  magnitude  of  about  3.5  (Ringdal  et  al,  2000). 

To  get  a  better  understanding  of  the  accident,  we  have  tried  to  estimate  the  relative  location 
between  these  two  events.  We  have  applied  a  signal  correlation  analysis  procedure  to  measure, 
with  high  accuracy,  the  time  difference  between  the  two  events  for  various  phases  recorded  at 
six  seismic  stations.  We  then  applied  the  master-event  location  technique  and  inverted  the  time 
differences  to  obtain  a  relative  horizontal  location  between  the  two  events.  The  application  of 
the  master-event  analysis  suggests  that  the  submarine  moved  about  145  m  to  the  north-west 
during  the  135.8  s  between  the  two  events.  The  azimuth  of  this  movement  is  estimated  at  302°. 
After  the  accident  not  only  the  exact  position  of  the  Kursk  submarine  became  known  but  also 
the  direction  in  which  the  submarine  was  lying  on  the  sea  bottom.  This  direction  was  reported 
as  288°.  This  is  in  good  agreement  with  the  results  in  the  study  about  the  relative  movement  of 
the  submarine  during  the  time  interval  between  Kursk- 1  and  Kursk-2.  An  estimate  of  the  uncer¬ 
tainty  involved  in  the  calculations  is  also  provided. 
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1  Summary 

This  report  describes  the  research  activities  carried  out  at  NORSAR  under  Contract  No.  F08650- 
Ol-C-0055  for  the  period  1  July  -  31  December  2001.  In  addition,  it  provides  summary  informa¬ 
tion  on  operation  and  maintenance  (O&M)  activities  at  the  Norwegian  National  Data  Center 
(NDC)  during  the  same  period.  Research  activities  described  in  this  report,  as  well  as  transmis¬ 
sion  of  selected  data  to  the  United  States  NDC,  are  funded  by  the  United  States  Department  of 
Defense.  The  O&M  activities,  including  operation  of  transmission  links  within  Norway  and  to 
Vienna,  Austria  are  being  funded  jointly  by  the  CTBTO/PTS  and  the  Norwegian  Government, 
with  the  understanding  that  the  funding  of  all  IMS-related  activities  will  gradually  be  transferred 
to  the  CTBTO/PTS.  The  O&M  statistics  presented  in  this  report  are  included  for  the  purpose  of 
completeness,  and  in  order  to  maintain  consistency  with  earlier  reporting  practice. 

The  seismic  arrays  operated  by  the  Norwegian  NDC  comprise  the  Norwegian  Seismic  Array 
(NOA),  the  Norwegian  Regional  Seismic  Array  (NORES),  the  Arctic  Regional  Seismic  Array 
(ARCES)  and  the  Spitsbergen  Regional  Array  (SPITS).  This  report  also  presents  statistics  for 
additional  seismic  stations  which  through  cooperative  agreements  with  institutions  in  the  host 
countries  provide  continuous  data  to  the  NORSAR  Data  Processing  Center  (NPDC).  These  sta¬ 
tions  comprise  the  Finnish  Regional  Seismic  Array  (FINES),  the  Hagfors  array  in  Sweden  and 
the  regional  seismic  array  in  Apatity,  Russia. 

The  NOA  Detection  Processing  system  has  been  operated  throughout  the  period  with  an  average 
uptime  of  100.00%.  A  total  of  2004  seismic  events  have  been  reported  in  the  NOA  monthly  seis¬ 
mic  bulletin  from  July  through  December  2001 .  On-line  detection  processing  and  data  recording 
at  the  NDC  of  NORES,  ARCES  and  FINES  data  have  been  conducted  throughout  the  period. 
Data  from  two  small-aperture  arrays  at  sites  in  Spitsbergen  and  Apatity,  Kola  Peninsula,  as  well 
as  the  Hagfors  array  in  Sweden,  have  also  been  recorded  and  processed.  Processing  statistics  for 
the  arrays  for  the  reporting  period  are  given. 

A  summary  of  the  activities  related  to  the  GSETT-3  experiment  and  experience  gained  at  the 
Norwegian  NDC  during  the  reporting  period  is  provided  in  Section  4.  Norway  is  now  contribut¬ 
ing  primary  station  data  from  two  seismic  arrays:  ARCES  and  NOA  and  one  auxiliary  array 
(SPITS).  TTiese  data  are  being  provided  to  the  IDC  via  the  global  communications  infrastructure 
(GCI).  Continuous  data  from  all  three  arrays  are  in  addition  being  transmitted  to  the  US  NDC. 
The  performance  of  the  data  transmission  to  the  US  NDC  has  been  satisfactory  during  the 
reporting  period. 

The  PrepCom  has  encouraged  states  that  operate  IMS-designated  stations  to  continue  to  do  so 
on  a  voluntary  basis  and  in  the  framework  of  the  GSETT-3  experiment  until  the  stations  have 
been  certified  for  formal  inclusion  in  IMS.  So  far  among  the  Norwegian  stations,  the  NOA  and 
the  ARCES  array  (PS27  and  PS28  respectively)  have  been  certified.  We  envisage  continuing  the 
provision  of  data  from  these  and  other  Norwegian  IMS-designated  stations  in  accordance  with 
current  procedures. 

Summaries  of  six  scientific  and  technical  contributions  are  presented  in  Chapter  6  of  this  report. 

Section  6.1  contains  a  study  of  the  estimated  global  and  regional  detection  capability  of  the  IMS 
primary  seismic  network.  The  study  makes  use  of  the  Threshold  Monitoring  (TM)  method.  The 
TM  method  is  capable  of  using  actual  seismic  data  for  a  given  time  interval  as  the  basis  for  the 
detectability  calculations.  In  cases  when  a  seismic  station  did  not  provide  data  during  the  time 
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period  under  study,  an  estimated  background  noise  level  can  be  assigned.  These  noise  estimates 
and  signal/noise  variances  can  be  based  on  results  from  earlier  studies,  or  they  can  be  taken  from 
stations  assumed  to  have  similar  noise  and  signal  characteristics. 

We  have  found  that  the  35  stations  of  the  primary  seismic  IMS  network  operational  as  of  July 
2001  have  a  three-station  detection  capability  that  is  quite  close  to  the  projected  performance  of 
the  49  stations  of  the  network.  The  largest  difference  of  about  0.5 1115  units  is  found  in  southwest 
Asia.  This  can  be  explained  by  the  fact  that  most  of  the  sensitive  array  stations  of  the  network  are 
already  in  place.  Another  uncertainty  is  the  noise  levels  assigned  to  the  planned  stations  and  to  the 
stations  that  did  not  provide  real  noise  data  during  the  period  under  investigation.  The  assignment 
of  realistic  noise  levels  for  planned  stations  is  therefore  a  topic  requiring  additional  studies. 

This  study  also  confirms  the  results  provided  by  Kvaema  and  Ringdal  (1999)  that  large  earth¬ 
quakes  and  the  corresponding  coda  energy  can  temporarily,  over  tens  of  minutes,  significantly 
reduce  the  detection  capability  of  the  IMS  network.  During  such  conditions  the  use  of  high-fre¬ 
quency  regional  data,  combined  with  less  stringent  event  formation  criteria  (1  or  2  station  only) 
will  be  important  for  CTBT  verification  purposes. 

Section  6.2  is  a  study  of  travel  times  and  attenuation  relations  of  regional  seismic  phases  in  the 
Barents  Sea  region.  A  database  containing  42  events  in  the  Barents  Sea  region  has  been  compiled 
and  analyzed  with  the  aim  of  evaluating  crustal  models,  travel  times  and  attenuation  relations  in 
the  context  of  performing  regional  detection  threshold  monitoring  of  this  region.  The  42  events 
are  mostly  located  around  the  circumference  of  the  study  area  due  to  the  virtually  aseismic  nature 
of  the  Barents  Sea  itself.  Regional  and  Sn  phases  were  observable  for  most  events  in  the  data¬ 
base,  while  Pg  and  Lg  phases  were  only  observable  for  events  with  ray  paths  within  continental 
crust.  This  corroborates  a  number  of  previous  observations  of  Lg- wave  blockage  across  the  Bar¬ 
ents  Sea.  Three  existing  velocity  models  were  evaluated,  with  a  model  having  slightly  lower  S- 
velocities  than  earlier  assumed  in  the  upper  mantle  giving  the  overall  best  fit  to  the  observed  arriv¬ 
als.  In  order  to  estimate  magnitudes,  short  term  average  (STA)  and  spectral  amplitude  values  were 
calculated  in  several  frequency  bands  for  all  phase  arrivals  in  the  data  base.  There  were  no  signif¬ 
icant  differences  between  spectral  and  STA  amplitudes,  so  the  latter  were  used  as  this  parameter  is 
more  efficient  to  calculate  in  real-time  processing.  An  inversion  was  performed  in  order  to  deter¬ 
mine  a  Pj,  and  Sn  attenuation  relation  specific  for  this  region.  The  resulting  magnitudes  based  on 
Pn  and  Sj,  phases  gave  an  internally  consistent,  reasonably  stable  set  of  values,  which  can  be  cali¬ 
brated  towards  any  existing  global  or  regional  scale.  An  attenuation  relation  was  also  determined 
for  the  Lg  phase,  but  the  low  number  of  amplitude  readings  in  this  case  renders  the  results  less 
reliable. 

Section  6.3  is  a  study  of  experimental  Site-Specific  Threshold  Monitoring  (SSTM)  applied  to  the 
Lop  Nor  test  site  in  China.  This  study  is  a  follow-up  on  similar  experimental  studies  applied  to 
Novaya  Zemlya  and  the  “Kursk”  accident  site  in  the  Barents  Sea.  These  two  previous  case  studies 
have  been  reported  upon  in  earlier  Semiannual  Reports,  and  have  made  use  of  seismic  stations  at 
regional  distances,  mostly  from  the  Fennoscandian  array  network.  In  the  case  of  Lop  Nor,  we 
have  appUed  a  combination  of  3-component  stations  and  arrays,  at  both  regional  and  teleseismic 
ranges.  Our  efforts  so  far,  as  reported  in  this  contribution,  comprise  mainly  a  study  of  available 
seismic  stations,  selected  among  those  stations  which  are  most  sensitive  to  seismic  events  in  the 
Lop  Nor  general  area,  and  tuning  of  the  signal  parameters  of  these  stations  so  as  to  prepare  pro¬ 
cessing  recipes  for  the  application  of  the  threshold  monitoring  tool. 
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We  have  selected  a  total  of  17  seismic  stations,  including  both  IMS  stations  (mostly  arrays  at 
teleseismic  distances)  and  other  available  stations  (mostly  of  the  three-component  type,  and  sit¬ 
uated  in  the  regional  distance  range).  It  appears  that  the  combined  threshold  monitoring  level  of 
this  station  network  is  close  to  or  better  than  m^  3.0  for  the  Lop  Nor  test  site.  This  is  quite 
encouraging,  taking  into  account  that  the  IMS  network  is  not  yet  weU  developed  in  central  and 
eastern  Asia.  Once  the  full  IMS  network  becomes  available,  this  capability  should  improve  con¬ 
siderably. 

Section  6.4  is  a  study  of  seismicity  of  the  Spitsbergen  archipelago.  Spitsbergen  and  the  adjacent 
areas  are  parts  of  a  geologically  complex  region  with  moderate  to  high  earthquake  activity.  The 
main  seismicity  in  the  area  is  associated  with  the  North- Atlantic  Ridge,  and  especially  the 
Knipovich  Ridge.  In  addition,  some  coal  mines  are  located  in  the  area  of  Spitsbergen,  causing 
occasional  induced  seismicity.  We  have  compared  maps  of  seismic  events  in  and  near  Spitsber¬ 
gen  for  the  period  1964  - 1998  as  taken  from  ISC  bulletins  with  a  corresponding  map  for  1998- 
2001  using  the  solutions  from  the  NORSAR  Reviewed  Regional  Seismic  Bulletin.  Even  though 
the  NORSAR  regional  bulletin  covers  a  much  shorter  time  interval  than  the  ISC  bulletin,  the 
observed  seismicity  patterns  in  the  two  plots  are  remarkably  similar.  In  particular,  a  segment  of 
the  mid-atlantic  ridge  to  the  west  can  be  observed,  and  pronounced  seismic  zones  in  Nordaust- 
landet  (northeast)  and  Heerland  (southeast)  can  also  be  identified  on  both  maps.  We  have  pur¬ 
sued  further  the  location  calibration  of  the  Spitsbergen  area,  and  have  adapted  the  Barents 
regional  velocity  model  to  this  area  by  adding  a  sedimentary  layer  of  2  km.  This  has  resulted  in 
encouraging  improvements  in  locating  events  with  known  (ground  truth)  locations. 

Section  6.5  is  a  comparison  of  various  location  procedures  applied  to  the  Kara  Sea  seismic  event 
of  16  August  1997.  This  event  has  been  the  subject  of  considerable  discussion,  and  in  particular 
it  has  been  difficult  to  obtain  a  reliable  focal  depth  estimate  fi'om  available  seismic  recordings. 
We  have  undertaken  a  sequence  of  location  experiments  to  compare  the  results  of  a)  using  differ¬ 
ent  velocity  models  to  describe  the  travel  times  of  the  phases  and  b)  to  make  a  comparison 
between  the  use  of  a  linearized  location  algorithm  (HYPOS AT  -  Schweitzer,  2001)  and  a  fully 
non-linear  scheme  (shakeNA  -  Sambridge  and  Kennett,  2001).  For  direct  comparisons  between 
the  two  methods  we  have  used  a  standard  least-squares  misfit  criterion,  but  we  have  also  exam¬ 
ined  the  influence  of  more  robust  choices  for  data  misfit  when  using  the  non-linear  location 
scheme. 

This  study  has  shown  both  the  importance  of  S  wave  information  in  assessing  the  depth  of 
regional  events,  and  the  need  to  get  a  good  regional  velocity  model  for  both  P  and  S  in  order  to 
place  the  strongest  constraints  on  the  location  of  the  event. 

The  conclusions  from  our  experiment  of  comparing  the  different  data  centre  solutions  with  our 
results  are  that  using  only  a  limited  data  set  but  an  adequate  travel-time  model  one  can  locate  the 
event  in  the  Kara  Sea  relatively  close  to  our  HYPOS AT-based  location.  However,  in  this  case 
there  is  very  little  depth  resolution.  The  relative  small  error  ellipses  are  a  problem  which  arises 
when  using  only  a  limited  number  of  data.  Then  the  data  errors  do  not  usually  follow  a  normal 
distribution  but  are  biased  in  one  direction  and  suggest  an  unreasonably  high  resolution  and 
accuracy. 

The  location  estimates  for  the  whole  data  set  from  the  different  techniques  agree  quite  well,  with 
some  overlap  of  the  estimated  confidence  regions.  The  study  concludes  that  the  event  cannot  be 
shallower  than  10  km  and  is  most  likely  in  the  lower  crust  around  20-30  km  depth.  Such  deep 
crustal  events  are  often  attributed  to  the  long-term  effects  of  ice-unloading  from  the  last  glacia- 
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tion  and  have  previously  been  observed  at  Novaya  Zemlya  (e.g.  Marshall  et  al,  1989).  The 
assertion  that  this  event  was  not  an  underwater  explosion  is  supported  by  the  occurrence  of  an 
aftershock  in  the  same  epicentral  area  about  four  hours  later  and  the  fact  that  the  observed  seis¬ 
mic  signals  do  not  show  bubble  pulse  reverberations,  typical  for  underwater  explosions. 

Section  6.6  presents  some  results  derived  from  studying  the  seismic  signals  generated  by  the 
accident  of  the  Russian  submarine  “Kursk”  on  12  August  2000.  Two  seismic  events,  which 
occurred  about  2  minutes  and  16  seconds  apart,  were  associated  with  this  accident.  The  first  of 
these  two  events  (Kursk- 1)  was  about  two  magnitude  units  smaller  than  the  second  one  (Kursk- 
2),  which  had  a  local  magnitude  of  about  3.5  (Ringdal  et  al,  2000). 

To  get  a  better  understanding  of  the  accident,  we  have  tried  to  estimate  the  relative  location 
between  these  two  events.  We  have  applied  a  signal  correlation  analysis  procedure  to  measure, 
with  high  accuracy,  the  time  difference  between  the  two  events  for  various  phases  recorded  at 
six  seismic  stations.  We  then  applied  the  master-event  location  technique  and  inverted  the  time 
differences  to  obtain  a  relative  horizontal  location  between  the  two  events.  The  application  of 
the  master-event  analysis  suggests  that  the  submarine  moved  about  145  m  to  the  north-west 
during  the  135.8  s  between  the  two  events.  The  azimuth  of  this  movement  is  estimated  at  302°. 
After  the  accident  not  only  the  exact  position  of  the  Kursk  submarine  became  known  but  also 
the  direction  in  which  the  submarine  was  lying  on  the  sea  bottom.  This  direction  was  reported 
as  288°.  This  is  in  good  agreement  with  the  results  in  this  study  about  the  relative  movement  of 
the  submarine  during  the  time  interval  between  Kursk- 1  and  Kursk-2.  An  estimate  of  the  uncer¬ 
tainty  involved  in  the  calculations  is  also  provided. 

Frode  Ringdal 
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2  Operation  of  International  Monitoring  System  (IMS)  Stations 
in  Norway 

2.1  PS27  —  Primary  Seismic  Station  NOA 

The  average  recording  time  was  100%,  the  same  as  for  the  previous  reporting  period. 

Monthly  uptimes  for  the  NORSAR  on-line  data  recording  task,  taking  into  account  aU  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as 
follows; 


July  2001 

100% 

August 

100% 

September 

100% 

October 

100% 

November 

100% 

December 

100% 

Fig.  2.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
NORSAR  data  in  our  tape  archive,  on  a  day-by-day  basis  for  the  reporting  period. 

J.  Torstveit 
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NOA  uptime  for  Jul.  2001  (This  month:  100.00  percent) 


Fig.  2.1.1.  The  figure  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability 
of  NOA  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period.  (Page  1  of  2, 
July-Sep  2001 ). 
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NOA  uptime  for  Oct.  2001  (This  month:  100.00  percent) 


NOA  uptime  for  Nov.  2001  (This  month:  100.00  percent) 


NOA  uptime  for  Dec.  2001  (This  month:  100-00  percent) 


Fig.  2.1.1.  (cont.)  (Page  2  of  2,  Oct-Dec  2001). 
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NOA  Event  Detection  Operation 

In  Table  2.1.1  some  monthly  statistics  of  the  Detection  and  Event  Processor  operation  are 
given.  The  table  lists  the  total  number  of  detections  (DPX)  triggered  by  the  on-line  detector,  the 
total  number  of  detections  processed  by  the  automatic  event  processor  (EPX)  and  the  total 
number  of  events  accepted  after  analyst  review  (teleseismic  phases,  core  phases  and  total). 


Total 

DPX 

Total 

EPX 

Accepted  Events 

Sum 

Daily 

P'phases 

Core 

Phases 

Jul  2001 

7,589 

1,077 

359 

62 

421 

13.6 

Aug 

10,045 

1,550 

250 

129 

379 

12.2 

Sep 

8,228 

832 

225 

78 

303 

10.1 

Oct 

10,728 

856 

249 

91 

340 

11.0 

Nov 

12,043 

942 

204 

70 

274 

9.1 

Dec 

11,867 

908 

228 

59 

287 

60,500 

6,165 

1,515 

489 

2,004 

10.9 

Table  2.1.1.  Detection  and  Event  Processor  statistics,  1  July  -  31  December  2001. 


NOA  detections 

The  number  of  detections  (phases)  reported  by  the  NORSAR  detector  during  day  182,  2001, 
through  day  365, 2001,  was  60,500,  giving  an  average  of  329  detections  per  processed  day  (184 
days  processed). 

B.  Paulsen 
U.  Baadshaug 
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2.2  PS28  —  Primary  Seismic  Station  ARCES 

The  average  recording  time  was  99.20%  as  compared  to  99.74%  for  the  previous  period. 
Table  2.2.1  lists  the  reasons  for  and  time  periods  of  the  main  downtimes  in  the  reporting  period. 


Date 

Time 

Cause 

11  Jul 

2241 

- 

Power  failure 

12Jul 

18  Jul 

1835 

-  1836 

Power  failure 

19  Jul 

20  Jul 

1634 

-  0221 

-  2357 

Power  failure 

Table  2.2.1.  The  main  interruptions  in  recording  of  ARCES  data  at  NDPC,  1  July  - 
31  December  2001. 


Monthly  uptimes  for  the  ARCES  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmission  lines,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows: 


July  2001 

95.18% 

August 

100% 

September 

100% 

October 

100% 

November 

100% 

December 

100% 
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ARCES  uptime  for  Oct.  2001  (This  month;  100.00  percent) 


ARCES  uptime  for  Nov.  2001  (This  month:  100.00  percent) 


ARCES  uptime  for  Dec.  2001  (This  month:  100.00  percent) 


Fig.  2.2.1  (cont.)  (Page  2  of  2,  Oct-Dec  2001). 
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Event  Detection  Operation 
ARCES  detections 

The  number  of  detections  (phases)  reported  during  day  182, 2001,  through  day  365, 2001,  was 
145,357,  giving  an  average  of  790  detections  per  processed  day  (184  days  processed). 

Events  automatically  located  by  ARCES 

During  days  182, 2001,  through  365,  2001, 11,215  local  and  regional  events  were  located  by 
ARCES,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average  of 
61.0  events  per  processed  day  (184  days  processed).  44%  of  these  events  are  within  300  km, 
and  76%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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2.3  AS72  —  Auxiliary  Seismic  Station  Spitsbergen 

The  average  recording  time  was  99.99%  as  compared  to  93.52%  for  the  previous  reporting 
period. 

Monthly  uptimes  for  the  Spitsbergen  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows: 


July  2001 

100% 

August 

100% 

September 

100% 

October 

99.97% 

November 

99.99% 

December 

99.97% 

J.  Torstveit 


NORSAR  Sci.  Rep.  1-2002 


February  2002 


13  3  4  5  6  7  6  •  16  11  12  13  t4  15  16  t7  16  16  20  21  22  23  2*  25  26  37  28  »  30  31 
SVA  Uptime  for  Jui.  2001  (This  month:  100.00  percent) 


SVA  uptime  for  Aug  2001  (This  month;  100.00  percent) 


SVA  uptime  for  Sep.  2001  (This  month:  100.00  percent) 


Fig.  23.1.  The  figure  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability 
of  Spitsbergen  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period.  (Page 
1  of  2,  Jul-Sep  2001). 
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SVA  uptime  for  Oct.  2001  (This  month:  99.97  percent) 


SVA  uptime  for  Nov.  2001  (This  month:  99.99  percent) 


I 


SVA  uptime  for  Oec.  2001  (This  month;  99.97  percent) 


Fig.  2.3.1  (cont.)  (Page  2  of  2,  Oct-Dec  2001). 
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Event  Detection  Operation 
Spitsbergen  array  detections 

The  number  of  detections  (phases)  reported  from  day  182, 2001,  through  day  365, 2001,  was 
400,191,  giving  an  average  of  2175  detections  per  processed  day  (184  days  processed). 

Events  automatically  located  by  the  Spitsbergen  array 

During  days  182, 2001,  through  365, 2001, 42,476  local  and  regional  events  were  located  by 
the  Spitsbergen  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an 
average  of  184.0  events  per  processed  day  (184  days  processed).  60%  of  these  events  are 
within  300  km,  and  82%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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2.4  AS73  —  Auxiliary  Seismic  Station  at  Jan  Mayen 

The  IMS  auxiliary  seismic  network  will  include  a  three-component  station  on  the  Norwegian 
island  of  Jan  Mayen.  The  station  location  given  in  the  protocol  to  the  Comprehensive  Nuclear- 
Test-Ban  Treaty  is  70.9°N,  8.7°W. 

The  University  of  Bergen  has  operated  a  seismic  station  at  this  location  since  1970.  An  invest¬ 
ment  in  the  new  station  at  Jan  Mayen  will  be  made  in  due  course  and  in  accordance  with  Prep- 
Com  program  and  budget  decisions.  In  the  meanwhile  data  from  the  existing  seismic  station  on 
Jan  Mayen  are  being  transmitted  to  the  NDC  at  Kjeller  and  to  the  University  of  Bergen  via  a 
VS  AT  link  installed  in  April  2000. 

S.  Mykkeltveit 

2.5  IS37  —  Infrasound  Station  at  Karasjok 

The  IMS  infrasound  network  will  include  a  station  at  Karasjok  in  northern  Norway.  The  coor¬ 
dinates  given  for  this  station  are  69.5°N,  25.5°E.  These  coordinates  coincide  with  those  of  the 
primary  seismic  station  PS28. 

A  site  survey  for  this  station  was  carried  out  during  June/July  1998  as  a  cooperative  effort 
between  the  Provisional  Technical  Secretariat  of  the  CTBTO  and  NORSAR.  Analysis  of  the 
data  collected  at  several  potential  locations  for  this  station  in  and  around  Karasjok  has  been 
completed.  The  results  of  this  analysis  have  lead  to  a  reconunendation  on  the  exact  location  of 
the  infrasound  station.  This  location  needs  to  be  surveyed  in  detail.  The  next  step  will  be  to 
approach  the  local  authorities  to  obtain  the  permissions  needed  to  establish  the  station.  Station 
installation  is  now  expected  to  take  place  in  the  year  2003. 

S.  Mykkeltveit 

2.6  RN49  —  Radionuclide  Station  on  Spitsbergen 

The  IMS  radionuclide  network  will  include  a  station  at  Longyearbyen  on  the  island  of  Spitsber¬ 
gen,  with  location  78.2°N,  16.4°E,  as  given  in  the  protocol  to  the  Comprehensive  Nuclear-Test- 
Ban  Treaty.  These  coordinates  coincide  with  those  of  the  auxiliary  seismic  station  AS72. 
According  to  PrepCom  decision,  this  station  will  also  be  among  those  IMS  radionuclide  sta¬ 
tions  that  will  have  a  capability  of  monitoring  for  the  presence  of  relevant  noble  gases  upon 
entry  into  force  of  the  CTBT. 

A  site  survey  for  this  station  was  carried  out  in  August  of  1999  by  NORSAR,  in  cooperation 
with  the  Norwegian  Radiation  Protection  Authority.  The  site  survey  report  to  the  PTS  con¬ 
tained  a  recommendation  to  establish  this  station  at  Plataberget,  some  20  km  away  from  the 
Treaty  location.  The  PrepCom  approved  the  corresponding  coordinate  change  in  its  meeting  in 
May  2000.  The  station  installation  was  part  of  PrepCom’s  work  program  and  budget  for  the 
year  2000.  The  infrastructure  for  housing  the  station  equipment  has  been  established,  and  a 
noble  gas  detection  system,  based  on  the  Swedish  “SAUNA”  design,  was  installed  at  this  site  in 
May  2001 .  A  particulate  station  (“ARAME”  design)  was  installed  at  the  same  location  in  Sep¬ 
tember  2001.  Currently,  the  two  systems  are  undergoing  testing  and  evaluation. 

S.  Mykkeltveit 
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3  Contributing  Regional  Seismic  Arrays 
3.1  NORES 

Average  recording  time  was  70.32%  as  compared  to  99.98%  for  the  previous  period. 

The  outage  between  16  August  and  9  October  was  due  to  a  thunderstorm  that  caused  major 
damage  to  the  equipment. 

Monthly  uptimes  for  the  NORES  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as 
follows: 


July  2001 

100% 

August 

50.70% 

September 

0% 

October 

72.55% 

November 

98.65% 

December 

100% 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
NORES  data  in  our  tape  archive  on  a  day-by-day  basis  for  the  reporting  period. 


J.  Torstveit 
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MORES  Uptime  for  Aug.  2001  (This  month:  50.70  percent) 
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Fig.  3.1.1.  The  figure  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability 
ofNORES  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period  (Page  1  of 
2.  Jul-Sep  2001). 
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NORES  uptime  for  Oct.  2001  (Thie  month;  72.55  percent) 


HIM 

11 

1 

MM 

11 

11 

1 

II 

HIM 

11 

1 

llllll 

11 

11 

MM 

II 

HIM 

11 

1 

llllll 

11 

11 

MM 

II 

HIM 

11 

1 

llllll 

II 

11 

MM 

II 

HIM 

11 

1 

llllll 

II 

11 

MM 

II 

HIM 

II 

1 

llllll 

11 

11 

MM 

II 

HIM 

11 

1 

llllll 

11 

11 

MM 

II 

HIM 

11 

1 

llllll 

11 

11 

MM 

II 

HIM 

11 

1 

llllll 

II 

11 

MM 

II 

HIM 

■■(■EkKk! 

11 

Hn 

1 

IfH 

llllll 

IWMLll.ftkJ 

11 

11 

IHtlH 

MM 

II 

MH 

NORES  uptime  for  Dec.  2001  (This  month:  100.00  percent) 


Fig.  3.1.1  (cont.)  (Page  2  of  2,  Oct-Dec  2001). 
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NORES  Event  Detection  Operation 
NORES  detections 

The  number  of  detections  (phases)  reported  from  day  182, 2001,  through  day  365, 2001,  was 
95,065,  giving  an  average  of  726  detections  per  processed  day  (131  days  processed;  the  pro¬ 
cessing  of  NORES  data  was  suspended  for  a  toM  of  53  days,  due  to  severe  clock  problems  fol¬ 
lowing  a  thunderstorm). 

Events  automatically  located  by  NORES 

During  days  182, 2001,  through  365, 2001, 2351  local  and  regional  events  were  located  by 
NORES,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average  of 
17.9  events  per  processed  day  (131  days  processed).  43%  of  these  events  are  within  300  km, 
and  75%  of  these  events  are  within  1(X)0  km. 

U.  Baadshaug 
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3.2  Hagfors  (IMS  Station  ASlOl) 

The  average  recording  time  was  99.91%  as  compared  to  99.90%  for  the  previous  reporting 
period. 

Table  3.2.1  lists  the  reasons  for  and  times  of  the  main  outages  in  the  reporting  period. 

Date  Time  Cause 

14  Aug  0832  -  1050  Installation  of  new  station 

Table  3.2.1.  The  main  interruptions  in  Hagfors  recordings  at  the  Norwegian 
NDC,  1  July  -  31  December  2001. 


The  Hagfors  array  was  refurbished  in  August  through  the  installation  of  entirely  new  equip¬ 
ment  deployed  at  new  site  locations.  The  “old”  array  continues  its  operation  in  parallel  with  the 
“new”  array.  All  data  reported  here  refer  to  the  “old”  Hagfors  array. 

Monthly  uptimes  for  the  Hagfors  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as 
follows: 


July  2001 

99.99% 

August 

99.55% 

September 

100% 

October 

100% 

November 

99.94% 

December 

99.98% 

J.  Torstveit 
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HFS  uptime  for  Jul.  2001  (This  month;  99.99  percent) 


HFS  uptime  for  Aug.  2001  (This  month;  99.55  percent) 


Fig.  3.2.1.  The  figure  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability 
of  Hagfors  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period  (Page  1 
of  2,  Jul-Sep  2001). 
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Hrs  uptime  for  Oct.  2001  (This  month:  100.00  percent) 


HFS  uptime  for  Nov.  2001  (This  month:  99.94  percent) 
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Hrs  uptime  for  Dec.  2001  (This  month:  99.98  percent) 


Fig.  3.2.1  (cont.)  (Page  2  of  2,  Oct-Dec  2(X)1). 
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Hagfors  Event  Detection  Operation 
Hagfors  array  detections 

The  number  of  detections  (phases)  reported  from  day  182, 2001,  through  day  365, 2001,  was 
103,399,  giving  an  average  of  562  detections  per  processed  day  (184  days  processed). 

Events  automatically  located  by  the  Hagfors  array 

During  days  182, 2001,  through  365, 2001, 2952  local  and  regional  events  were  located  by  the 
Hagfors  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  16.0  events  per  processed  day  (184  days  processed).  56%  of  these  events  are  within  300  km, 
and  83%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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3.3  FINES  (IMS  station  PS17) 

The  average  recording  time  was  99.99%  as  compared  to  99.51%  for  the  previous  reporting 
period. 

Monthly  uptimes  for  the  FINES  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows: 


July  2001 

99.99% 

August 

100% 

September 

100% 

October 

100% 

November 

99.96% 

December 

100% 

Fig.  3.3.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
FINES  data  in  our  tape  archive  on  a  day-by-day  basis  for  the  reporting  period. 

J.  Torstveit 
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PINES  uptime  for  Jul.  2001  (This  month:  99.99  percent) 


FINES  uptime  for  Aug.  2001  (This  month;  100.00  percent) 
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8. 

R. 


FINES  uptime  for  Sep.  2001  (This  month:  100.00  percent) 


Fig.  3.3.1.  The  figure  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability 
of  FINES  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period  (Page  1  of 
2,Jul-Sep2001). 
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FINES  Event  Detection  Operation 
FINES  detections 

The  number  of  detections  (phases)  reported  during  day  182, 2001,  through  day  365,  2001,  was 
41,957,  giving  an  average  of  228  detections  per  processed  day  (184  days  processed). 

Events  automatically  located  by  FINES 

During  days  182, 2001,  through  365, 2001, 2381  local  and  regional  events  were  located  by 
FINES,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average  of  12.9 
events  per  processed  day  (184  days  processed).  79%  of  these  events  are  within  300  km,  and 
88%  of  these  events  are  within  1000  km. 

U.  Baadshaug 


29 


NORSAR  Sci.  Rep.  1-2002 


February  2002 


3.4  Apatity 

The  average  recording  time  was  98.3 1  %  in  the  reporting  period  compared  to  99.92%  during  the 
previous  period. 

Table  3.4.1  lists  the  reasons  for  and  times  of  the  main  outages  during  the  reporting  period  . 


Date 

Time 

Cause 

02Jul 

1528  - 

1724 

Stop  in  Apatity 

01  Oct 

0723  - 

0958 

Stop  in  Apatity 

01  Oct 

1153  - 

Stop  in  Apatity 

02  Oct 

- 

0414 

02  Oct 

1120  - 

Stop  in  Apatity 

03  Oct 

- 

1347 

05  Nov 

1207  - 

1458 

Stop  in  Apatity 

28  Nov 

1326  - 

Stop  in  Apatity 

29  Nov 

- 

0848 

29  Nov 

0925  - 

1008 

Stop  in  Apatity 

29  Nov 

1245  - 

1300 

Stop  in  Apatity 

Table  3.4.1.  The  main  interruptions  in  Apatity  recordings  at  the  NDC,  1  July  - 
31  December  2001. 


Monthly  uptimes  for  the  Apatity  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows: 


July  2001 

99.74% 

August 

100% 

September 

100% 

October 

94.06% 

November 

96.09% 

December 

99.95% 

J.  Torstveit 
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APA  uptime  for  Jul.  2001  (This  month:  99.74  percent) 


Fig.  3.4.1.  The  figure  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability 
of  Apatity  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period  ( Page  1  of 
2,  Jul-Sep  2001 ). 
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Fig.  3.4.1  (cont.)  (Page  2  of  2,  Oct- Dec  2001 ) 
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Apatity  Event  Detection  Operation 
Apatity  array  detections 

The  number  of  detections  (phases)  reported  from  day  182, 2001,  through  day  365, 2001,  was 
229,833,  giving  an  average  of  1200  detections  per  processed  day  (184  days  processed). 

As  described  in  earlier  reports,  the  data  from  the  Apatity  array  are  transferred  by  one-way  (sim¬ 
plex)  radio  links  to  Apatity  city.  The  transmission  suffers  from  radio  disturbances  that  occa¬ 
sionally  result  in  a  large  number  of  small  data  gaps  and  spikes  in  the  data.  In  order  for  the 
conununication  protocol  to  correct  such  errors  by  requesting  retransmission  of  data,  a  two-way 
radio  link  would  be  needed  (duplex  radio).  However,  it  should  be  noted  that  noise  from  cultural 
activities  and  from  the  nearby  lakes  cause  most  of  the  unwanted  detections.  These  unwanted 
detections  are  “filtered”  in  the  signal  processing,  as  they  give  seismic  velocities  that  are  outside 
accepted  limits  for  regional  and  teleseismic  phase  velocities. 

Events  automatically  located  by  the  Apatity  array 

During  days  182, 2001,  through  365, 2001,  2246  local  and  regional  events  were  located  by  the 
Apatity  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  12.2  events  per  processed  day  (184  days  processed).  42%  of  these  events  are  within  300  km, 
and  73%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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3.5  Regional  Monitoring  System  Operation  and  Analysis 

The  Regional  Monitoring  System  (RMS)  was  installed  at  NORSAR  in  December  1989  and 
was  operated  at  NORSAR  from  1  January  1990  for  automatic  processing  of  data  from  ARCES 
and  NORES.  A  second  version  of  RMS  that  accepts  data  from  an  arbitrary  number  of  arrays 
and  single  3-component  stations  was  installed  at  NORSAR  in  October  1991,  and  regular  oper¬ 
ation  of  the  system  comprising  analysis  of  data  from  the  4  arrays  ARCES,  NORES,  FINES  and 
GERES  started  on  15  October  1991.  As  opposed  to  the  first  version  of  RMS,  the  one  in  current 
operation  also  has  the  capability  of  locating  events  at  teleseismic  distance. 

Data  from  the  Apatity  array  were  included  on  14  December  1992,  and  from  the  Spitsbergen 
array  on  12  January  1994.  Detections  from  the  Hagfors  array  were  available  to  the  analysts  and 
could  be  added  manually  during  analysis  from  6  December  1994.  After  2  February  1995,  Hag¬ 
fors  detections  were  also  used  in  the  automatic  phase  association. 

Since  24  April  1999,  RMS  has  processed  data  from  all  the  seven  regional  arrays  ARCES, 
NORES,  FINES,  GERES  (until  January  2000),  Apatity,  Spitsbergen,  and  Hagfors.  Starting 
19  September  1999,  waveforms  and  detections  from  the  NORSAR  array  have  also  been  avail¬ 
able  to  the  analyst. 

Phase  and  event  statistics 

Table  3.5.1  gives  a  summary  of  phase  detections  and  events  declared  by  RMS.  From  top  to  bot¬ 
tom  the  table  gives  the  total  number  of  detections  by  the  RMS,  the  number  of  detections  that 
are  associated  with  events  automatically  declared  by  the  RMS,  the  number  of  detections  that 
are  not  associated  with  any  events,  the  number  of  events  automatically  declared  by  the  RMS, 
and  finally  the  total  number  of  events  worked  on  interactively  (in  accordance  with  criteria  that 
vary  over  time;  see  below)  and  defined  by  the  analyst. 

New  criteria  for  interactive  event  analysis  were  introduced  from  1  January  1994.  Since  that 
date,  only  regional  events  in  areas  of  special  interest  (e.g,  Spitsbergen,  since  it  is  necessary  to 
acquire  new  knowledge  in  this  region)  or  other  significant  events  (e.g,  felt  earthquakes  and 
large  industrial  explosions)  were  thoroughly  analyzed.  Teleseismic  events  of  special  interest 
are  also  analyzed. 

To  further  reduce  the  workload  on  the  analysts  and  to  focus  on  regional  events  in  preparation 
for  Ganuna-data  submission  during  GSETT-3,  a  new  processing  scheme  was  introduced  on  2 
February  1995.  The  GBF  (Generalized  Beamforming)  program  is  used  as  a  pre-processor  to 
RMS,  and  only  phases  associated  to  selected  events  in  northern  Europe  are  considered  in  the 
automatic  RMS  phase  association.  All  detections,  however,  are  still  available  to  the  analysts 
and  can  be  added  manually  during  analysis. 
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Jul 

01 

Aug 

01 

Sep 

01 

Oct 

01 

Nov 

01 

Dec 

01 

Total 

Phase  detections 

194113 

180434 

163008 

165445 

135244 

159879 

998123 

-  Associated  phases 

7043 

6003 

4657 

5588 

3252 

5029 

31572 

-  Unassociated  phases 

187070 

174431 

158351 

159857 

131992 

154850 

966551 

Events  automatically 
declared  by  RMS 

1718 

1329 

1050 

1270 

561 

1127 

7055 

No.  of  events  defined  by 
the  analyst 

80 

104 

103 

62 

98 

80 

527 

Table  3.5.1,  RMS  phase  detections  and  event  summary. 

U.  Baadshaug 
B.  Paulsen 
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4  NDC  and  Field  Activities 

4.1  NDC  Activitities 

NORSAR  is  functioning  as  the  Norwegian  National  Data  Center  (NDC)  for  CTBT  verification. 
Six  monitoring  stations,  comprising  altogether  119  field  instruments,  will  be  located  on  Nor¬ 
wegian  territory  as  part  of  the  future  IMS  as  described  elsewhere  in  this  report.  The  four  seis¬ 
mic  IMS  stations  are  all  in  operation  today,  and  three  of  them  are  currently  providing  data  to 
the  IDC.  The  radionuclide  station  at  Spitsbergen  is  currently  operating  in  a  testing  mode, 
whereas  the  infrasound  station  in  northern  Norway  will  need  to  be  established  within  the  next 
few  years.  Data  recorded  by  the  Norwegian  stations  is  being  transmitted  in  real  time  to  the  Nor¬ 
wegian  NDC,  and  provided  to  the  IDC  through  the  Global  Communications  Infrastructure 
(GCI).  Norway  is  connected  to  the  GCI  with  a  frame  relay  link  to  Vienna. 

Operating  the  Norwegian  IMS  stations  will  require  increased  resources  and  additional  person¬ 
nel  both  at  the  NDC  and  in  the  field.  It  will  require  establishing  new  and  strictly  defined  proce¬ 
dures  as  well  as  increased  emphasis  on  regularity  of  data  recording  and  timely  data 
transmission  to  the  IDC  in  Vieima.  Anticipating  these  requirements,  a  new  organizational  unit 
has  been  established  at  NORSAR  to  form  a  core  group  for  the  future  Norwegian  NDC  for 
treaty  monitoring.  The  NDC  will  carry  out  all  the  technical  tasks  required  in  support  of  Nor¬ 
way’s  treaty  obligations.  NORSAR  will  also  carry  out  assessments  of  events  of  special  interest, 
and  advise  the  Norwegian  authorities  in  technical  matters  relating  to  treaty  compliance. 

Verification  functions;  information  received  from  the  IDC 

After  the  CTBT  enters  into  force,  the  IDC  will  provide  data  for  a  large  number  of  events  each 
day,  but  will  not  assess  whether  any  of  them  are  likely  to  be  nuclear  explosions.  Such  assess¬ 
ments  will  be  the  task  of  the  States  Parties,  and  it  is  important  to  develop  the  necessary  national 
expertise  in  the  participating  countries.  An  important  task  for  the  Norwegian  NDC  will  thus  be 
to  make  independent  assessments  of  events  of  particular  interest  to  Norway,  and  to  communi¬ 
cate  the  results  of  these  analyses  to  the  Norwegian  Ministry  of  Foreign  Affairs. 

Monitoring  the  Arctic  region 

Norway  will  have  monitoring  stations  of  key  importance  for  covering  the  Arctic,  including 
Novaya  Zemlya,  and  Norwegian  experts  have  a  unique  competence  in  assessing  events  in  this 
region.  On  several  occasions  in  the  past,  seismic  events  near  Novaya  Zemlya  have  caused  polit¬ 
ical  concern,  and  NORSAR  specialists  have  contributed  to  clarifying  these  issues. 

International  cooperation 

After  entry  into  force  of  the  treaty,  a  number  of  countries  are  expected  to  establish  national 
expertise  to  contribute  to  the  treaty  verification  on  a  global  basis.  Norwegian  experts  have  been 
in  contact  with  experts  from  several  countries  with  the  aim  to  establish  bilateral  or  multilateral 
cooperation  in  this  field.  One  interesting  possibility  for  the  future  is  to  establish  NORSAR  as  a 
regional  center  for  European  cooperation  in  the  CTBT  verification  activities. 
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NORSAR  event  processing 

The  automatic  routine  processing  of  NORSAR  events  as  described  in  NORSAR  Sci.  Rep.  No. 
2-93/94,  has  been  running  satisfactorily.  The  analyst  tools  for  reviewing  and  updating  the  solu¬ 
tions  have  been  continuously  modified  to  simplify  operations  and  improve  results.  NORSAR  is 
currently  applying  teleseismic  detection  and  event  processing  using  the  large-aperture 
NORSAR  array  as  well  as  regional  monitoring  using  the  network  of  small-aperture  arrays  in 
Fennoscandia  and  adjacent  areas. 

Certification  ofPS28 

On  8  November  2001  the  IMS  station  PS28-ARCES  was  formally  certified.  PTS  personnel  vis¬ 
ited  the  station  during  the  winter  of  2000,  and  carried  out  a  detailed  technical  evaluation.  As  a 
result  of  this  inspection  and  subsequent  discussions  between  NORSAR  and  the  PTS,  it  was 
concluded  that  PS28  needed  only  one  enhancement  in  order  to  be  certified:  to  install  a  central¬ 
ized  authentication  process  at  the  central  array  recording  facility.  After  this  was  done  during 
the  fall  of  2001  and  subsequently  verified  by  the  PTS,  station  certification  was  granted. 

Communication  topology 

Norway  has  elected  to  use  the  option  for  an  independent  subnetwork,  which  will  connect  the 
IMS  stations  AS72,  AS73,  PS28,  IS37  and  RN49  operated  by  NORSAR  to  the  GCI  at 
NOR_NDC.  A  contract  has  been  concluded  and  VSAT  antennas  have  been  installed  at  each 
station  in  the  network.  Under  the  same  contract,  VSAT  antennas  for  6  of  the  PS27  subarrays 
have  been  installed  for  intra-array  communication.  The  seventh  subarray  is  connected  to  the 
central  recording  facihty  via  a  leased  land  line.  The  central  recording  facility  for  PS27  is  con¬ 
nected  directly  to  the  GCI  (Basic  Topology).  All  the  VSAT  communication  is  functioning  satis¬ 
factorily. 

The  Norwegian  NDC  has  been  cooperating  with  institutions  in  other  countries  for  transmission 
of  IMS  data  to  the  Prototype  EDC  during  GSETT-3.  Details  on  this  can  be  found  in  Section  4.2. 

Jan  Fyen 


4.2  Status  Report:  Norway’s  Participation  in  GSETT-3 
Introduction 

This  contribution  is  a  report  for  the  period  July  -  December  2001  on  activities  associated  with 
Norway’s  participation  in  the  GSETT-3  experiment,  which  provides  data  to  the  International 
Data  Centre  (IDC)  in  Vienna  on  an  experimental  basis  until  the  participating  stations  have  been 
conunissioned  as  part  of  the  International  Monitoring  System  (IMS)  network  defined  in  the 
protocol  to  the  Comprehensive  Nuclear-Test-Ban  Treaty.  This  report  represents  an  update  of 
contributions  that  can  be  found  in  previous  editions  of  NORSAR’s  Semiannual  Technical 
Summary.  It  is  noted  that  as  of  31  December  2001,  two  out  of  the  three  Norwegian  seismic  sta¬ 
tions  providing  data  to  the  IDC  have  been  formally  certified  and  thus  commissioned  as  part  of 
the  IMS  network. 
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Norwegian  GSETT-3  stations  and  communications  arrangements 

During  the  reporting  interval  1  July  -  31  December  2001,  Norway  has  provided  data  to  the 
GSETT-3  experiment  from  the  three  seismic  stations  shown  in  Fig.  4.2.1.  The  NORSAR  array 
(PS27,  station  code  NOA)  is  a  60  km  aperture  teleseismic  array,  comprised  of  7  subarrays,  each 
containing  sbc  vertical  short  period  sensors  and  a  three-component  broadband  instrument. 
ARCES  is  a  25-element  regional  array  with  an  aperture  of  3  km,  whereas  the  Spitsbergen  array 
(station  code  SPITS)  has  9  elements  within  a  1-km  aperture.  ARCES  and  SPITS  both  have  a 
broadband  three-component  seismometer  at  the  array  center. 

The  intra-array  communication  for  NOA  has  been  achieved  by  a  land  line  for  subarray  NC6 
and  VSAT  hnks  based  on  TDMA  technology  for  the  other  6  subarrays.  The  central  recording 
facihty  of  NOA  is  at  NOR_NDC. 

Continuous  ARCES  data  have  been  transmitted  from  the  ARCES  site  to  NOR_NDC  using  a 
64  kbits/s  VSAT  satellite  link,  based  on  BOD  technology. 

Continuous  SPITS  data  have  been  transmitted  to  NOR_NDC  via  a  VSAT  terminal  located  at 
Platiberget  in  Longyearbyen  (which  is  the  site  of  the  IMS  radionuclide  monitoring  station 
RN49  installed  during  2(X)1). 

Seven-day  station  buffers  have  been  established  at  the  ARCES  and  SPITS  sites  and  at  all  NOA 
subairay  sites,  as  well  as  at  NOR_NDC  for  ARCES,  SPITS  and  NOA  (central  array  station 
buffer). 

The  NOA  and  ARCES  arrays  are  primary  stations  in  the  GSETT-3  network  and  the  IMS,  which 
implies  that  data  from  these  stations  are  transmitted  continuously  to  the  receiving  international 
data  center.  Since  October  1999,  these  data  have  been  transmitted  (from  NOR_NDC)  via  the 
Global  Communications  Infrastructure  (GCI)  to  the  IDC  in  Vienna,  whereas  transmission  of 
the  same  data  to  the  Prototype  International  Data  Center  (PIDC)  in  Arlington,  VA,  was  discon¬ 
tinued  on  7  February  2000.  The  SPITS  array  is  an  auxiliary  station  in  GSETT-3  and  the  IMS, 
and  the  SPITS  data  have  been  available  to  both  the  IDC  and  the  PIDC  throughout  the  reporting 
period  on  a  request  basis  via  use  of  the  AutoDRM  protocol  (Kradolfer,  1993;  Kradolfer,  1996). 
The  Norwegian  stations  are  thus  participating  in  GSETT-3  with  the  same  status  (primary/auxil¬ 
iary  seismic  stations)  they  have  in  the  IMS  defined  in  the  protocol  to  the  Comprehensive 
Nuclear-Test-Ban  Treaty.  In  addition,  continuous  data  from  all  three  arrays  are  being  transmit¬ 
ted  to  the  US  NDC. 

Uptimes  and  data  availability 

Figs.  4.2.2  -  4.2.3  show  the  monthly  uptimes  for  the  Norwegian  GSETT-3  primary  stations 
ARCES  and  NOA,  respectively,  for  the  period  1  July  -  31  December  2001,  given  as  the  hatched 
(taller)  bars  in  these  figures.  These  barplots  reflect  the  percentage  of  the  waveform  data  that  are 
available  in  the  NOR_NDC  tape  archives  for  these  two  arrays.  The  downtimes  inferred  from 
these  figures  thus  represent  the  cumulative  effect  of  field  equipment  outages,  station  site  to 
NOR_NDC  communication  outage,  and  NOR_NDC  data  acquisition  outages. 

Figs.  4.2.2-4.2.3  also  give  the  data  availability  for  these  two  stations  as  reported  by  the  PIDC  in 
the  PIDC  Station  Status  reports.  The  main  reason  for  the  discrepancies  between  the 
NOR_NDC  and  PIDC  data  availabilities  as  observed  from  these  figures  is  the  difference  in  the 
ways  the  two  data  centers  report  data  availability  for  arrays:  Whereas  NOR_NDC  reports  an 
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array  station  to  be  up  and  available  if  at  least  one  channel  produces  useful  data,  the  PEDC  uses 
weights  where  the  reported  availability  (capability)  is  based  on  the  number  of  actually  operat¬ 
ing  channels.  The  PIDC  receives  its  ARCES  and  NOA  data  via  the  IDC  in  Vienna. 

Use  of  the  AutoDRM  protocol 

NOR_NDC’s  AutoDRM  has  been  operational  since  November  1995  (Mykkeltveit  &  Baads- 
haug,  1996).  The  PIDC  started  actively  and  routinely  using  NOR_NDC’s  AutoDRM  service 
after  SPITS  changed  its  station  status  from  primary  to  auxiliary  on  1  October  1996.  The 
monthly  number  of  requests  by  the  PIDC  for  SPITS  data  for  the  period  July  -  December  2001 
is  shown  in  Fig.  4.2.4. 

NDC  automatic  processing  and  data  analysis 

These  tasks  have  proceeded  in  accordance  with  the  descriptions  given  in  Mykkeltveit  and 
Baadshaug  (1996).  For  the  period  July  -  December  2001,  NOR_NDC  derived  information  on 
524  supplementary  events  in  northern  Europe  and  submitted  this  information  to  the  Finnish 
NDC  as  the  NOR_NDC  contribution  to  the  joint  Nordic  Supplementary  (Ganuna)  Bulletin, 
which  in  turn  is  forwarded  to  the  PIDC.  These  events  are  plotted  in  Fig.  4.2.5. 

Data  forwarding  for  GSETT-3  stations  in  other  countries 

NOR_NDC  continued  to  provide  conununications  for  the  GSETT-3  auxiliary  station  at  Nilore, 
Pakistan,  through  a  VSAT  satellite  link  between  NOR_NDC  and  Pakistan’s  NDC  in  Nilore. 
The  PIDC  as  well  as  the  IDC  obtain  data  from  the  Hagfors  array  (HFS)  in  Sweden  through 
requests  to  the  AutoDRM  server  at  NOR_NDC  (in  the  same  way  requests  for  Spitsbergen  array 
data  are  handled,  see  above).  Fig.  4.2.6  shows  the  monthly  number  of  requests  for  HFS  data 
from  the  two  PIDC  accounts  “pipeline”  and  “testbed”. 

Current  developments  and  future  plans 

NOR_NDC  is  continuing  the  efforts  towards  improving  and  hardening  all  critical  data  acquisi¬ 
tion  and  data  forwarding  hardware  and  software  components,  so  as  to  meet  future  requirements 
related  to  operation  of  IMS  stations  to  the  maximum  extent  possible. 

The  PrepCom  has  tasked  its  Working  Group  B  with  overseeing,  coordinating,  and  evaluating 
the  GSETT-3  experiment.  The  PrepCom  has  also  encouraged  states  that  operate  IMS- 
designated  stations  to  continue  to  do  so  on  a  voluntary  basis  and  in  the  framework  of  the 
GSETT-experiment  until  such  time  that  the  stations  have  been  certified  for  formal  inclusion  in 
IMS.  The  NOA  array  was  formally  certified  by  the  PTS  on  28  July  2000,  and  a  contract  with 
the  PTS  in  Vienna  currently  provides  partial  funding  for  operation  and  maintenance  of  this  sta¬ 
tion.  The  ARCES  array  was  formally  certified  by  the  PTS  on  8  November  2001.  It  is  expected 
that  a  contract  will  be  signed  with  the  PTS  for  operation  and  maintenance  of  this  station,  with 
an  effective  date  of  1  January  2002.  Provided  that  adequate  funding  continues  to  be  made  avail¬ 
able  (from  the  PTS  and  the  Norwegian  Ministry  of  Foreign  Affairs),  we  envisage  continuing 
the  provision  of  data  from  all  Norwegian  IMS-designated  seismic  stations  without  intermption 
to  the  IDC  in  Vienna. 
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Fig.  4.2.1.  The  figure  shows  the  locations  and  configurations  of  the  three  Norwegian  seismic  array 
stations  that  have  provided  data  to  the  GSETT-3  experiment  during  the  period  1  July  - 
31  December  2001.  The  data  from  these  stations  are  transmitted  continuously  and  in  real  time 
to  the  Norwegian  NDC  (NOR_NDC).  The  stations  NOA  and  ARCES  have  participated  in 
GSETT-3  as  primary  stations,  whereas  SPITS  has  contributed  as  an  auxiliary  station. 
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ARCES  data  availability  at  NDC  and  PIDC 


Jul  Aug  Sep  Oel  Nov  Dec 


Fig.  4.2.2.  The  figure  shows  the  monthly  availability  of  ARCES  array  data  for  the  period  June  - 

December  2001  at  NOR_NDC  and  the  PIDC.  See  the  text  for  explanation  of  differences  in  def¬ 
inition  of  the  term  “data  availability”  between  the  two  centers.  The  higher  values  (hatched 
bars)  represent  the  NOR_NDC  data  availability. 


NOA  data  availability  at  NDC  and  PIDC 


Fig.  4.2.3.  The  figure  shows  the  monthly  availability  of  NORSAR  array  data  for  the  period  June  - 
December  2001  at  NOR_NDC  and  the  PIDC.  See  the  text  for  explanation  of  differences  in  def¬ 
inition  of  the  term  “data  availability”  between  the  two  centers.  The  higher  values  (hatched 
bars)  represent  the  NORJNDC  data  availability. 
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Fig.  4.2.4.  The  figure  shows  the  monthly  number  of  requests  received  by  NOR_NDC  from  the  PIDC 
for  SPITS  waveform  segments  during  July  -  December  2001. 
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Reviewed  Supplementary  events 


Fig.  4.2.5.  The  map  shows  the  524  events  in  and  around  Norway  contributed  by  NOR_NDC  during 
July  -  December  2001  as  supplementary  (Gamma)  events  to  the  PIDC,  as  part  of  the  Nordic 
supplementary  data  compiled  by  the  Finnish  NDC.  The  map  also  shows  the  seismic  stations 
used  in  the  data  analysis  to  define  these  events. 
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AutoDRM  HFS  requests  received  by  NOR_NDC  from  pipeline  and  testbed 


Fig.  4.2.6.  The  figure  shows  the  monthly  number  of  requests  received  by  NOR_NDCfrom  the  PIDC 
for  HFS  waveform  segments  during  July  -  December  2001. 
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4.3  Field  Activities 

The  activities  at  the  NORSAR  Maintenance  Center  (NMC)  at  Hamar  currently  includes  work 
related  to  operation  and  maintenance  of  the  following  IMS  seismic  stations:  the  NOA  teleseis- 
mic  array  (PS27),  the  ARCES  array  (PS28)  and  the  Spitsbergen  array  (AS72).  Some  prepara¬ 
tory  work  has  also  been  carried  out  in  connection  with  the  seismic  station  on  Jan  Mayen 
(AS73),  the  infrasound  station  at  Karasjok  (IS37)  and  the  radionuclide  station  at  Spitsbergen 
(RN49).  NORSAR  also  acts  as  a  consultant  for  the  operation  and  maintenance  of  the  Hagfors 
array  in  Sweden  (AS  101). 

In  addition  to  the  above  activities,  which  are  directly  related  to  the  International  Monitoring 
System,  NORSAR ’s  field  staff  are  continuing,  within  available  resources,  to  maintain  the 
small-aperture  NORES  array,  which  is  co-located  with  NOA  subarray  06C.  These  efforts  are 
given  low  priority,  since  there  is  no  requirement  for  specific  uptimes  at  NORES. 

NORSAR  carries  out  the  field  activities  relating  to  IMS  stations  in  a  manner  generally  consis¬ 
tent  with  the  requirements  specified  in  the  appropriate  IMS  Operational  Manuals,  which  are 
currently  being  developed  by  Working  Group  B  of  the  Preparatory  Commission.  For  seismic 
stations  these  specifications  are  contained  in  the  Operational  Manual  for  Seismological  Moni¬ 
toring  and  the  International  Exchange  of  Seismological  Data  (CTBTAVGB/TL-l  1/2),  currently 
available  in  a  draft  version. 

All  regular  maintenance  on  the  NORSAR  field  systems  is  conducted  on  a  one-shift-per-day, 
five-day-per-week  basis.  The  maintenance  tasks  include: 

•  Operating  and  maintaining  the  seismic  sensors  and  the  accociated  digitizers,  authentication 
devices  and  other  electronics  components. 

•  Maintaining  the  power  supply  to  the  field  sites  as  well  as  backup  power  supplies  . 

•  Operating  and  maintaining  the  VSATs,  the  data  acquisition  systems  and  the  intra-array  data 
transmission  systems. 

•  Assisting  the  NDC  in  evaluating  the  data  quality  and  making  the  necessary  changes  in  gain 
settings,  frequency  response  and  other  operating  characteristics  as  required. 

•  Carrying  out  preventive,  routine  and  emergency  maintenance  to  ensure  that  all  field  sys¬ 
tems  operate  properly. 

•  Maintaining  a  computerized  record  of  the  utilization,  status,  and  maintenance  history  of  all 
site  equipment. 

•  Providing  appropriate  security  measures  to  protect  against  incidents  such  as  intrusion,  theft 
and  vandalism  at  the  field  installations. 

Details  of  the  daily  maintenance  activities  are  kept  locally.  As  part  of  its  contract  with  CTBTO/ 
PTS,  NORSAR  submits,  when  applicable  problem  reports,  outage  notification  reports  and 
equipment  status  reports.  Die  contents  of  these  reports,  and  the  circumstances  under  which 
they  will  be  submitted,  are  specified  in  the  draft  Operational  Manual. 

P.W.  Larsen 
K.A.  Lpken 
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6  Summary  of  Technical  Reports  /  Papers  Published 

6.1  Estimating  global  and  regional  IMS  detection  capability 

Introduction 

The  primary  seismic  network  of  the  International  Monitoring  System  (IMS)  for  verifying  com¬ 
pliance  with  the  Comprehensive  Nuclear  Test  Ban  Treaty  (CTBT)  consists  of  49  stations  (see 
Fig.  6.1.1),  out  of  which  35  are  installed  and  operational  as  of  July  2001.  These  stations  are  the 
key  element  of  the  IMS  as  they  are  used  for  detecting  events  that  might  be  violations  of  the 
CTBT.  We  have  in  this  study  used  the  threshold  monitoring  (TM)  method  (Kvaema  and 
Ringdal,  1999,  Ringdal  and  Kvaema,  1989, 1992)  to  assess  the  detection  capability  of  the  IMS  . 
seismic  network.  The  TM  method  is  capable  of  using  actual  seismic  data  for  a  given  time  inter¬ 
val  as  the  basis  for  the  detectability  calculations.  In  cases  when  a  seismic  station  did  not  pro¬ 
vide  data  during  the  time  period  under  study,  an  estimated  background  noise  level  can  be 
assigned.  These  noise  estimates  can  be  based  on  results  from  earlier  studies,  or  they  can  be 
taken  from  stations  assumed  to  have  similar  noise  characteristics. 


Fig.  6.1.1.  Station  configuration  of  the  full  IMS  primary  seismic  network. 
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The  detection  capability  of  the  full  IMS  primary  seismic  network 

We  have  used  the  time  interval  10:20-1 1:00  on  29  June  2001  as  the  basis  for  estimating  the 
detection  capability  of  the  full  IMS  primary  seismic  network.  This  time  interval  does  not  con¬ 
tain  any  major  seismic  events  and  most  of  the  stations  have  relatively  quiet  noise  conditions. 
Fig.  6.1.2  shows  the  short-term-averages  (STAs)  representing  the  noise  levels  at  each  of  the 
IMS  primary  seismic  stations  operational  as  of  29  June  2001.  For  stations  not  providing  data 
during  this  time  interval,  a  typical  constant  noise  level  has  been  assigned.  These  stations  have 
the  label  const  at  the  top  of  each  panel. 

The  noise  levels  assigned  for  the  planned  primary  seismic  stations  are  shown  in  Fig.  6.1.3. 
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Fig.  6.1.2.  Real  or  estimated  noise  levels  of  current  IMS  primary  seismic  stations  for  the  time 
period  10:20-11:00  on  29  June  2001. 
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Fig.  6.1.3.  Noise  levels  assigned  for  the  planned  IMS  primary  seismic  stations. 


We  have  used  the  average  noise  levels  shown  in  Figs.  6.1.2  and  6.1.3  to  calculate  the  three-sta¬ 
tion  network  detection  capability  at  the  90%  confidence  level  of  the  full  IMS  primary  seismic 
network.  The  results  are  shown  in  Fig.  6.1.4.  Very  good  detection  capability  is  found  in  north¬ 
ern  Europe  and  North  America  where  several  high  performance  seismic  arrays  are  installed. 
The  results  shown  in  Fig.  6.1.4  agree  well  with  the  statistical  simulations  provided  by  the 
NetSim  program  (Sereno  et.  al,  1990). 
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Fig.  6.1.4.  Three  station  detection  capability  of  the  full  IMS  primary  network. 


The  detection  capability  of  the  IMS  primary  seismic  network  as  of  July  2001 

It  is  of  particular  interest  to  compare  the  detection  capability  of  the  current  IMS  network  with 
the  projected  performance  of  the  49  stations  of  the  full  network.  As  of  7  July  2001  the  primary 
seismic  IMS  network  consisted  of  35  stations,  but  several  of  these  stations  did  not  provide  any 
data.  We  have  in  Fig.  6.1.5  plotted  the  background  noise  levels  for  a  time  interval  without  any 
significant  seismic  signals,  and  we  see  that  several  stations  are  down  or  providing  data  gaps. 
The  corresponding  three-station  detection  capability  is  shown  in  Fig.  6.1.6,  and  we  again  notice 
the  good  detectability  associated  with  the  high  performance  arrays  in  northern  Europe,  North 
America  and  Australia. 

We  have  in  Fig.  6.1.7  plotted  the  difference  in  detectability  between  the  current  IMS  primary 
seismic  and  the  projected  performance  of  the  full  network.  The  average  difference  is  only  0.06 
magnitude  units,  and  the  largest  differences  are  found  in  southwest  Asia  with  a  maximum  value 
of  0.47.  The  constant  noise  levels  assigned  to  the  planned  or  non-operational  stations  are  taken 
from  typical  noise  levels  at  stations  assumed  to  have  similar  noise  characteristics.  The  rela¬ 
tively  large  uncertainty  associated  with  the  assigned  constant  noise  levels  is  also  reflected  in  the 
detectability  estimates. 
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Fig.  6. 1 .5.  Noise  levels  of  current  IMS  primary  seismic  stations  on  7  July  2001. 
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Fig.  6.1.6.  Three-station  detection  capability  of  IMS  primary  stations  on  7  July  2001. 
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Fig.  6.1.7.  Difference  in  detection  capability  between  the  current  operational  primary  seismic  net 
work  (Fig.  6.1.6)  and  full  network  (Fig  6.1.4). 
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The  detection  capability  of  the  IMS  primary  seismic  network  in  the  coda  of  a  large 
earthquake 

At  09:38:43.5  on  7  July  2001  a  m^,  5.7,  7.6  earthquake  occurred  near  the  coast  of  Peru.  Fig. 

6.1.8  shows  the  average  three-station  network  detection  capability  for  a  five  minute  time  inter¬ 
val  that  starts  1  minute  and  20  seconds  after  the  origin  time  of  the  event.  Notice  the  reduced 
detectability  for  the  areas  around  the  epicenter. 


Fig.  6.1.8.  Three-station  detection  capability  of  IMS  primary  stations  1  minute  and 
20  seconds  after  the  occurrence  of  a  m^  5.7  earthquake  near  the  coast  of  Peru. 


In  Fig.  6.1.9  we  have  plotted  the  difference  in  detectability  between  the  time  interval  with  the 
earthquake  signals  (Fig.  6.1.8)  as  compared  to  the  detectability  during  quiet  background  noise 
conditions  (Fig.  6.1.6).  For  large  regions,  including  South  America  and  adjacent  areas,  the 
detection  performance  is  reduced  by  more  than  1  magnitude  unit.  The  maximum  difference  is 
1 .925  units  in  the  vicinity  of  the  epicenter.  As  time  passes,  the  seismic  signals  will  propagate  to 
longer  distances  and  reduce  the  detectability  for  larger  regions  of  the  Earth.  However,  due  to 
the  signal  attenuation  the  degradation  of  the  detection  capability  will  be  less  than  for  regions 
closer  to  the  event. 
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Fig.  6.1.9.  Difference  in  detection  capability  between  the  levels  given  in 
Fig.  6.1.8  and  6.1.6. 


Conclusions 

We  have  found  that  the  35  stations  of  the  primary  seismic  IMS  network  operational  as  of  July 
2001  have  a  three-station  detection  capability  that  is  quite  close  to  the  projected  performance  of 
the  49  stations  of  the  network.  The  largest  difference  of  about  0.5  m^,  units  is  found  in  south¬ 
west  Asia.  This  can  be  explained  by  the  fact  that  most  of  the  sensitive  array  stations  of  the  net¬ 
work  are  already  in  place.  Another  uncertainty  is  the  noise  levels  assigned  to  the  planned 
stations  and  to  the  stations  that  did  not  provide  real  noise  data  during  the  period  under  investi¬ 
gation.  The  assignment  of  realistic  noise  levels  for  planned  stations  is  therefore  a  topic  requir¬ 
ing  additional  studies. 

This  study  also  confirms  the  results  provided  by  Kvaema  and  Ringdal  (1999)  that  large  earth¬ 
quakes  and  the  corresponding  coda  energy  can  temporarily,  over  tens  of  minutes,  significantly 
reduce  the  detection  capability  of  the  IMS  network.  During  such  conditions  the  use  of  high-fre¬ 
quency  regional  data,  combined  with  less  stringent  event  formation  criteria  (1  or  2  station  only) 
will  be  important  for  CTBT  verification  purposes. 
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6.2  Travel  times  and  attenuation  relations  for  regional  phases  in  the 
Barents  Sea  region 

This  research  is  conducted  under  contract  DTRAOl-OO-C-0107. 

A  database  containing  42  events  in  the  Barents  Sea  region  has  been  compiled  and  analyzed 
with  the  aim  of  evaluating  crustal  models,  travel  times  and  attenuation  relations  in  the  context 
of  performing  regional  detection  threshold  monitoring  of  this  region.  The  42  events  are  mostly 
located  around  the  circumference  of  the  study  area  due  to  the  virtually  aseismic  nature  of  the 
Barents  Sea  itself.  Regional  and  S^  phases  were  observable  for  most  events  in  the  database, 
while  Pg  and  Lg  phases  were  only  observable  for  events  with  ray  paths  within  continental  crust. 
This  corroborates  a  number  of  previous  observations  of  Lg-wave  blockage  across  the  Barents 
Sea.  Three  existing  velocity  models  were  evaluated,  with  a  model  having  slightly  lower  S- 
velocities  than  earlier  assumed  in  the  upper  mantle  giving  the  overall  best  fit  to  the  observed 
arrivals.  In  order  to  estimate  magnitudes,  short  term  average  (STA)  and  spectral  ^plitude  val¬ 
ues  were  calculated  in  several  frequency  bands  for  all  phase  arrivals  in  the  data  base.  There 
were  no  significant  differences  between  spectral  and  STA  amplitudes,  so  the  latter  were  used  as 
this  parameter  is  more  efficient  to  calculate  in  real-time  processing.  An  inversion  was  per¬ 
formed  in  order  to  determine  a  Pj,  and  S„  attenuation  relation  specific  for  this  region.  The 
resulting  magnitudes  based  on  P^  and  Sn  phases  gave  an  internally  consistent,  reasonably  stable 
set  of  vdues,  which  can  be  calibrated  towards  any  existing  global  or  regional  scale.  An  attenu¬ 
ation  relation  was  also  determined  for  the  Lg  phase,  but  the  low  number  of  amplitude  readings 
in  this  case  renders  the  results  less  reliable. 

Introduction 

The  Barents  Sea  region  is  an  area  which  is  of  particular  interest  in  the  context  of  the  compre¬ 
hensive  nuclear  test  ban  treaty  (CTBT),  as  it  contains  the  former  Soviet  Union  nuclear  test  sites 
on  Novaya  Zemlya.  For  this  reason,  it  is  of  interest  to  perform  regional  seismic  detection 
threshold  monitoring  in  order  to  continuously  assess  the  upper  limit  magnitude  of  events  that 
could  go  undetected,  and  also  to  provide  automatic  locations  for  seismic  events  in  the  region 
with  the  best  possible  precision.  The  Threshold  Monitoring  (TM)  method  (e.g.  Ringdal  and 
Kvaema,  1989;  1992;  Kvaema  and  Ringdal,  1999)  uses  continuous  seismic  data  from  a  network 
of  stations  to  calculate  a  threshold  magnitude  for  each  point  in  a  grid,  for  which  an  event  would 
have  90%  probability  of  detection  by  the  network.  It  is  thus  able  to  take  varying  noise  levels, 
interfering  signals  from  teleseismic  events,  data  gaps,  particularly  favorable  site-station  trans¬ 
mission  properties  etc.  into  account,  which  traditional,  static  capability  maps  are  incapable  of 
doing  (Ringdal  and  Kvaema,  1992).  However,  for  the  threshold  monitoring  to  be  accurate,  reli¬ 
able  travel  time  curves  and  attenuation  relations  are  required  for  the  target  area.  Ideally,  calibra¬ 
tion  events  should  be  available  for  all  station/phase  combinations  at  each  target  point.  However, 
in  practice,  particularly  in  areas  of  low  seismic  activity  such  as  the  Barents  Sea,  it  is  necessary 
to  rely  on  travel  time  curves  calculated  from  a  regional  velocity  model  and  regional  attenuation 
relations.  Still,  a  good  coverage  of  seismic  events  is  a  requirement  in  order  to  determine  mean 
regional  values  of  such  relations.  In  order  to  provide  the  necessary  data  for  this  region,  we  have 
compiled  a  database  of  recent  events,  and  have  attempted  to  extract  the  necessary  information 
in  spite  of  the  less  than  optimal  coverage. 


57 


NORSAR  Sci.  Rep.  2-2001 


December  2001 


Seisnwtectonic  setting 

The  Barents  Sea  is  an  epicontinental  sea,  bordering  on  the  Precambrian  and  Caledonian  crust 
of  northern  Fennoscantha,  the  Kola  Peninsula  and  northern  Siberia  to  the  south,  and  young  pas¬ 
sive  margins  to  the  north  and  west,  formed  during  the  Cenozoic  opening  of  the  Eurasian  Basin 
and  the  Norwegian-Greenland  Sea  respectively  (e.g.  Faleide  et  ai,  1993).  Large  parts  of  the 
western  margin  and  the  oceanic  crust  to  the  west  is  covered  by  a  huge  sediment  wedge  depos¬ 
ited  after  the  opening  of  the  Norwegian-Greenland  Sea,  primarily  during  the  Pliocene-Pleis¬ 
tocene  age  (Eidvin  et  al,  1993).  The  interior  of  the  Barents  Sea  is  underlain  by  large 
thicknesses  of  Upper  Paleozoic  to  Cenozoic  sediments,  with  accentuated  Moho  relief,  leaving 
as  little  as  ~4  km  of  crystalline  basement  under  the  deepest  basins  (Jackson  et  ai,  1990;  Faleide 
etai,  1993). 

The  internal  Barents  Sea  has  very  little  seismic  activity,  implying  a  stable  tectonic  situation 
(Bungum  and  Lindholm,  1996).  This  has  also  been  interpreted  as  a  consequence  of  a  relatively 
weak  ridge-push  force  counteracted  by  tensional  stresses  related  to  Pliocene-Pleistocene  ero- 
sional  unloading  (Fiedler  and  Faleide,  1996).  The  activity  is  significantly  higher  along  the 
western  Barents  Sea  margin,  which  is  under  weak  oblique  (NW-SE)  compression  from  the 
north-Atlantic  ridge,  with  the  large  PUocene-Pleistocene  sediment  load  likely  being  the  most 
important  stress-generating  mechanism  in  these  areas  (Byrkjeland  etai,  2000). 

Northern  Fennoscandia  and  the  Kola  Peninsula  are  areas  of  intermediate  seismic  activity,  tend¬ 
ing  mainly  toward  earthquakes  in  the  upper  half  of  the  crust.  There  appears  to  be  some  activity 
connected  to  a  postglacial  reverse  and  strike-slip  fault  system  in  northern  Fennoscandia  (Bun- 
gum  and  Lindholm,  1996),  the  stress  tensor  appears  to  have  NW-SE  compression  approxi¬ 
mately  corresponding  to  the  expected  ridge-push  direction  (Hicks  et  al.,  2000).  There  is  also 
some  seismic  activity  along  the  northern  coast  of  the  Kola  Peninsula,  while  onshore  northern 
Siberia  has  lower  levels.  The  seismic  activity  in  northern  Fennoscandia  is  most  likely  tied  to  a 
combination  of  the  tectonic  ridge  push  force  and  constructive  postglacial  uplift  stresses  (Bun- 
gum  and  Lindholm,  1996).  Novaya  Zemlya  and  the  Kara  Sea  to  the  east  has  some  sporadic 
earthquake  activity,  as  evidenced  by  some  interesting  events  in  recent  years  (e.g.  Marshall  et 
ai,  1989;  Ringdal  etai,  1997;  Bowers  etai,  2001;  Schweitzer  and  Kennett,  2002). 

Database  of  seismic  events 

A  total  of  42  seismic  events  in  the  Barents  Sea  and  northern  Fennoscandia  were  used  as  the 
basis  for  this  study.  The  events  were  selected  to  provide  the  best  possible  ray  path  coverage  of 
the  cmst,  although  since  large  areas  of  the  Barents  Sea  are  in  practice  aseismic,  there  is  a  con¬ 
centration  of  activity  on  the  outer  parts:  Svalbard,  western  Barents  Sea,  northern  Fennoscandia, 
the  Kola  Peninsula  and  Novaya  Zemlya.  The  database  consists  of  earthquakes,  mining  blasts 
and  other  explosions  (both  chemical  and  nuclear),  and  some  events  of  unknown  origin.  Areas 
with  large  numbers  of  seismic  events  with  similar  locations  were  represented  with  a  single 
event,  in  order  to  avoid  introducing  any  bias  towards  these  areas  in  the  subsequent  inversion. 

In  addition  to  waveform  data  from  seismic  array  stations  on  mainland  Norway  (ARCES, 
NORES),  we  have  also  used  available  array  data  from  Finland  (FINES),  Svalbard  (SPITS),  and 
the  Russian  Kola  peninsula  (APA).  Data  from  the  Amderma  (AMD)  station  in  northern  Russia 
were  available  for  15  of  the  events,  waveform  data  were  also  retrieved  from  the  IRIS  consor¬ 
tium  for  the  Global  Seismographic  Network  (GSN)  seismic  stations  in  Kevo,  Northern  Finland 
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(KEV),  Lovozero  on  the  Kola  peninsula  (LVZ)  and  Ny  Alesund  on  the  Svalbard  archipelago 
(KBS).  The  location  of  these  nine  stations  used  are  shown  along  with  the  locations  of  the  42 
selected  seismic  events  in  Fig.  6.2.1.  The  waveform  data  were  re-analyzed,  with  emphasis 
placed  on  consistent  phase  identification  and  onset  time  picking  of  observable  phases.  Only 
phases  that  could  be  clearly  observed  and  identified  were  analyzed  and  used.  A  preliminary 
relocation  of  the  events  was  performed  during  the  analysis,  using  the  ‘Fennoscandia’  crustal 
model  (Mykkeltveit  and  Ringdal,  1981)  which  is  the  model  routinely  used  by  NORSAR  for 
seismic  event  location  in  northwestern  Europe  and  the  Barents  Sea. 

An  investigation  of  which  phases  were  observable  on  the  different  stations  from  various  source 
areas  was  performed,  in  order  to  evaluate  the  geographical  coverage  for  each  phase.  As 
expected,  crustal  phases  (Pg  and  Lg)  are  in  general  only  observable  for  paths  that  travel  more  or 
less  exclusively  within  the  shield  areas,  i.e.  from  onshore  or  coastal  events  in  Fennoscandia  and 
Svalbard.  This  confirms  previous  observations  indicating  the  blockage  of  Lg  phases  that  have 
paths  crossing  large  sedimentary  structures  such  as  those  encountered  in  the  Barents  Sea  (e.g. 
Zhang  and  Lay,  1994;  Baumgardt,  2001 ;  Bowers  et  al,  2001).  Pn  and  S^  phases  are  observable 
for  most  events  with  distances  greater  than  2-3  degrees.  Fig.  6.2.2  shows  phase  maps  for  the 
ARCES  array  in  northern  Norway. 

Crustal  models  and  travel  times 

In  order  to  be  able  to  predict  phase  arrival  times  for  a  given  origin  in  a  threshold  monitoring 
application,  an  accurate  crustal  velocity  model  is  required.  An  evaluation  of  several  available 
crustal  models  was  therefore  performed.  Observed  and  theoretical  travel  times  for  two  models 
are  shown  in  Fig.  6.2.3,  plotted  according  to  epicentral  distance  after  relocation  with  the  corre¬ 
sponding  model  using  HYPOSAT  (Schweitzer,  2(X)1).  All  depths  were  fixed  at  10  km  for  this 
comparison.  The  Fennoscandia  model  (Mykkeltveit  and  Ringdal,  1981),  shown  in  Table  6.2.1, 
gave  a  reasonably  good  travel  time  residuals,  but  there  were  some  systematic  discrepancies  as 
clearly  visible  in  Fig.  6.2.3.  The  Lg-phase  arrivals  have  consistent  and  quite  large  negative 
residuals  on  average,  while  the  phases  trend  towards  positive  residuals.  The  P^  arrivals  have 
quite  large  negative  residuals,  in  particular  around  the  6  - 12  degree  distance  range.  Two  mod¬ 
els  from  Schweitzer  and  Kennett  (2002),  named  BAREY  and  BAREZ  (Table  1)  were  also 
tested,  with  the  BAREY  model  giving  the  lowest  travel  time  residuals  overall  (Fig.  6.2.3). 
These  two  models  were  adapted  from  a  model  developed  for  the  Barents  and  Kara  Seas  by 
Kremenetskaya  et  al.  (2001)  by  making  minor  P-velocity  adjustments  and  varying  the  P/S  ratio 
in  the  upper  mantle  (Schweitzer  and  Kennett,  2002).  The  main  differences  between  these  two 
models  are  the  S  velocities  in  the  upper  mantle,  which  of  course  have  the  greatest  effect  on 
regional  S  phases  propagating  at  these  depths  (S^).  The  BAREY  model  has  the  lowest  S  veloc¬ 
ities  in  the  upper  mantle  while  BAREZ  has  the  highest.  The  Fennoscandia  model  has  upper 
mantle  S  velocities  between  the  other  two  models.  As  shown  in  Table  6.2.1,  the  velocity  pro¬ 
files  for  both  P  and  S  within  the  crust  are  virtually  identical,  with  the  Conrad  discontinuity  at  16 
km  and  the  Mohorovicic  discontinuity  at  40  km  (Fennoscandia)  and  41  km  (BAREY/BAREZ) 
depth. 

The  small  differences  in  upper  mantle  P  velocities  (Table  6.2. 1)  serve  to  reduce  the  observed  Pj, 
residuals  after  relocation.  The  resulting  small  differences  in  epicenter  locations  for  the  different 
models  also  have  an  effect  on  the  travel  times  of  crustal  phases,  most  visible  for  the  Lg  onsets 
(Fig.  6.2.3).  The  BAREY  model  appears  to  be  particularly  accurate  for  paths  from  Novaya 
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Zemlya  to  Fennoscandia,  which  is  the  setting  for  which  it  was  developed.  The  BAREZ  model 
is  slightly  better  for  paths  crossing  the  Barents  Sea  from  Novaya  Zemlya  to  the  north-west 
(towards  Svalbard  and  Bj0m0ya),  and  also  to  the  south  (towards  Amderma)  (Schweitzer  and 
Kennett,  2(X)2).  The  BAREY  model  however  provided  the  smallest  overall  travel  time  residuals 
for  all  paths  in  our  data  set.  The  geographic  distribution  of  the  travel  time  residuals  for  the  Fen¬ 
noscandia  and  BAREY  models  is  shown  in  Fig.  6.2.4,  showing  the  generally  lower  residuals,  in 
particular  for  the  Sn  arrivals,  obtained  using  the  BAREY  model. 

Attenuation  relations  and  magnitudes 

In  order  to  provide  threshold  magnitudes  for  a  given  location,  it  is  essential  to  be  able  to  con¬ 
vert  observed  amplitudes  to  event  magnitudes  at  the  target  location.  To  this  end  an  attenuation 
relation  from  Jenkins  etal.  (1998),  based  on  Sereno  (1990)  is  used: 

where  A  is  amplitude  in  nanometers,  A  is  the  epicentral  distance  in  km,/is  the  logarithmic  cen¬ 
ter  frequency  of  the  passband  at  which  the  amplitude  reading  is  taken,  while  a  and  6  are  phase- 
dependent  constants.  The  value  200  represents  a  reference  distance  where  the  geometrical 
spreading  changes  from  spherical  to  a  more  complicated  function  based  on  the  phase  analyzed 
(Sereno,  1990).  Due  to  the  fact  that  the  decrease  in  seismic  amplitude  with  epicentral  distance 
is  a  combined  effect  of  geometrical  spreading,  anelasticity  and  scattering,  care  must  be  taken  in 
relating  these  coefficients  to  physical  properties  of  the  medium  (cf.  Alsaker  eta/.,  1991).  How¬ 
ever,  this  is  not  necessary  for  our  purpose,  as  long  as  the  relation  as  a  whole  is  able  to  describe 
the  amplitude/distance  relation  in  a  consistent  and  reasonably  accurate  manner.  The  relation 
was  developed  using  spectral  amplitudes,  however  for  continuos  processing  it  is  more  efficient 
to  use  time-domain  short  term  average  (STA)  amplitudes.  A  comparison  of  STA  amplitudes 
and  spectral  amplitudes  in  the  same  frequency  bands  for  the  selected  events  showed  that  they 
are  in  practice  equivalent  with  regard  to  magnitude  calculation  (Asja  ~  Aspec)-  STA  ampli¬ 
tudes  were  calculated  in  three  frequency  bands  (2-4, 3-6  and  4-8  Hz)  corresponding  to  three  of 
the  frequency  ranges  used  in  Jenkins  et  al.  (1998),  using  a  moving  window  with  a  step  length 
of  one  sample  and  window  lengths  of  2  s  length  for  P„  and  5  s  for  and  Lg.  The  maximum 
STA  amplitude  within  7.5  s  (?„),  10  s  (Sn)  and  15  s  (Lg)  after  the  observed  phase  onset  was 
then  selected.  Only  amplitude  readings  with  a  minimum  signal-to-noise  ratios  (SNR)  of  3.0  or 
greater  for  and  S„  and  2.5  or  greater  for  Lg,  compared  to  a  noise  window  of  5  s  length  start¬ 
ing  10  s  ahead  of  the  observed  onset  time  have  been  used.  The  SNR  for  Lg  is  almost  always 
lower  than  for  Sn  or  as  the  Lg  phase  arrives  within  the  coda  of  the  earlier  onsets.  Similarly, 
Sn  arrivals  generally  show  lower  SNRs  than  Pn  arrivals,  although  this  can  also  depend  on 
source  effects  and  the  relative  noise  situation  at  the  seismic  station. 

Magnitudes  calculated  using  equation  (1)  did  reveal  some  inconsistencies  between  magnitudes 
calculated  from  different  frequency  bands,  phases  and  stations,  as  shown  in  Fig.  6.2.5.  The  a 
and  b  coefficients  used  in  this  care  were  determined  using  data  from  eastern  North  America, 
central  Asia  and  Australia  by  Jenkins  et  al.  (1998)  (Table  6.2.2),  also  applied  by  Bowers  et  al. 
(2001)  for  the  Barents  Sea.  For  ARCES  it  is  clear  that  magnitudes  calculated  from  STA  values 
in  the  2-4  Hz  passband  are  systematically  higher,  by  about  0.3  magnitude  units  on  average, 
than  magnitudes  calculated  in  the  3-6  Hz  passband,  which  again  are  almost  0.2  magnitude  units 
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higher  on  average  than  magnitudes  calculated  from  the  4-8  Hz  passband.  Similar  discrepancies 
were  observed  at  the  other  stations.  As  Fig.  6.2.5  also  shows,  there  is  a  quite  large,  systematic 
offset  between  the  Pj,  and  phase  magnitudes  calculated  within  the  same  frequency  band  at 
the  same  station.  In  addition,  there  are  large  discrepancies  between  magnitudes  calculated  from 
the  same  phase/frequency  range  at  different  stations.  This  implies  that  the  relation  does  not 
accurately  represent  the  attenuation  of  seismic  waves  for  this  set  of  reference  events. 

In  order  to  provide  an  attenuation  relation  that  better  explains  the  measured  amplitudes  in  the 
database,  and  also  to  determine  station  corrections  for  the  nine  stations  used,  an  inversion  of 
the  available  amplitude  data  was  performed.  Using  equation  (1),  the  STA  amplitude  for 
phase  i  from  event  j  measured  at  station  k  may  be  expressed  as: 

where  5,^.  is  the  station  correction  for  phase  i  at  station  k,  and  Mj  is  the  event  magnitude  for 
event y.  The  parameters  Aj  and  bi  are  the  coefficients  from  equation  (1)  for  phase  i.  Equation  (2) 
represents  a  set  of  linear  equations,  and  can  be  written  in  matrix  form: 


Ax  =  D  (3) 

where  A  represents  the  Jacobi  matrix,  x  the  vector  of  unknowns  and  D  is  the  data  vector.  We 
used  standard,  least-squares  techniques  to  solve  equation  (2). 

A  total  of  863  amplitude  readings  with  acceptable  SNR  (>3.0)  from  (548  readings)  and  Sjj 
(315  readings)  phases  from  the  40  seismic  events  that  had  sufficient  data  to  be  used,  were 
inverted  in  order  to  yield  the  two  relation  coefficients  and  nine  station  corrections  for  each  of 
the  two  phases,  in  addition  to  relative  magnitudes  for  the  40  events.  The  individual  events  had 
between  four  and  40  amplitude  readings  contributing  to  the  solution.  The  a,-  and  £»,•  values  deter¬ 
mined  through  this  inversion  are  given  in  Table  6.2.2,  while  the  corresponding  station  correc¬ 
tions  Sii^  are  listed  in  Table  6.2.3.  New  magnitudes  calculated  using  these  results  are  plotted  in 
Fig.  6.2.6,  showing  that  the  mean  offset  was  reduced  in  all  cases  compared  to  the  original  cal¬ 
culations  (Fig.  6.2.5).  The  scatter,  expressed  as  the  standard  deviation,  was  marginally  higher 
for  some  combinations,  and  reduced  for  others. 

The  magnitudes  for  the  P^  and  Sq  phases  calculated  using  the  inversion  results  are  internally 
consistent,  and  can  thus  be  tied  in  to  any  other  magnitude  scale.  To  calibrate  the  scale  in  this 
case  we  have  used  five  events  in  the  Novaya  Zemlya  region  with  body  wave  magnitudes 
ranging  from  m^  2.4  to  5.6  (Table  4)  that  have  previously  been  studied  in  great  detail  (e.g. 
Ringdal  and  Kremenetskaya,  1999),  giving  a  constant  offset  of  -0.33.  Applying  this  offset 
value  and  the  estimated  phase-dependant  station  corrections  Si  (listed  in  Table  6.2.3)  to 
equation  (1),  the  final  magnitude  relation  for  P^  and  Sj,  STA  amplitudes  yields: 


M 


=  log.4+log((A)“'^  + 


L{Pn,SA 


0.33 


(4) 


where  A  is  observed  STA  amplitude  within  a  frequency  band  with  logarithmic  center  frequency 
/,  A  is  epicenfral  distance,  and  a,  and  bi  are  phase-dependant  constants  given  in  Table  6.2.2.  The 
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Pn  and  Sj,  magnitudes  calculated  using  equation  (4)  were  in  general  comparable  to  previously 
reported  magnitudes  for  other,  well  studied  events,  as  shown  in  Table  6.2.4. 

A  separate  inversion  was  attempted  in  order  to  determine  coefficients  for  the  Lg  phase  for 
events  and  stations  located  within  Fennoscandia,  as  Lg  observations  are  hmited  to  this  area. 
The  inversion  gave  seemingly  acceptable  results  when  applied  to  our  small  group  of  events, 
however,  attenuation  was  extremely  low  for  distances  over  around  500  km.  Given  this  contrast 
to  the  general  Ml  relation  for  Norway  from  Alsaker  et  al.  (1991),  the  inversion  in  this  case 
appears  to  be  unable  to  adequately  constrain  the  required  parameters.  Data  from  more  event/ 
station  combinations  in  Fennoscandia  are  needed  in  order  to  provide  a  valid  Lg  STA-based 
magnitude  relation  for  Fennoscandia.  Lg  magnitude  relations  are  anyhow  of  little  interest  with 
regard  to  magnitude  estimation  for  events  in  the  Barents  and  Kara  Seas  due  to  the  strong  atten¬ 
uation  or  blocking  of  Lg  energy  observed  for  events  in  this  area. 

Fig.  6.2.7  shows  comparisons  of  corrected  network  magnitudes  as  calculated  from  individual 
phases.  There  are  still  clear  discrepancies  between  and  magnitudes  for  different  source 
areas,  which  seem  to  be  the  result  of  different  ray  paths.  Events  that  predominantly  have  ray 
paths  within  Fennoscandia  have  larger  Sn  magnitudes  (relative  to  ?„),  while  the  opposite  is  true 
for  events  that  have  ray  paths  crossing  the  sediment  basins  of  the  Barents  Sea  (Novaya  Zemlya/ 
Kara  Sea  and  the  western  Barents  Sea/Mid- Atlantic  ridge  areas).  Table  4  lists  the  date,  location 
and  P„,  Sn,  and  overall  network  magnitudes  with  number  of  readings  and  RMS  values  for  the 
42  events  in  the  database.  Magnitudes  (Ml)  estimated  by  NORSARs  GBF  system  (Ringdal  & 
Kvaema,  1989)  and  reference  magnitudes  (m^)  for  the  five  Novaya  Zemlya  events  from  Ringdal 
and  Kremenetskaya  (1999)  are  also  given. 

Data  from  the  Amderma  station 

The  seismic  station  in  Amderma  (AMD)  has  a  unique  and  important  location  with  regard  to 
studying  seismic  activity  in  the  eastern  Barents  Sea/Kara  Sea  area  (Fig.  6.2.1).  Data  from  this 
station  have  not  been  generally  available,  but  waveform  data  for  a  number  of  events  in  the  Bar¬ 
ents  Sea  region  were  collected  as  part  of  a  joint  project  between  NORSAR  and  the  the  Kola 
Regional  Seismological  Centre  in  Apatity.  These  data  covered  a  total  of  15  of  the  42  events  in 
our  database,  of  which  phase  arrivals  could  be  identified  for  13  of  the  events.  These  arrivals 
represent  an  important  contribution  to  the  evaluation  of  the  models.  The  epicentral  distances 
for  the  majority  of  the  events  in  the  database  are  larger  than  for  the  other  stations  used  in  this 
study,  so  only  nine  events  had  onsets  with  SNRs  large  enough  to  be  included  in  the  inversion  of 
the  attenuation  relations.  Magnitudes  calculated  from  amplitude  observations  at  AMD  using 
the  coefficients  determined  by  the  magnitude  inversion  in  this  study  appear  to  correlate  on 
average  quite  well  to  the  other  stations,  although  the  scatter  is  somewhat  higher.  The  higher 
scatter  is  most  likely  an  effect  of  the  significantly  different  location  (and  thereby  differences  in 
ray  paths),  and  also  the  relatively  low  number  of  good  amplimde  readings  for  this  station  com¬ 
pared  to  the  other  stations  used  in  this  study. 

Conclusions 

The  following  main  conclusions  can  be  drawn  from  this  study  of  regional  phases  in  the  Barents 
Sea  region: 
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The  pattern  of  Lg  arrivals  and  associated  amplitudes  support  the  previously  published  indica¬ 
tions  that  the  deep  sediment  basins  and  Moho  topography  under  the  Barents  Sea  efficiently 
blocks  Lg  wave  energy  from  crossing. 

The  ‘BAREY’  model  (Table  6.2.1)  from  Schweitzer  and  Kenned:  (2002),  based  on  a  model  for 
the  Barents  Sea  area  from  Kremenetskaya  et  al.  (2(X)1),  provides  the  smallest  overall  travel 
time  residuals  when  locating  events  within  the  vicinity  of  the  Barents-  and  Kara  Seas. 

The  attenuation  in  the  Barents  Sea  region  differs  somewhat  from  that  observed  in  other  stable 
tectonic  regions,  as  evidenced  by  the  fact  that  the  coefficients  given  by  Jenkins  et  al.  (1998)  for 
such  regions  do  not  give  consistent  magnitudes  across  frequencies,  phases  and  stations  for  our 
amplitude  observations  from  the  events  in  the  Barents  Sea  region. 

Amplitude  inversion  has  been  used  in  this  study  to  resolve  new  coefficients  and  station  correc¬ 
tions  for  estimating  Ml  magnitudes  from  STA  amplitude  observations  for  and  S^  phases  in 
the  Barents  Sea  region  (Tables  6.2.2  and  6.2.3).  The  number  of  Lg  amplitude  observations  for 
raypaths  within  northern  Fennoscandia  was  too  low  to  provide  a  usable  magnitude  relation  for 
Lg.  A  futme  study  using  a  greater  number  of  continental  events  could  most  likely  provide  a 
more  stable  relation  for  STA  based  Lg  magnitudes. 

The  magnitudes  determined  using  the  relation  in  equation  (4)  are  self-consistent,  and  can  be 
tied  in  to  any  other  global,  regional  or  local  magnitude  scale.  The  relation  given  in  equation  (4) 
has  been  adjusted  in  order  to  comply  with  body  wave  magnitudes  for  five  events  at  Novaya 
Zemlya  from  Ringdal  and  Kremenetskaya  (1999). 

The  seismic  station  in  Amderma  can  be  tied  in  to  the  regional  network  in  Fennoscandia  and  on 
the  Svalbard  archipelago  using  an  appropriate  crustal  model,  and  is  able  to  provide  important 
information  regarding  the  location  of  events  in  the  eastern  parts  of  the  Barents  Sea  and  the  Kara 
Sea  (Schweitzer  and  Kennett,  2002).  Magnitudes  calculated  at  this  station  are  on  the  whole 
quite  consistent  with  the  other  observations. 

The  BAREY  velocity  model  (Schweitzer  and  Kennett,  2002),  and  attenuation  relation  given  in 
equation  (4),  currently  provide  the  best  possible  locations  and  magnitudes  for  events  in  the  Bar¬ 
ents  Sea  region,  and  can  be  used  for  a  Threshold  Monitoring  implementation  for  this  area. 
When  more  and  better  quality  data  become  available  it  may  be  possible  to  further  improve 
these  results,  although  the  aseismic  nature  of  the  Barents  Sea  will  continue  to  pose  problems 
with  regard  to  a  detailed  study  of  seismic  wave  propagation  within  the  Barents  Sea  itself. 

E.C.  Hicks 
T.  Kvaerna 
S.  Mykkeltveit 
J.  Schweitzer 
R  Ringdal 
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Fennoscandia 

Depth  Vp  Vg 

(km) 

Depth 

(km) 

BAREY 

Vp  Vs 

Depth 

(km) 

BAREZ 

Vp  Vs 

0.0 

6.20 

3.58 

0.0 

6.20 

3.58 

0.0 

6.20 

3.58 

16.0 

6.20 

3.58 

16.0 

6.20 

3.58 

16.0 

6.20 

3.58 

16.0 

6.70 

3.87 

16.0 

6.70 

3.87 

16.0 

6.70 

3.87 

40.0 

6.70 

3.87 

41.0 

6.70 

3.87 

41.0 

6.70 

3.87 

40.0 

8.15 

4.705 

41.0 

8.10 

4.58 

41.0 

8.10 

4.71 

70.0 

8.23 

4.65 

70.0 

8.23 

4.78 

95.0 

8.15 

4.705 

95.0 

8.25 

4.763 

210.0 

8.30 

4.792 

210.0 

8.26 

4.67 

210.0 

8.26 

4.80 

Table  6.2.1.  The  Fennoscandia  (Mykkeltveit  and  Ringdal,  1981 ),  BAREY  and  BAREZ  (Schweitzer 
and  Kennett,  2002)  crustal  models.  The  main  differences  between  the  models  are  the  slightly 
different  velocities  in  the  uppermost  mantle. 


Pn 

a  b 

s„ 

a  b 

Jenkins  et  al.  (1998) 

0.08 

1.44 

0.12 

1.85 

This  study 

0.36±0.0 

2 

0.88±0.0 

9 

0.4110.0 

2 

0.6310.1 

1 

Table  6.2.2.  The  inversion  results  for  the  aandb  coefficients  (±la)for  P„  and  S„  phases  used  in  the 
attenuation  relation  (equations  4).  The  values  from  Jenkins  et  al.  (1998)  are  shown  for 
comparison.  Note  that  the  new  coefficients  have  a  stronger  frequency  dependency  (higher  a 
values). 
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Station 

Pn 

Sn 

AMD 

0.05+0.07 

0.01+0.07 

APA 

0.15±0.04 

-0.16±0.05 

ARCES 

0.13±0.04 

0.0310.04 

FINES 

0.13±0.05 

-0.0110.05 

KBS 

-0.07±0.04 

-0.0210.07 

KEV 

-0.09±0.04 

-0.3710.05 

LVZ 

0.09±0.04 

-0.2910.05 

NORES 

0.13±0.05 

0.0910.07 

SPITS 

0.05+0.04 

0.1210.09 

Table  6.2.3.  .Station  corrections  (±lo)  in  magnitude  units  obtained  as  part  of  the  inversion. 
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Pn(ML) 


1988-235:16.20.00.13 

1990-167:12.43.26.91 

1990- 297:14.58.06.49 

1991- 157:12.46.11.32 

1991- 236:10.56.29.65 

1992- 366:09.29.25.84 

1993- 005:10.19.34.40 
1995-013:04.34.08.58 
1995-021:01.58.07.21 
1995-054:21.50.00.15 
1995-063:18.29.04.36 
1995-101:20.07.23.02 
1995-102:08.18.52.49 
1995-133:22.38.51.04 
1995-161:18.45.34.11 
1995-162:19.27.15.07 
1995-164:19.22.38.41 
1995-184:12.49.32.76 
1995-185:03.26.24.87 
1995-241:22.12.19.06 
1995-261:03.26.06.18 
1995-313:01.10.23.59 

1995- 329:19.41.26.68 

1996- 013:17.17.23.57 
1996-021:02.16.32.03 
1996-218:20.04.38.23 
1996-272:06.08.47.35 
1996-301:23.55.17.21 

1996- 361:04.44.15.67 

1997- 122:07.31.17.23 
1997-228:02.11.00.36 
1997-279:12.33.27.04 
1997-279:21.17.31.53 
1997-279:21.29.18.26 
1999-229:04.44.37.19 

1999- 290:12.07.16.70 

2000- 225:07.30.41.31 
2000-316:02.01.06.65 
2000-341:00.34.40.15 

2000- 360:03.50.28.47 

2001- 090:11.30.55.35 
2001-122:15.59.43.93 


Mag. 


5.10  1 
3.75  2 

5.65  1 

3.18  2 
3.48  1 
3.02  1 

3.97  4 

3.45  8 
2.99  5 

2.18  2 

3.65  6 

3.67  5 
2.35  4 

3.54  9 

2.55  3 
3.29  3 

3.51  4 

2.90  3 

3.71  2 
3.75  4 
2.86  6 

2.91  6 
3.65  6 

2.46  2 

3.72  7 
3.13  7 
1.95  3 
3.45  4 

2.98  4 
2.86  6 
3.42  5 
3.77  4 
3.26  5 
3.83  4 

3.99  6 

3.55  7 

3.41  8 

3.52  6 
3.86  3 

3.68  8 
3.90  8 
2.93  5 


S„(Ml) 

Network  (Ml) 

Mag.  no  RMS 

Mag.  no  RMS 

GBF(ML)|Ref. 
Mag.  no  (nib) 


3.98  2 

5.16  1 

3.41  2 
3.65  1 
2.71  1 

4.16  3 
3.44  3 

3.34  4 

2.34  1 


4.16  6 


3.86  4  0.31 

5.41  2  0.25 

3.29  4  0.20 

3.57  2  0.09 

2.87  2  0.16 

4.05  7  0.26 

3.45  11  0.25 

3.14  9  0.22 

2.23  3  0.15 

3.65  11  0.19 

3.75  9  0.22 

2.39  5  0.24 

3.51  13  0.14 

2.56  6  0.12 

3.31  5  0.22 

3.52  9  0.21 

2.98  4  0.17 

3.72  5  0.10 

3.67  6  0.15 

2.93  8  0.20 

2.96  9  0.11 

3.68  11  0.26 

2.44  4  0.20 

3.78  9  0.29 

3.12  11  0.24 

1.95  3  0.21 

3.36  5  0.28 

2.99  7  0.12 

2.78  10  0.22 

3.43  7  0.23 

3.8  7  0.12 

3.17  6  0.25 

3.75  5  0.20 

3.96  9  0.22 

3.59  8  0.22 

3.34  11  0.32 

3.68  10  0.35 

3.71  4  0.29 

3.66  13  0.18 

4.01  14  0.24 

2.97  12  0.12 


3.51  13 
2.56  6 
3.31  5 

3.52  9 
2.98  4 
3.72  5 

3.67  6 
2.93  8 
2.96  9 

3.68  11 


2.78  10 


3.59  8 
3.34  11 


2.25  1 
3.61  7 

2.99  10 


3.90  13 

3.79  13 

2.80  8 
3.46  4 
2.94  4 
3.39  9 
4.00  13 
3.07  10 


Table  6.2.4.  Origin  time  and  locations  for  the  42  seismic  events  after  relocation  using  the  BAREY 
model.  Magnitudes  (M^)  were  calculated  using  the  attenuation  parameters  determined 
through  inversion,  using  amplitude  observations  in  the  3-6  Hz  frequency  band.  Mean 
magnitudes,  no.  of  observations  and  RMS  of  residuals  are  listed  for  P„  and  S„  phases,  and 
also  network  magnitudes,  which  are  the  mean  of  all  reported  station  magnitudes  (P  and  S  } 
for  each  event.  Automatic  network  magnitudes  (M^J  from  NORSARs  GBF  system  (Ringdal\ 
Kvcema,  1989)  and  the  reference  magnitudes  used  to  calibrate  the  magnitude  scale 
(Ringdal  &  Kremenetskaya,  1999)  are  also  shown.  The  two  events  marked  in  italic  (1988-235 
and  1995-162)  could  not  be  used  as  input  to  the  inversion  determining  the  attenuation  relation 
for  and  S„. 
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Fig.  6.2.1.  Events  (circles)  and  seismic  stations  used  for  deriving  wave  propagations  characteristics 
of  the  Barents  Sea  and  surrounding  areas.  Array  stations  are  shown  as  squares,  while  3C  sta¬ 
tions  are  shown  as  triangles.  The  event  locations  shown  are  after  analysis  and  relocation 
using  the  BAREY  model  (Schweitzer  and  Kennett,  2002).  The  symbol  size  for  the  events  are 
proportional  to  the  network  magnitudes  in  Table  6.2.3. 
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Fig.  6.2.2.  Phase  observations  (dark  circles)  for  (upper  left),  P  (lower  left),  S„  (upper  right)  and 

Lg  (lower  right)  by  the  ARCES  array  (white  square).  The  white  circles  represent  events  where 
rfe  phase  in  question  is  not  observed,  but  at  least  one  other  phase  is  observable. 
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Fig.  6.2.3.  Mykkeltveit  and  Ringdal  (1981)  and  BAREY  (right)  (Schweitzer  and  Kennett,  2002)  mod¬ 
els.  For  the  Fennoscandia  model,  P„  and  Lg  arrivals  have  mostly  negative  residuals  while  the 
S„  arrivals  have  positive  residuals.  The  observations  are  distributed  in  a  fairly  symmetrical 
manner,  and  residuals  are  smaller  (particularly  for  S  phases)  when  the  BAREY  model  is  used. 
The  travel  times  are  reduced  using  8.0  km/sfor  P  and  4.5  km/sfor  S.  Note  that  the  vertical 
scales  are  different  for  the  P  and  S  plots. 
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Travel  Time  Residuals  (Fennoscandia)  Travel  Time  Residuals  (BAREY) 


RMS  (s)  RMS  (s) 

Fig.  6.2.4.  RMS  of  travel  time  residuals  for  P  (top)  and  S  (bottom)  phases  for  individual  events 
located  using  the  Fennoscandia  (left)  and  BAREY  (right)  models,  darker  symbols  indicate 
higher  residuals.  Events  with  larger  residuals  than  4.0  (P  phases)  or  6.0  (S phases)  seconds 
are  blacLThe  residuals  are  generally  smaller  for  locations  using  the  BAREY  model,  in  partic¬ 
ular  for  the  S  phases. 
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2345  2345 
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Fig.  6.2.5.  Magnitudes  calculated  using  the  attenuation  relation  of  Jenkins  et  al.  (1998).  Top:  Com¬ 
parison  between  different  frequency  bands  at  ARCES.  Note  the  systematic  trend  towards 
lower  magnitudes  with  increasing  frequency.  Lower  left:  Comparison  between  magnitudes 
calculated  for  P„  and  S„  phases  in  the  3-6  Hz  band  at  ARCES.  The  systematic  offset  is  almost 
0.5  magnitude  units.  Lower  right:  Comparison  between  P„  magnitudes  in  the  3-6  Hz  band  at 
ARCES  and  SPITS.  The  scatter  (a  =  0.423)  is  quite  large,  with  differences  of  up  to  a  full  mag¬ 
nitude  unit  for  some  events. 
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Fig.  6.2.6.  Magnitudes  calculated  using  the  attenuation  coefficients  and  station  corrections  from  the 
inversion.  Top:  Comparison  between  different  frequency  bands  at  ARCES.  The  systematic  fre¬ 
quency  dependency  as  visible  in  Fig.  6.2.5  is  significantly  reduced.  Lower  left:  Comparison 
between  magnitudes  calculated  for  P„  and  S„  phases  in  the  3-6  Hz  frequency  band  at  ARCES. 
There  is  no  obvious  systematic  offset,  although  the  scatter  is  marginally  increased  compared 
to  Fig.  6.2.5.  Lower  right:  Comparison  between  P„  magnitudes  in  the  3-6  Hz  frequency  band 
at  ARCES  and  SPITS.  The  average  offset  and  scatter  is  reduced  compared  to  Fig.  6.2.5. 
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Pn-S„ 


Fig.  6.2.7.  Network  event  magnitude  comparisons  and  maps  of  the  geographical  distribution  of  the 
magnitude  differences  for  vs.  5„.  Note  that  S„  magnitudes  are  overestimated  compared  to 

Pnfor  events  that  have  paths  predominantly  within  the  Baltic  shield,  while  events  with  paths 
that  cross  the  Barents  Sea  have  larger  magnitudes  compared  to  S„. 
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6.3  Site-Specific  Threshold  Monitoring  (SSTM)  applied  to  the  Lop  Nor 
test  site 

Introduction 

Continuous  seismic  threshold  monitoring  (TM)  is  a  technique  that  has  been  developed  at  NOR¬ 
SAR  over  the  past  decade  to  monitor  a  geographical  area  continuously  in  time.  Data  from  a 
network  of  arrays  and  single  stations  are  combined  and  “steered”  toward  a  specific  area  to  pro¬ 
vide  a  continues  assessment  of  the  upper  magnitude  limit  of  seismic  events  that  might  have 
occurred  in  that  area.The  basic  principles  have  been  described  by  Ringdal  and  Kvaema  (1989, 
1992),  who  showed  that  this  method  could  be  useful  as  a  supplement  to  event  detection  analy¬ 
sis.  The  usefulness  of  the  TM  method  for  global  network  capability  estimation  has  been  shown 
by  Kvsma  and  Ringdal  (1999).  Examples  of  site-specific  threshold  monitoring  (SSTM) 
applied  to  the  former  Soviet  test  site  at  Novaya  Zemlya,  and  the  Indian  and  Pakistani  test  areas 
have  been  demonstrated  by  Kvaema  et  al.  (2002a,b).  We  have  also  recently  reported  an  applica¬ 
tion  of  this  technique  to  the  site  of  the  “Kursk”  accident  in  the  Barents  Sea. 

The  main  purpose  of  the  threshold  monitoring  technique  is  to  highlight  instances  when  a  given 
threshold  magnitude  is  exceeded,  thereby  helping  the  analyst  to  focus  on  those  events  truly  of 
interest  in  a  monitoring  situation.The  analyst  can  then  apply  traditional  tools  in  detecting, 
locating  and  identifying  the  source  of  the  disturbance. 

In  this  paper  we  apply  the  SSTM  technique  to  the  Lop  Nor  test  site  in  China.  The  emphasis  will 
be  put  on  detection  (and  location)  of  small  seismic  events  with  mb  <  4.0,  and  the  purpose  is  to 
evaluate  the  SSTM  method  as  a  potential  monitoring  tool.  In  contrast  to  most  previous  case 
studies,  which  have  been  based  on  recordings  by  seismic  arrays  at  regional  distances,  we  will 
in  this  study  apply  a  combination  of  3-component  stations  and  arrays,  at  both  regional  and  tele- 
seismic  ranges. 

Our  efforts  so  far,  as  reported  in  this  contribution,  comprises  mainly  a  study  of  available  seis¬ 
mic  stations,  selection  of  those  stations  which  are  most  sensitive  to  seismic  events  in  the  Lop 
Nor  general  area,  and  tuning  of  the  signal  parameters  of  these  stations  so  as  to  prepare  process¬ 
ing  recipes  for  the  application  of  the  threshold  monitoring  tool. 

Development  of  processing  recipes 

For  the  successful  implementation  of  SSTM,  beams  filtered  in  optimal  frequency  bands  must 
be  steered  from  the  individual  arrays  towards  the  Lop  Nor  test  site,  and  amplitude  calibration 
constants  developed  from  older  events  are  applied  to  facilitate  the  calculation  of  continuous 
magnitude  thresholds  for  the  site.  For  single  (or  3-component)  stations,  the  vertical  component 
was  filtered  in  the  optimal  frequency  band. 

For  the  purposes  of  this  study,  23  stations  and  arrays  were  initially  chosen  as  shown  in  Fig. 
6.3.1.  For  each  of  the  stations  data  were  collected  from  known  nuclear  tests  at  Lop  Nor.  The 
explosions  used  for  the  calibration  is  summarized  in  Table  6.3.1,  however,  for  most  stations 
only  a  few  data  sets  with  explosions  were  available  for  the  calibration. 
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Event 

Event 

Event 

Event 

1990-146  M=5.5 

1993-278  M=5.9 

1995-135  M=6.0 

1996-211  M=4.7 

1992-142  M=6.5 

1994-161  M=5.8 

1995-229  M=5.9 

1992-269  M=5.0 

1994-280  M=5.9 

1996-160  M=5.8 

Table  6.3.1.  Lop  Nor  nuclear  test  explosions  used  for  the  station  calibration  purpose  (year,  Julian 
day  and  magnitude  of  the  explosion). 


Fig.  6.3. 1.  Stations  selected  for  SSTM  calibration  for  the  Lop  Nor  site  ( indicated  with  red  star). 


The  calibration  procedure 

The  development  of  the  calibration  parameters  followed  a  3 -step  procedure  for  the  arrays  listed 
in  Table  6.3.2: 

•  Firstly  the  average  azimuth  and  apparent  velocity  was  obtained  for  the  arrays.  Generally,  a 
higher  power  in  the  FK  domain  indicated  also  appropriate  frequency  bands. 

•  Through  beamforming,  the  STA  trace  aimed  at  Lop  Nor  was  computed  for  different  fre¬ 
quency  bands,  and  the  optimal  frequency  band  was  defined  based  on  maximum  SNR. 

•  For  the  optimal  beamforming  parameters,  the  maximum  corrected  STA  amplitude  for  the 
phase  was  computed  and  the  corresponding  magnitude  corrections  were  computed. 

The  process  is  exemplified  in  Fig.  6.3.2,  and  shows  the  distribution  of  FK  power  with  azimuth 
and  apparent  velocity  in  different  frequency  bands  for  the  array  CMAR.  The  final  parameters 
extracted  (and  later  applied  in  the  beamsteering)  from  the  data  are  apparent  velocity  and  azi¬ 
muth  (for  arrays),  optimal  frequency  filtering  band,  phase  travel  time,  phase  amplitude  and  tol¬ 
erance  bands  for  these  parameters.  Table  6.3.2  summarizes  some  of  the  parameters  obtained  for 
the  different  stations  and  arrays. 
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Station 

No.  of 
records 

Delta 

P- 

travel- 

time 

Calib. 

con¬ 

stant 

Calib. 
St.  devi¬ 
ation 

Optimal  Fre¬ 
quency  Band 

;  AAK  -v 

■11^9 

3.32607 

0.122042 

1.5 

abkt 

317.9 

3.20032 

0.131820 

0.5 

ARCES 

7 

42.47 

ARU 

6 

0.210574 

0.5 

2.0 

ASAR 

6 

0.209677 

1.5 

4.0 

BGCA 

1 

71.7 

687.2 

2.24374 

- 

1.0 

4.0 

BRVK 

4 

16.87 

239.6 

0.112567 

1.0 

4.0 

CMAR 

5 

24.64 

4.0 

FINES 

6 

472.4 

4.02934 

0.208297 

1.0 

4.0 

GBA 

^  29.60 

371.6 

4.25127 

0.178033 

1.0 

4.0 

GERES 

10 

51.41 

548.6 

4.42041 

0.241027 

2.0 

HFS 

3 

47.93 

0.063838 

1.0 

ILAR 

8 

65.57 

645.9 

3.69785 

0.112615 

0.5 

2.0  i 

INGN 

29.5 

- 

2.0 

4.0 

KURK 

1 

11.47 

165.7 

3.19758 

- 

1.0 

4.0 

KZA 

2 

9.92 

0.102790 

2.0 

LPAZ 

5  :;■> 

147.8 

1190.8 

2.83462 

0.050651 

1.5 

mmSH 

467.9 

2.45095 

0.991537 

1.5 

M^o: 

'■  2 

38.50 

443.3 

2.37779 

0.058921 

2.0 

lUlllli 

NIL 

4 

14.38 

0.067611 

1.0 

4.0 

NNA  ■ 

6 

147.9 

1189.8 

3.70520 

0.108466 

0.5 

2.0 

NORES 

9 

48.84 

528.5 

3.50340 

0.245675 

1.0 

TKM2 

1 

9.68 

146.0 

1.96420 

- 

0.5 

2.0 

ULHL 

2 

9.18 

0.082691 

1.0 

4.0 

ULN 

1 

14.64 

216.4 

1.60884 

- 

0.5 

2.0 

USP 

3 

10.51 

156.6 

2.44410 

0.189674 

1.0 

4.0 

YKA 

1 

74.6 

702.4 

4.77741 

- 

2.0 

6.0 

Table  6.3.2.  The  stations  calibrated  for  the  ACD  Lop  Nor  experiment  with  main  calibration  values. 
For  the  sJtaded  stations  calibration  parameters  were  developed,  but  these  have  not  yet  been 
implemented. 


78 


NORSAR  Sci.  Rep.  1-2002 


February  2002 


Frequency  Band 


Frequency  Band 


Frequency  Band 


Fig.  6.3.2.  Distribution  ofFK power, 
azimuth  and  apparent  velocity 
in  different  frequency  bands. 
Example  of  analysis  for  the 
CMAR  array  is  shown. 


Preliminary  Results 

Two  performance  tests  were  carried  out.  First,  data  for  one  day  (September  10,  2001)  was  col¬ 
lected,  and  in  these  data  the  recordings  of  4  explosions  were  scaled  and  embedded  at  some  (but 
not  all  of)  the  stations  (see  Table  6.3.3).  Secondly,  a  10  day  test  period  with  data  from  August  2 
through  11,  2001  was  selected  and  data  were  collected  for  the  stations.  Fig.  6.3.3  shows  exam¬ 
ples  of  network  results  of  these  two  performance  tests. 


Event 

Time  used 

Scaling 

magn. 

Stations  without  embedded,  scaled 
data 

1995/05/15 

09/10/2001  05:01:19 

mb  3.0 

1995/05/15 

09/10/2001  09:15:13 

mb  3.5 

1996/07/29 

09/10/2001  06:37:43 

mb  3.5 

ULN,  KURK 

1996/07/29 

09/10/2001  12:17:14 

mb  3.0 

ULN,  KURK 

1999/01127 

09/10/2001  13:25:08 

mb  3.5 

FINES,  GERES,  HFS,  NORES,  ULHL 

1999/01/30 

09/10/2001  14:23:03 

mb  3.5 

FINES,  GERES,  HFS,  NORES,  ULHL 

1999/01/27 

09/10/2001  16:59:57 

mb  3.0 

FINES,  GERES,  HFS,  NORES,  ULHL 

1999/01/30 

09/10/2001  20:59:58 

mb  3.0 

FINES,  GERES,  HFS,  NORES,  ULHL 

Table  6.3.3.  The  four  events/explosions  scaled  and  embedded  as  eight  events  in  the  data  for  day  253 
as  shown  in  Figs.  6.3.3-6.3.6. 
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Figs.  6.3.4  and  6.3.5  shows  the  SSTM  threshold  traces  for  single  stations  and  arrays  respec¬ 
tively.  In  both  figures,  the  data  cover  the  one  day  (10  September  2001)  which  contained  the 
scaled,  embedded  events,  and  the  combined  network  trace  (using  all  stations)  is  shown  on  top 
of  each  figure.  It  can  be  observed  that  the  arrays  generally  have  thresholds  between  magnitude 
3.2  and  3.8,  and  that  the  network  threshold  is  close  to  or  better  than  3.0  (based  on  both  arrays 
and  single  stations).  Some  stations  feature  frequent  data  problems. 

Fig.  6.3.6  shows  a  close-up  figure  of  the  detection  parameters  at  ARCES  for  an  hour  where 
there  is  also  an  embedded  event  (13:25).  It  should  be  observed  that  the  slowness  and  azimuth 
observed  from  the  event  are  within  the  bands  that  are  indicative  for  an  event  from  the  Lop  Nor 
area  (yellow  shaded  area). 

From  these  tests  the  following  preliminary  observations  can  be  stated: 

•  Out  of  the  eight  embedded  explosions  on  day  253  six  were  flagged,  and  the  two  non-flagged 
could  be  clearly  seen  in  the  data. 

•  The  quiet  day  (214)  did  not  have  any  events  flagged. 

•  On  day  217  the  SSTM  method  triggered  on  two  large  teleseismic  events. 

•  In  a  monitoring  situation,  all  peaks  exceeding  the  threshold  could  be  analyzed.  The  maxi¬ 
mum  number  of  such  peaks  during  one  day  was  10-15  for  the  days  processed.  If  such  peaks 
were  analyzed,  all  of  the  scaled  events  on  day  2001-253  would  have  been  found. 

The  above  observations  should  be  evaluated  under  the  perspective  that  the  calibration  is  a  pre¬ 
liminary  one,  where  the  majority  of  the  stations  could  only  be  calibrated  with  one  or  two 
events/explosions,  and  where  only  some  of  the  calibrated  stations  were  included  in  the  perfor¬ 
mance  test.  Furthermore,  the  day  with  embedded  events  did  not  include  embedded  data  for  all 
the  stations  (see  Table  6.3.3),  and  this  did  introduce  a  bias  in  the  threshold  monitoring  results. 

Conclusions 

From  the  above  we  conclude  that  the  Site-Specific  Threshold  Monitoring  performance  tests  for 
the  Lop  Nor  test  site  were  successful.  It  is  expected  that  these  initial  tests  will  be  followed  by 
more  detailed  studies,  where  in  particular  the  calibration  parameters  will  be  more  firmly  estab¬ 
lished. 


Conrad  Lindholm 
Tormod  Kvaema 
Johannes  Schweitzer 
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Fig.  6.3.3.  The  network  traces  from  three  different  days.  Upper  trace  from  day  253  with  embedded 
Lop  Nor  event  recordings  scaled  to  magnitudes  3.0  and  3.5.  Red  circles  indicate  that  the 
SSTM  analysis  picked  the  events,  while  blue  circles  indicate  that  the  events  were  not  picked. 
Lower  traces  are  from  days  214  and  217.  Day  214  is  silent  with  no  detections,  while  217 
show  two  detections  from  large  earthquakes  (teleseismic  events  far  from  Lop  Nor). 
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10  September  2001  Day  253 


Fig.  6.3.4.  Network  and  array  TM  traces  for  one  day  with  embedded,  scaled,  events.  Red  circles 
indicate  that  the  explosions  were  picked  by  SSTM,  blue  circles  indicate  that  they  were  not 
picked. 
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10  September  2001  Day  253 

Fig.  6.3.5.  Network  and  single  station  TM  traces  for  one  day  with  embedded,  scaled,  events.  Red 
circles  indicate  that  the  explosions  were  picked  by  SSTM,  blue  circles  indicate  that  they  were 
not  picked. 
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Fig.  6.3.6.  Network  and  ARCES  TM  traces  for  one  day  and  one  hour  with  one  embedded,  scaled, 
event.  The  lower  traces  provide  additional  information  for  the  ARCES  station.  The  yellow 
shaded  areas  indicate  the  azimuth  and  slowness  bands  expected  from  a  Lop  Nor  event. 
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6.4  Monitoring  the  seismicity  of  the  Spitsbergen  Archipelago 

6.4.1  Introduction 

The  work  described  in  this  paper  is  a  part  of  the  KRSC  -  NORSAR  cooperative  activity  aimed 
at  a  detailed  study  of  seismicity  in  the  Spitsbergen  region.  Part  of  the  motivation  for  the  study  is 
to  improve  the  quality  and  availability  of  well-located  reference  events  (“ground  truth  data”) 
for  location  calibration  purposes,  and  to  develop  local  velocity  models  for  improved  location 
accuracy. 

Spitsbergen  and  the  adjacent  areas  are  parts  of  a  geologically  complex  region  with  moderate  to 
high  seismicity.  The  main  seismicity  in  the  area  is  associated  with  the  North- Atlantic  Ridge, 
and  especially  the  Knipovich  Ridge  situated  at  a  distance  less  than  400  km  from  the  archipel¬ 
ago  (Sundvor  and  Eldholm,  1979,  Mitchell  et.  al.,  1990).  In  addition,  some  coal  mines  are 
located  in  the  area  of  Spitsbergen,  causing  occasional  induced  seismicity.  In  this  study  we 
describe  some  of  our  observation  of  unusual  features  of  seismic  events  and  wave  propagation 
in  the  area  and  our  preliminary  attempts  to  develop  local  travel  time  models. 

6.4.2  Station  Network  and  Geology 

The  Spitsbergen  (Svalbard)  Archipelago  is  characterized  by  a  very  complicated  geological 
structure.  The  simplified  map  in  Fig.  6.4.1  shows  the  delineation  of  the  four  major  geologic 
units: 

-  local  Tertiary  basin  sediments; 

-  Carboniferous  through  Cretaceous  platform  cover  sequence; 

-  Devonian  basin  sediments; 

-  metamorphic  basement  rocks; 

The  complexity  of  geological  structure  and  underwater  topography  causes  difficulties  in  loca¬ 
tion  and  interpretation  of  seismic  events.  We  will  return  to  some  illustrations  of  this  problem 
later. 

Local  seismicity  has  been  observed  in  Spitsbergen  for  a  considerable  time  by  a  seismic  net¬ 
work  comprising  digital  and  analog  seismic  stations  (see  Fig.  6.4.1).  The  analog  stations  BRB 
and  PYR  were  operational  from  about  1976  to  1990,  and  were  then  closed  down.  Recently,  a 
digital  station  has  been  reinstalled  at  the  BRB  location.  The  station  KBS  operated  as  an  analog 
station  until  the  1980s,  and  is  now  a  digital  broad-band  station.  The  array  SPI  is  part  of  the 
International  Monitoring  System,  and  has  been  in  operation  since  1992. 

Maps  of  seismic  events  in  and  near  Spitsbergen  for  the  period  1964  -  2001  are  shown  in  Fig. 
6.4.2.  The  locations  for  years  1964-1998  are  taken  from  ISC  bulletins  whereas  those  for  1998- 
2001  are  taken  from  NORSAR  Reviewed  Regional  Seismic  Bulletin.  Even  though  the  NOR¬ 
SAR  regional  bulletin  covers  a  much  shorter  time  interval  than  the  ISC  bulletin,  the  observed 
seismicity  patterns  in  the  two  plots  are  remarkably  similar.  In  particular,  a  segment  of  the  mid- 
atlantic  ridge  to  the  west  can  be  observed,  and  pronounced  seismic  zones  in  Nordaustlandet 
(northeast)  and  Heerland  (southeast)  can  also  be  identified  on  both  maps. 
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The  temporary  seismic  network  shown  in  Fig.  6.4.1  was  installed  in  Spitsbergen  in  1979  (Bun- 
gum  et.  al,  1982)  and  operated  during  that  year.  A  map  of  seismic  events  detected  and  located 
by  this  network  is  shown  in  Fig.  6.4.3. 

6.4.3  Travel  time  model 

In  2001  a  digital  seismic  station  GBV  316B  was  installed  in  Barentsburg  (Kremenetskaya  et. 
al,  2001b).  A  large  number  of  mining  rockbursts  and  explosions  have  been  recorded  (distances 
2-5  km  from  the  station).  Based  on  these  events  we  have  estimated  P  and  S  velocities  in  the 
upper  layer  and  obtained  values  of  4.54  and  2.52  km/sec,  respectively.  Such  values  are  usual 
for  sediments. 

We  used  the  calibration  event  28.01.2001  with  known  coordinates  (78.066  N,  14.324  E)  to 
extend  the  model  to  larger  distances.  This  event,  which  has  been  described  by  Kremenetskaya 
et.  al  (2001b),  was  well-recorded  by  the  SPI  and  KBS  stations  (Fig  6.4.4). 

We  tried  to  modify  the  BARENTS  travel  time  model  adding  a  layer  of  sediments  with  the 
velocities  mentioned  above  and  fitting  the  depth  of  this  sedimentary  layer  as  well  as  the  Moho 
depth.  As  a  preliminary  result,  which  seems  to  provide  the  best  fit  to  the  data,  we  obtained  the 
following: 

Table  6.4.1.  Velocity  model  SPITSO 


Depth  (km) 

Vp  (km/s) 

Vs  (km/s) 

0-2 

4.54 

2.52 

2-10 

6.20 

3.44 

10-30 

6.70 

3.72 

30-55 

8.10 

4.50 

55-210 

8.23 

4.57 

>210 

Same  as  lASPEI  91 

This  model  is  similar,  although  not  identical  to  the  model  presented  in  Kremenetskaya  et.  al. 
(2001a).  We  re-located  the  event  using  the  BARENTS  model  and  the  SPITSO  model.  The 
results  are  shown  in  Table  6.4.2,  and  demonstrate  that  the  location  error  decreases  from  16  km 
using  the  Barents  model  to  less  than  1  km  using  the  SPITSO  model.  While  this  kind  of  perfor¬ 
mance  cannot  be  expected  on  a  general  basis,  it  is  nevertheless  encouraging. 

Table  6.4.2.  Re-locating  the  calibration  event  28.01.2001 


Model 

Latitude 

Longitude 

Error  (km) 

Ground  Truth 

78.066 

14.324 

0.0 

SPITSO 

78.072 

14.317 

0.7 

BARENTS 

77.93 

14.045 

16 
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6.4.4  Earthquakes  in  Heer  Land 

The  Norwegian  Government  has  recently  approved  the  operation  of  a  new  coal  mine  in  the 
Svea  area  close  to  the  Heerland  seismic  zone.  In  a  report  to  the  Norwegian  Parliament  (JD 
Report  No.  9  (1999-2000))  it  was  stated  that  “rockfalls”  occurred  several  times  in  the  Svea 
mine  during  experimental  operation  in  1997-1998.  To  clarify  the  situation  regarding  local  seis¬ 
micity  we  relocated  a  set  of  seismic  events  that  occurred  in  this  area  using  the  travel  time  model 
mentioned  above.  The  results  are  shown  in  Fig.  6.4.5  and  Table  6.4.3. 

Table  6.4.3.  List  of  relocated  Heer  Land  earthquakes 


Number 

DATE 

TIME 

Latitude 

Longitude 

ML 

1 

1997/04/28 

06.09:48 

77.7 

17.1 

1.8 

2 

1997/05/29 

06.42:41 

77.7 

17.8 

2.0 

3 

1998/02/03 

05.40:33 

77.8 

17.6 

1.5 

4 

1998/05/17 

19.55:57 

77.8 

17.6 

1.5 

5 

1998/09/10 

21.57:37 

77.8 

18.1 

1.8 

6 

1998/10/13 

11.10:45 

77.8 

18.0 

1.8 

7 

1998/11/11 

07.37:37 

77.8 

17.9 

1.8 

8 

1998/11/26 

05.50:42 

77.9 

17.9 

3.0 

9 

1998/11/26 

14.12:36 

77.7 

17.7 

1.9 

10 

2001/07/31 

11.59:20 

77.8 

17.1 

2.1 

11 

2002/01/20 

08.53:21 

77.9 

16.9 

1.6 

Although  the  location  errors  here  could  amount  to  several  kilometers  (see  below)  it  is  possible 
to  see  that  most  of  the  earthquakes  occurred  in  a  compact  zone.  To  contour  this  zone  more 
exactly  it  would  be  desirable  to  install  a  set  of  additional  seismic  stations.  This  would  also  be 
useful  in  recording  possible  rockbursts  associated  with  future  mining  operations. 

6.4.5  Limits  of  the  model  application  and  possible  sources  of  location  errors 

There  are  some  indirect  arguments  in  favor  of  the  SPITSO  model  in  regions  of  Svalbard  cov¬ 
ered  by  sediments.  Thus,  some  geologists  estimate  the  thickness  of  sediments  to  be  up  to  2  km. 
In  addition,  the  estimates  of  apparent  velocities  computed  by  sensors  of  the  SPI  array  for  very 
close  events  (2-5  km  from  the  array)  are  in  the  range  4.4  -  5.2  km/sec  which  is  more  or  less  in 
agreement  with  the  model. 

It  is  clear  that  in  a  place  with  such  complicated  geological  structure  as  Svalbard  one  should  be 
very  careful  when  using  1-D  models.  Thus,  it  is  impossible  to  use  the  1-D  model  SPITSO  to 
locate  Svalbard  events  with  more  distant  stations  like  ARCESS.  An  example  is  shown  in  Fig 
6.4.6. 
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One  more  striking  example  of  location  problems  is  the  earthquake  which  occurred  on 
29.12.2001  at  14.24  GMT  in  Nordaustlandet.  It  was  recorded  by  KBS  and  SPI,  and  the  record¬ 
ings  are  shown  in  Fig.  6.4.7. 

Multiple  onsets  of  P  waves  are  clearly  seen,  especially  for  the  KBS  recording,  where  one  can 
identify  up  to  4  different  P  onsets.  In  contrast  to  typical  Pn  and  Pg  onsets,  say,  for  Khibiny 
explosions  recorded  at  ARCESS,  the  amplitudes  of  the  P  onsets  for  this  event  are  very  differ¬ 
ent. 

To  explain  such  features,  careful  modelling,  probably  3-dimensional,  is  required.  We  might 
expect  that  if  the  P  wave  has  several  onsets,  then  the  S  wave  should  have  the  same,  too.  But  we 
have  not  been  able  to  find  them  in  the  waveforms.  The  reason  might  be  that  first  S  onsets  are 
masked  by  the  coda  of  the  P  waves.  It  means  that  the  error  of  estimation  of  time  difference 
between  P  and  S  onsets  could  amount  to  as  much  as  5-10  seconds  or  more.  And,  therefore,  the 
location  error  by  a  single  station  could  amount  to  many  tens  of  kilometers. 

The  result  of  locating  this  event  is  shown  in  Fig.  6.4.8.  Lines  of  location  do  not  cross  at  the 
same  point.  The  backazimuth  calculated  by  the  SPI  array  using  beamforming  is  significantly 
different  from  that  corresponding  to  the  estimated  location. 

For  weaker  events  these  first  weak  onsets  of  P  waves  might  be  unnoticeable,  which  of  course 
could  influence  the  location  accuracy. 

After  having  looked  through  recordings  for  events  in  typical  seismic  regions  of  the  Svalbard 
Archipelago,  we  have  found  that  such  weak  first  P  onsets  are  typical  for  almost  all  areas,  except 
some  events  occurring  in  the  ridge  zone  to  the  West  of  Svalbard. 

6.4.6  Anomalies  in  P/S  ratios 

The  complex  geological  conditions  in  the  Svalbard  area  manifest  themselves  in  extremely  high 
variability  of  P/S  ratios.  We  have  noticed  this  in  particular  for  events  occurring  in  the  ridge 
zone  and  around  Bear  Island. 

Thus,  for  example,  KBS  recordings  of  two  earthquakes  occurring  very  closely  to  each  other  are 
shown  in  Fig.  6.4.9.  (3.12.2001  at  1.45:49  GMT,  78.33  N,  8.24  E  and  3.01.2002  at  22.50:40 
GMT,  78.4  N,  7.6  E).  The  P/S  ratios  for  these  two  events  differ  by  about  a  full  order  of  magni¬ 
tude,  even  though  the  events  are  very  close  together  and  of  similar  magnitudes. 

6.4.7  Conclusions 

1 .  A  preliminary  1 -dimensional  travel  time  model  has  been  developed  for  Spitsbergen. 

2.  Taking  into  account  the  very  complicated  geology  and  underwater  topography  in  and 
near  Spitsbergen,  it  is  necessary  to  use  a  3 -dimensional  model  or  some  kind  of  regionali¬ 
zation  for  improving  locations; 

3.  Multiple  onsets  of  P  and  probably  S  waves  can  strongly  increase  location  errors.  Such 
onsets  are  not  unconmion  for  Spitsbergen  earthquakes; 

4.  The  complexity  of  the  geology  manifests  itself  also  in  a  great  variability  of  P/S  ratios  of 
earthquakes  in  the  Spitsbergen  region.  Study  of  this  variability  could  help  to  “calibrate” 
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the  P/S  ratio  as  a  discriminant  between  earthquakes  and  explosions  or  find  conditions 
under  which  such  discriminants  fail; 

5.  Study  of  all  the  unusual  phenomena  mentioned  above  would  require  the  development  of 
more  advanced  techniques  for  analysis  of  seismic  events. 

E.  Kremenetskaya 
S.  Baranov 

V.E.  Asming 

F.  Ringdal 
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Fig.  6.4.1.  Geology  of  the  Spitsbergen  Archipelago.  Location  of  seismic  stations  are  indicated. 
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Fig.  6.4.2  Seismicity  of  Spitsbergen  during  1964-1998  using  ISC  data  (top)  and  during  1998-2001 
using  NORSAR  Reviewed  Bulletins  (bottom). 
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Fig.  6.4.3.Seismic  events  detected  and  located  by  the  temporary  network  in  1979. 


p  28.01.2001  11.01:23.2  aT=10  sec,  1=0.025  sec/suple  11.02:43.2-) 


28-01.2001  CALIBRATION  EVIHT 


Fig.  6.4.4.  Recordings  of  the  calibration  event  28.01.2001. 
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Fig.  6.4.5.  Location  of  Heerland  earthquakes. 


Fig.  6.4.6.  Example  of  difficulty  in  obtaining  location  of  event  south  of  Spitsbergen  using  the  SPITSO 
model  (isolines  do  not  intersect). 
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r  29.12.2001  14.23:89.9  OT=10  sec,  1=0.025  set/sanple _ 14.25:09.9  -\ 


NORTH  EAST  EVENT 


Fig.  6.4.7.  Recording  of  the  29.12.2001  earthquake  in  Nordaustlandet. 
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8  12  16  20  24  28 


LOCATION  OF  THE  NORTH  EAST  EVENT 


Fig.  6.4.8.  Example  of  difficulty  in  obtaining  location  of  an  event  in  Nordaustlandet  using  the 
SPITSO  model  (isolines  do  not  intersect). 
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Fig.  6.4.9.  Example  of  differences  in  P/S  ratios  for  earthquakes  located  closely  together.  (3.12.2001 
at  01.45.49  GMT;  78.33N  8.24E;  and  3.01.2002  at  22.50.40  GMT;  78.4N  7.6E). 
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6.5  Comparison  of  location  procedures  -  The  Kara  Sea  event  of  16  August 
1997 

Introduction 

A  seismic  event  in  the  Kara  sea  to  the  east  of  Novaya  Zemlya  on  16  August  1997  has  been  the 
subject  of  considerable  discussion,  because  the  very  limited  set  of  stations  available  to  the 
International  Seismological  Centre  (ISC)  and  the  Reviewed  Event  Bulletin  (REB)  of  the  Exper¬ 
imental  International  Data  Centre  (EIDC)  are  inadequate  to  constrain  the  depth  effectively. 
Could  this  event  have  been  an  underwater  explosion? 

A  more  comprehensive  data  set  has  been  assembled  by  repicking  all  available  records  from  the 
event,  including  the  station  AMD  at  the  northern  end  of  the  Urals  which  lies  the  closest  to  the 
event  at  about  400  km  range.  The  distribution  of  seismic  stations,  which  observed  this  event  is 
shown  in  Fig.  6.5.1,  the  station  coordinates  are  listed  in  Table  6.5.1,  and  the  full  phase  list  with 
arrival  time,  slowness  and  azimuth  readings  are  presented  in  Table  6.5.2. 

Many  of  the  stations  lie  in  Fennoscandia,  but  azimuthal  control  is  improved  by  a  set  of  stations 
in  the  western  Barents  Sea  and  in  Russia.  The  total  data  set  contains  many  S  readings;  even 
where  these  have  significant  reading  uncertainties  (see  Table  6.5.2),  they  can  help  to  provide 
control  on  position  and  depth  through  the  sensitivity  of  the  differential  S-P  times  to  distance. 

A  sequence  of  location  experiments  have  been  undertaken  to  compare  the  results  of  using  both 
different  velocity  models  to  describe  the  travel  times  of  the  phases  and  also  to  make  a  compari¬ 
son  between  the  use  of  a  linearized  location  algorithm  (HYPOS AT  -  Schweitzer,  2001)  and  a 
fully  non-linear  seheme  (shakeNA  -  Sambridge  and  Kennett,  2001).  For  direct  comparisons 
between  the  two  methods  we  have  used  a  standard  least-squares  misfit  criterion,  but  have  also 
examined  the  influence  of  more  robust  choices  for  data  misfit  when  using  the  non-linear  loca¬ 
tion  scheme. 

Observed  data  and  tested  models 

We  have  employed  the  global  reference  models  JB  (Jeffreys  and  BuUen,  1940),  PREM 
(Dziewonski  and  Anderson,  1981),  IASP91  (Kennett  and  Engdahl,  1991),  and  AKl 35  (Kennett 
et  ai,  1995),  and  also  considered  two  regional  models  derived  from  the  Barents  model  of 
Kremenetskaya  et  al.  (2001).  The  barey  model  is  very  similar  to  the  Barents  model  with  a  P-  to 
S- velocity  ratio  close  to  1.78,  while  the  barez  model  has  identieal  P  velocities  as  barey  model 
but  faster  S  wavespeed  with  a  P-  to  S-velocity  ratio  of  1.72  (see  Fig.  6.5.2).  We  have  found  that 
the  barey  model  gives  a  good  representation  of  paths  from  the  Kara  Sea  to  Fennoscandia  but 
that  the  barez  model  is  more  suitable  for  paths  from  the  Kara  Sea  to  Spitsbergen,  Bear  Island, 
and  to  Northern  Siberia. 

All  of  the  listed  observations  in  Table  6.5.2  were  used  during  our  study.  However,  the  final 
inversions  were  made  without  some  of  these  data:  The  S  onsets  at  KTKl  and  TRO  were  not 
used  because  we  had  no  direct  P  observation  of  these  in  trigger  mode  operated  stations.  The  P 
onset  at  MOR8  was  not  used  because  we  suspect  an  absolute  timing  error  at  MOR8  (for  all 
inversion  results  the  residuum  of  this  P  onset  is  anomalously  large).  The  ISC  used  a  reading  for 
NB2,  which  is  a  reported  analysis  on  the  beam  calculated  for  the  whole  NORSAR  array. 
Instead  of  recalculating  this  beam,  we  analyzed  the  onsets  individually  on  all  available  3C 
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broadband  traces  (NAOOl,  NB201,  NC303,  and  NC602)  of  the  NORSAR  array.  However, 
NC602  is  co-located  with  the  centre  element  (NRAO)  of  the  NORES  array  and  therefore  we 
skipped  the  NC602  reading  for  our  final  analyses.  The  pIDC  published  in  the  REB  an  Sn-onset 
time  for  station  ARU.  We  requested  the  ARU  data  from  the  pIDC  to  reanalyze  also  these 
records.  However,  after  carefully  checking  the  3C  data  around  the  presumed  and  reported  onset 
time,  we  could  not  confirm  any  Sn  onset  from  the  theoretically  backazimuth  and  rejected  this 
reading.  However,  this  spurious  onset  was  also  listed  in  the  ISC  bulletin  but  never  used  as 
defining  observation  at  the  EIDC  or  at  the  ISC.  Our  final  inversion  therefore  employs  23  P- 
phase  and  13  S-phase  readings  from  23  different  seismic  stations  at  a  distance  range  of  about 
about  3°  to  21°. 

Results  with  the  Neighbourhood  Algorithm  shakeNA 

The  non-linear  location  procedure  using  the  Neighbourhood  Algorithm  (NA)  is  based  on  an 
exploration  of  the  4-D  hypocentral  parameter  space  to  find  models  with  good  fit  to  the  data.  At 
each  stage  the  algorithm  makes  use  of  all  the  prior  information  to  define  a  partition  of  the 
parameter  space  into  Voronoi  cells  around  each  sampled  point  describing  the  region  which  lies 
closest  to  that  point.  To  find  effective  location  estimates  the  algorithm  is  used  in  a  fairly 
focused  way  with  an  initial  sampling  of  9  points  in  the  entire  specified  volume  and  then  9  new 
points  in  each  iteration  randomly  sampling  the  current  two  “best”  Voronoi  cells,  i.e.  those  con- 
taimng  the  points  with  the  least  misfit.  In  this  mode  the  progress  of  the  non-linear  inversion 
resembles  the  convergence  of  a  cloud  of  points  towards  single  goal. 

The  NA  scheme  was  applied  to  a  zone  16  deg  across  in  longitude  and  8  deg  across  in  latitude 
and  extending  to  1 20  km  deep,  centred  in  the  Kara  Sea,  and  up  to  20  s  variation  in  origin  time 
was  also  allowed.  Despite  the  large  initial  domain  convergence  is  rapid  and  low  misfits  can 
readily  be  achieved  in  about  30  iterations.  The  misfit  levels  for  the  two  global  reference  models 
AK135,  IASP91  is  comparable  but  AK135  provides  a  slightly  better  fit.  A  variance  reduction  of 
about  30%  can  be  achieved  by  using  the  region-specific  models,  barey  and  barez. 

In  Fig.  6.5.3  we  show  the  central  portion  of  the  parameter  space  (a  20  km  square  about  the 
same  reference  point)  for  the  barez  model.  The  NA  algorithm  has  been  extended  here  to  60 
iterations  to  provide  a  further  exploration  of  the  zone  of  better  misfit,  the  group  of  black  models 
have  a  very  similar  level  of  misfit  and  define  a  ‘consistency  region’  indicated  with  the  gray 
polygons  which  indicates  the  variability  allowed  in  the  event  location.  This  consistency  region 
is  defined  with  the  aid  of  an  auxiliary  weighting  function,  based  on  analogues  from  statistical 
physics.  It  is  easier  to  derive  a  threshold  scheme  for  the  weighting  function  than  for  the  misfit 
itself,  because  the  shape  of  the  weighting  function  is  known  and  only  scaling  varies.  There  is 
only  a  marginal  difference  between  the  levels  of  misfit  for  the  solutions  indicated  with  black 
symbols.  The  solution  with  best  fit  to  data  is  indicated  with  the  open  diamond,  but  a  more  use¬ 
ful  quantity  is  generally  the  ensemble  weighted  average  indicated  with  the  gray  cross. 

We  believe  that  the  robust  statistic  tends  to  suppress  the  influence  of  outUers  and  in  doing  so 
finds  an  alternative  region  of  good  fit  with  a  different  combination  of  depth  and  origin  time. 
The  sensitivity  of  the  data  misfit  to  subareas  of  the  hypocenter  space  and  indeed  groups  of  solu¬ 
tions  with  very  similar  levels  of  misfit  can  exist.  This  is  illustrated  with  Fig.  6.5.4,  where  we 
show  the  central  portion  of  the  parameter  space  for  model  IASP9L  Clearly  subregions  of 
equally  good  misfit  levels  can  be  identified,  notably  in  depth. 
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The  epicentral  areas  of  all  four  models  tested  with  the  NA  algorithm  are  plotted  in  Fig.  6.5.5 
and  Table  6.5.3  lists  the  corresponding  numerical  results.  The  polygons  give  always  the  consis¬ 
tency  region  of  the  solution  for  each  model.  In  addition  two  plotted  two  solutions;  once  the  one 
with  the  absolute  minimal  misfit  as  open  symbols  and  once  the  mean  of  the  whole  assemble  of 
solutions  defining  the  ‘consistency  region’  as  small  filled  symbols.  The  shakeNA  locations  for 
the  standard  Earth  models  IASP91  and  AK135  are  very  close  and  clearly  separated  by  about  20 
km  from  the  solutions  using  the  regional  models  harey  and  barez,  the  latter  show  somewhat 
smaller  consistency  regions. 

For  a  quantitative  measure  of  the  quality  of  an  event  location  we  have  calculated  two  different 
parameters:  the  root-mean  square  (RMS)  of  all  defining  onset  times  and  a  weighted  mean  mis¬ 
fit  of  all  observed  parameters  (WMF).  The  WMF  values  are  calculated  by  calculating  the  mean 
of  the  residuals  weighted  by  the  standard,  deviation  of  the  observation.  RMS  and  WMF  are  also 
listed  for  each  hypoc^nter  determination  in  Table  6.5.3. 

The  presence  of  sets  of  solutions  with  very  similar  levels  of  fit  poses  problems  when  using  a 
linearized  inverse  scheme  because  the  algorithm  will  tend  to  be  pulled  towards  one  particular 
epicentral  subregion  or  depth  range  (e.g.  depending  on  where  it  was  initiated),  and  unhke  the 
non-linear  scheme  will  not  be  able  to  escape  to  find  a  better  point.  If  we  are  unlucky  there  is  the 
possibility  of  approaching  the  minimum  from  a  direction  such  that  the  next  iteration  overshoots 
its  target  and  then  returns  to  near  the  earlier  point  at  the  subsequent  iteration.  Such  oscillations 
are  difficult  to  tame. 

In  both,  the  non-hnear  and  the  linearized  inversions,  we  found  that  the  patterns  of  residuals 
when  using  either  the  harey  or  barez  models  reveal  an  inconsistency  for  S  waves  between  paths 
to  Spitsbergen  and  Siberia  one  side  and  Fennoscandia  on  the  other  side.  Either  one  group  or  the 
other  group  of  paths  is  well  fitted,  but  always  with  significant  residuals  for  the  second  group. 
This  suggests  the  use  of  different  models  for  the  two  sets  of  paths. 

Discussion  of  the  HYPOSAT  results 

To  get  the  results  of  the  HYPOSAT  inversions  comparable  with  the  shakeNA  results,  the 
HYPOSAT  inversions  were  restricted  to  a  strict  least-squares  scheme  without  using  special  sta¬ 
tions  corrections  or  travel-time  differences  between  the  phases  observed  at  one  stations.  All 
error  ellipses  for  the  epicenter  were  calculated  for  a  90%  significance  level  after  projecting  the 
corresponding  uncertainties  of  source  time  and  event  depth  in  consideration  of  the  covariance 
matrix  into  the  epicentral-error  space.  The  iteration  process  was  stopped  when  two  sequential 
solutions  were  closer  than  0.5  km  to  each  other.  To  force  a  stable  result,  the  program  itself 
removes  some  observations  during  the  iterations  of  the  inversion  when  the  residuals  became 
too  large.  Therefore  the  number  of  defining  parameters  changes  slightly  for  the  different  solu¬ 
tions. 

We  tested  the  same  models  as  for  the  shakeNA  scheme.  The  epicenters  are  plotted  together 
with  the  shakeNA  results  in  Fig.  6.5.5,  and  Table  6.5.3  again  lists  the  corresponding  numbers. 
The  HYPOSAT  epicenters  are  plotted  with  large  filled  symbols  and  their  corresponding  error 
ellipses  as  broken  lines.  The  error  ellipses  and  epicenters  for  the  standard  Earth  models  IASP91 
and  AK135  overlap  very  well  with  the  shakeNA  results.  However,  because  of  the  inadequate 
velocity  models,  for  both  models  the  inversions  were  oscillating  between  a  source  depth  in  the 
middle  crust  and  below  the  Mohorovicic  discontinuity  and  the  program  fixed  the  depth. 
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Because  of  this  source  parameter  reduction  we  also  get  a  dimension  reduction  of  the  covariance 
matrix  and  consequently  relatively  small  error  ellipses,  although  the  RMS  and  WMF  values  are 
relatively  large. 

For  the  models  barey  and  barez  the  data  were  consistent  enough  that  HYPOSAT  inverted  also 
for  the  event’s  depth.  As  in  case  of  shakeNA,  the  models  can  only  explain  a  subset  of  data, 
depending  on  the  actual  ray  paths.  For  barey  the  S-phase  observations  on  Spitsbergen  and  on 
the  Bear  Island  could  not  be  modeled  and  the  error  ellipse  become  very  large.  The  hypocenter 
itself  lies  within  the  shakeNA  results.  For  barez  all  residuals  of  the  S  observations  in  Fennos- 
candia  became  very  large  and  were  refused  as  defining  parameters.  Consequently,  the  data  set 
became  more  homogeneous  and  the  error  ellipse  smaller.  The  agreement  with  the  shakeNA 
solutions  for  barez  is  worse  but  both  the  consistency  region  and  the  error  ellipse  show  some 
overlap. 

The  data  set  for  this  Kara  Sea  event  shows  a  bimodal  distribution  in  observed  travel  times. 
Therefore,  the  original  HYPOSAT  code  was  extended  such  that  it  can  now  jointly  apply  two 
different  travel-time  tables  after  choosing  the  best  travel-time  table  for  each  path  (see  Table 
6.5.2).  For  a  final  solution  —  here  called  comb  (Fig.  6.5.5,  in  red)  —  we  applied  all  possibili¬ 
ties  of  HYPOSAT  (i.e.  station  elevation  corrections,  inverting  for  travel-time  diflferences,  and 
choosing  the  travel-time  table  with  respect  to  the  paths)  and  the  result  was  a  location  with  a 
very  small  error  ellipse  and  very  small  data  residuals. 

To  find  out,  how  sensitive  our  comb  location  is  for  data  from  the  nearest  station  AMD,  we  ran 
several  tests.  Starting  with  a  depth  of  0  km,  we  inverted  the  data  with  different  combinations  of 
AMD  contributions,  i.e.  without  any  AMD  observations,  only  with  the  Pn  or  only  with  the  Sn 
onset,  and  with  both  onsets  with  or  without  also  inverting  the  travel-time  difference  Sn-Pn.  The 
result  of  these  tests  was  that  only  if  we  include  both  the  Pn  and  the  Sn  onset  time  (in  any  com¬ 
bination)  at  AMD,  we  could  resolve  the  event  depth  as  about  20  km. 

Finally,  we  tested  if  the  source  depth  as  estimated  for  solution  comb  is  depending  on  the  start 
depth  of  the  inversion.  Therefore  a  start  solution  close  to  the  epicentral  area  with  start  depth  of 
0, 10,  20,  and  30  km  was  chosen  for  the  inversion.  In  all  cases  the  inversion  converged  to  a 
depth  between  15  and  20  km  and  using  the  start  depth  as  fixed  depth  the  result  for  20  km  had 
the  smallest  residuals. 

In  conclusion,  the  location  estimations  from  the  different  techniques  agree  quite  well  and  show 
some  overlap  of  their  confidence  regions  or  error  ellipses.  From  this  we  conclude  that  the  event 
cannot  be  shallower  than  10  km  and  is  most  likely  in  the  lower  crust  at  about  20  km  depth. 

Comparison  with  data  centre  solutions 

The  international  data  centres  pIDC  (REB),  NEIC  (PDE),  and  ISC  located  the  discussed  Kara 
Sea  event.  They  all  used  slightly  different  observations  and  published  different  solutions.  All 
these  solutions  are  listed  in  Table  6.5.4,  and  Fig.  6.5.6  shows  the  epicenters  (REB:  diamond, 
PDE:  hexagon,  ISC:  point)  with  their  published  error  ellipses  in  black  together  with  our  comb 
solution  for  all  available  data  in  red  (see  Table  6.5.3).  The  number  of  defining  observations  for 
the  ISC  location  is  relatively  large  because  the  ISC  used  in  its  inversion  the  data  from  the  arrays 
SPITS,  FINES,  and  NORES  twice  as  defining:  once  under  the  array-beam  name  and  once 
under  the  name  of  the  central  site  of  the  array  (SPAO,  FIAO,  and  NRAO). 
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The  precise  travel-time  curves,  location  algorithms  and  applied  data  uncertainties  of  the  data 
centre  solutions  are  unknown  to  us.  Therefore,  it  is  difficult  to  compare  our  solutions  with  the 
data  centre  solutions.  We  choose  the  following  approach  for  our  experiment:  to  use  only  those 
of  our  readings,  which  were  made  at  the  same  stations  as  listed  in  the  ISC  bulletin  (all  entries 
marked  with  a  ‘C’  in  Table  6.5.2).  Then  we  relocated  the  event  (see  Fig.  6.5.6  and  Table  6.5.1) 
with  HYPOS  AT  for  the  different  standard  Earth  models  JB  (Jeffreys  and  Bullen,  1940,  gray, 
point),  PREM  (Dziewonski  and  Anderson,  1981,  blue,  hexagon),  IASP91  (Kennett  and  Eng- 
dahl,  1991,  red,  triangle),  and  A  JO 35  (Kennett  et  al,  1995,  blue,  square)  and  for  the  three  mod¬ 
els  derived  for  the  Barents  Sea  region  {barey  (rhomb),  barez  (triangle),  and  comb  (point),  all  in 
green),  where  comb  is  again  the  combined  usage  of  barey  and  barez. 

Because  this  subset  of  data  have  no  resolution  for  the  event’s  depth,  we  used  here  a  fixed  depth 
of  18  km  which  came  from  our  comb  solution  for  the  whole  data  set.  In  addition  we  listed  in 
Table  6.5.4  the  number  of  defining  parameters  used  for  the  locations,  the  RMS  values  of  all  as 
defining  used  onset  times,  and  the  epicentral  distance  to  o\xt  comb  solution  from  table  Table 
6.5.3.  To  get  some  feeling  about  the  quality  of  the  different  solution  we  calculated  again  the 
WMF  values.  For  the  data  centre  solutions  we  used  the  data  uncertainties  as  given  in  Table 
6.5.2  but  the  residuals  as  published  in  the  bulletins.  In  the  REB  as  well  as  in  the  ISC  bulletin  an 
Sn  reading  from  station  ARU  is  listed  but  was  not  used  to  locate  the  event.  To  calculate  the 
WMF  value  for  the  REB  solution,  this  residual  of  -20.8  s  was  ignored. 

All  relocations  for  the  three  global  models  JB,  IASP91,  and  AK135  are  very  close  and  from  the 
data  centre  solutions  only  the  REB  solution  lies  within  their  definitely  larger  error  ellipses.  The 
relatively  small  error  ellipse  of  the  REB  solution  is  mainly  because  the  Sn  onsets  with  their 
large  residuals  were  not  used  as  defining  for  the  solution.  In  our  relocations  (including  PREM) 
all  onsets  were  used  and  contribute  with  their  residuals  to  the  larger  error  ellipses.  In  all  these 
solutions  the  P-type  onsets  are  more  or  less  equally  well  explained  but  not  the  Sn  onsets. 
Because  this  could  be  a  depth  effect  we  also  relocated  the  event  with  a  fixed  depth  at  0  km  and 
at  10  km.  For  all  four  global  models  (JB,  PREM,  IASP91,  and  AA'735)  the  residuals  were  small¬ 
est  for  the  fixed  depth  at  18  km. 

The  smallest  residuals  were  obtained  for  the  new  regional  models  (see  the  green  locations  on 
the  map).  The  worst  of  these  solutions  is  for  model  barez,  of  which  we  know  that  it  explains 
best  the  readings  from  stations  in  the  western  Barents  Sea,  at  AMD,  and  at  NRIS.  The  subset  of 
data,  used  in  this  experiment,  contains  such  readings  only  from  the  SPITS  array.  All  other  read¬ 
ings  are  from  Fennoscandia  or  the  Kola  Peninsula.  For  these  data  model  barey  explains  better 
the  Sn-onset  times.  Consequently  the  location  with  model  barey  has  smaller  residuals.  The 
residuals  and  the  error  ellipse  for  model  comb  are  the  smallest  and  this  solution  is  the  closest  to 
our  comb  solution.  This  is  a  clear  effect  of  using  path  specific  velocity  models. 

Conclusion 

This  study  has  shown  both  the  importance  of  S  wave  information  in  assessing  the  depth  of 
regional  events,  and  the  need  to  get  a  good  regional  velocity  model  for  both  P  and  S  in  order  to 
place  the  strongest  constraints  on  the  location  of  the  event. 

The  conclusions  from  our  experiment  of  comparing  the  different  data  centre  solutions  with  our 
results  are  that  using  only  a  limited  data  set  but  an  adequate  travel-time  model  one  can  locate 
the  event  in  the  Kara  Sea  relatively  close  to  our  comb  location.  In  this  case  the  depth  resolution 
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is  of  course  negligible.  The  relative  small  error  ellipses  are  a  problem  which  arises  when  using 
only  a  limited  number  of  data.  Then  the  data  errors  do  not  usually  follow  a  normal  distribution 
but  are  biased  in  one  direction  and  suggest  high  resolution  and  accuracy. 

The  location  estimates  for  the  whole  data  set  from  the  different  techniques  agree  quite  well, 
with  some  overlap  of  the  estimated  confidence  regions.  The  event  cannot  be  shallower  than  10 
km  and  is  most  likely  in  lower  crust  around  20-30  km  depth.  Such  deep  crustal  events  are  often 
attributed  to  the  long-term  effects  of  ice-unloading  from  the  last  glaciation  and  were  previously 
observed  around  Novaya  Zemlya  {e.g.  Marshall  et  ai,  1989).  That  this  event  was  not  an  under¬ 
water  explosion  is  also  confirmed  by  the  occurrence  of  an  aftershock  in  the  same  epicentral 
area  about  four  hours  later  (Ringdal  et  ai,  1997)  and  the  fact  that  the  observed  seismic  signals 
do  not  show  bubble  pulse  reverberations,  typical  for  underwater  explosions. 


Johannes  Schweitzer 

Brian  Kennett,  RSES,  Australian  National  University 
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Table  6.5.1.  Stations  and  their  coordinates  used  in  this  study  (see  also  Fig.  6.5.1).  A 
star  indicates  if  the  reading(s)  of  this  station  were  used  or  listed  by 
international  agencies  or  in  this  study. 


Station 

Latitude 

Longitude 

Elevation 

[“] 

["] 

[m] 

AMD 

69.7667 

61.6833 

200.0 

APAO 

67.6061 

32.9931 

200.0 

ARU 

56.4302 

58.5625 

250.0 

BJO 

74.5055 

19.1883 

18.0 

FINES 

61.4436  ^ 

1 

26.0771 

150.0 

HFS 

60.1342 

13.6956 

265.0 

JOF 

62.9182 

31.3124 

180.0 

KAF 

62.1128 

26.3061 

205.0 

KBS 

78.9256 

11.9417 

74.0 

KEF 

62.1672 

24.8703 

215.0 

KEV 

69.7553 

27.0067 

81.0 

KJN 

64.0853 

27.7119 

250.0 

KTKl 

69.0117 

23.2371 

340.0 

MOR8 

66.1713 

14.4411 

445.0 

NAOOl 

60.8442 

10.8865 

426.0 

NB201 

61.0495 

11.2939 

613.0 

NB2 

61.0397 

11.2148 

717.0 

NC303 

61.2251 

11.3690 

401.0 

NC602 

60.7353 

11.5414 

305.0 

NORES 

60.7353 

11.5414 

302.0 

NRIS 

69.0061 

87.9964 

498.0 

NUR 

60.5090 

24.6514 

102.0 

PKK 

60.0052 

24.5169 

10.0 

PVF 

60.5451 

25.8616 

45.0 

SDF 

67.4203 

26.3936 

276.5 

SPITS 

78.1777 

16.3699 

323.0 

SUF 

62.7192 

26.1506 

185.0 

TRO 

69.6325 

18.9281 

15.0 

VAF 

63.0422 

22.6715 

55.0 

PIDC  NEIC  ISC 
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Table  6.5.2.  All  reread  parameters  as  onset  times,  backazimuths  (BAZ),  and  ray 
parameters  (P)  and  their  standard  deviations  o  of  first  P-  and  S-type 
onsets.  Epicentral  distances  (A)  and  azimuths  (AZI)  were  calculated  with 
respect  to  the  comb  solution  of  Table  6.5.3.  To  find  this  and  all  other 
solutions  of  this  table,  all  with  an  ‘I’  marked  data  were  inverted.  For  the 
comb  solutions,  we  used  model  barey  for  all  travel-time  relevant 
parameters  except  for  the  observations  indicated  with  a  ‘2’  -  in  this  case 
model  barez  was  applied.  A  ‘C’  marks  data  used  for  the  comparison  of 
the  results  of  this  study  with  the  solutions  published  by  international 
data  centres  (see  Table  6.5.4). 
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Table  6.5.3.  The  new  hypocenter  solutions  for  the  16  August  1997  Kara  Sea  event  as 
plotted  in  Fig.  6.5.6.  Listed  are  the  source  parameters  and  the  number  of 
data  used  by  HYPOSAT  as  defining  for  onset  times  (OT),  backazimuths 
(BAZ),  ray  parameters  (RP),  and  travel-time  differences  (TTD).  In  the 
header  line,  the  maximum  possible  number  of  these  parameters  is  given 
in  parenthesis.  WMF  is  the  weighted  mean  misfit  for  afi  available  travel¬ 
time  observations  and  RMS  was  calculated  from  all  defining  onset  times. 
An  ‘F’  together  with  the  estimated  depth  indicates  that  the  depth  was 
fixed  by  the  program.  ‘M’  indicates  the  used  method:  NA-B  -  the 
shakeNA  with  smallest  misfit,  NA-M  -  the  mean  of  all  solutions  in  the 
consistency  region,  and  H  -  the  HYPOSAT  solution. 


Model 

Lat 

n 

Lon 

V] 

D 

[km] 

oD 

[km] 

Time 

OT 

(36) 

BAZ 

(7) 

RP 

(7) 

IVMF 

(36) 

RMS 

[s] 

EE 

[km^] 

M 

AK135 

72.382 

57.356 

30.4 

29.1  -31.6 

02:11:02.85 

36 

B 

7 

2.63 

4.61 

- 

NA-B 

AK135 

72.391 

57.371 

29.4 

29.1  -31.6 

02:11:02.65 

36 

B 

7 

2.64 

4.60 

- 

NA-M 

AK135 

72.366 

57.367 

34.4  F 

13.2  -  51.7 

02:11:03.26 

36 

5 

1 

2.64 

4.58 

252 

H 

IASP91 

72.365 

57.342 

35.0 

32.1  -  36.0 

02:11:03.37 

36 

B 

7 

2.83 

4.97 

- 

NA-B 

IASP91 

72.373 

57.360 

34.2 

32.1-36.0 

02:11:03.21 

36 

B 

7 

2.84 

4.99 

- 

NA-M 

IASP91 

72.288 

57.246 

37.3  F 

18.5  -  62.2 

02:11:04.12 

29 

5 

1 

2.89 

2.69 

247 

H 

BAREY 

72.383 

57.740 

30.3 

27.8-31.0 

02:11:03.76 

36 

B 

7 

1.21 

1.80 

- 

NA-B 

BAREY 

72.389 

57.753 

29.7 

27.8  -  31.0 

02:11:03.66 

36 

B 

7 

1.21 

1.79 

- 

NA-M 

BAREY 

72.394 

57.802 

39.8 

10.1 

02:11:04.88 

30 

B 

1 

1.35 

1.27 

1858 

H 

BAREZ 

72.427 

57.868 

21.5 

21.1-24.5 

02:11:02.91 

36 

B 

D 

1.66 

3.67 

- 

NA-B 

BAREZ 

72.427 

57.874 

22.2 

21.1-24.5 

02:11:03.01 

36 

B 

7 

1.68 

3.66 

- 

NA-M 

BAREZ 

72.342 

57.644 

15.9 

7.6 

02:11:02.80 

28 

B 

1 

1.56 

1.32 

402 

H 

COMB 

72.331 

57.616 

17.6 

5.2 

02:11:03.01 

35 

5 

1 

0.90 

1.27 

238 

H 
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Table  6.5.4.  Hypocenter  solutions  for  the  16  August  1997  event  in  the  Kara  Sea  as 
published  by  international  data  centres  and  the  event  locations 
performed  for  comparison  with  our  results.  Listed  are  source 
parameters  including  the  feed  depth,  the  number  of  as  defining  used  P- 
^d  S-onset  times,  backazimuths  (BAZ),  ray  parameters  (P),  and  travel¬ 
time  (^erences  (TTD).  In  parenthesis  the  maximum  possible  number  of 
for  this  solution  available  data  is  given,  respectively,  such  values  are  not 
available  for  NEIC  and  ISC.  RMS  is  calculated  from  the  defining  onset 
times,  WMF  is  the  weighted  misfit  for  all  possible  travel-time 
observations  (in  case  of  data  centre  solutions  for  all  onsets,  for  which 
residuals  were  published),  and  A  gives  the  epicentral  distance  between 
this  solution  and  the  comb  solution  of  Table  6.53.  EE  lists  the  size  of  the 
error  ellipse  as  plotted  in  Fig.  6.5.6.  For  further  details  see  text. 


Model 

Lat 

V] 

Lon 

f] 

D 

[km] 

Time 

D 

S 

BAZ 

D 

TTD 

RMS 

[s] 

WMF 

A 

[km] 

EE 

[km^] 

REB 

72.6484 

57.3517 

0 

02:10:59.9 

5(5) 

-(4) 

3(5) 

3(5) 

-(3) 

0.20 

1.77 

36.5 

179 

PDE 

72.835 

57.225 

10 

02:10:59.77 

7(7) 

-(3) 

-(?) 

-(?) 

-(3) 

1.42 

0.93 

57.8 

%3 

ISC 

72.6171 

56.9353 

10 

02:10:59.18 

13(13) 

-(9) 

-(?) 

-(?) 

-(6) 

1.30 

1.08 

39.3 

358 

J-B 

12.6272 

57.1404 

18 

02:10:57.27 
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Figure  6.5.1:  Seismic  stations  (triangles  and  diamonds)  which  observed  the  Kara  Sea  event  of 
16  August  1997  (red  point).  Triangles  (diamonds)  were  used  for  stations  for  which  model 
barey  (barez)  gave  the  best  results  when  calculating  travel  times  for  the  joint  model 
‘comb  Blue  symbols  show  stations,  which  contributed  with  parameter  data  to  the  buller 
tins  of  the  international  data  centres,  and  green  symbols  represent  stations  with  addition¬ 
ally  analyzed  data  (this  study).  The  circles  show  distances  from  the  event  in  steps  of  500 
km. 
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Barents  Sea  Velocity  Models 


Figure  6.5.2:  Plot  of  regional  velocity  models  used  in  this  study.  The  models  barey  (red)  and 
barez  (green)  were  both  derived  from  the  barents  model  (blue,  Kremenetskaya  et  al, 
2001). 
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Figure  6.5.3:  Progress  of  the  shakeNA  program  to  a  location  estimate  for  the  barez  model 

shown  as  a  map  view  and  depth  slices  around  the  central  portion  of  the  sampled  region.. 
The  misfits  are  colour  coded  and  the  region  of  acceptable  misfits  is  indicated  with  black 
symbols.  The  estimate  of  the  consistency  region  associated  with  these  acceptable  models 
is  indicated  with  grey  shading.  The  right  hand  panel  shows  the  misfit  distribution  and  the 
approach  to  convergence.  New  models  are  generated  as  the  iterations  proceed  and  not 
all  are  as  good  as  earlier  models. 

The  ensemble  properties  indicated  by  the  grey  marker  and  cross  give  a  good  indication 
of  the  location  and  its  uncertainities. 
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Figure  6.5.5:  Shows  the  NA  (polygons)  and  the  HYPOSAT  (error  ellipses)  inversion  results  for 
the  full  data  set  and  the  different  velocity  models  IASP91  (black),  AK135  (magenta), 
barez  (blue),  barey  (green),  and  ‘comb’  (red).  For  more  details  see  text  and  Table  6.5.3. 
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6.6  Some  results  derived  from  the  seismic  signals  of  the  accident  of  the 
Russian  submarine  Kursk 

Introduction 

In  several  studies  different  authors  have  used  the  observed  seismic  signals  from  the  accident  of 
the  Russian  submarine  Kursk  to  investigate  in  detail  the  concomitant  circumstances  of  this 
tragedy.  Especially  the  location  capabilities  of  the  seismic  networks  and  the  yield  of  explo¬ 
sions,  which  presumably  destroyed  the  submarine,  were  investigated  in  detail  (Ringdal  et  al, 
2000;  Koper  et  al,  2001a,  b;  Savage  and  Helmberger,  2001;  Northrop,  2001).  The  results  of 
these  studies  were  used  to  launch  different  theories  on  the  presumed  cause  and  sequence  of 
events  of  the  accident.  In  this  short  note,  focus  is  on  some  new  aspects  resulting  from  a  study  of 
the  seismic  signals. 

Relative  location  of  the  two  seismic  events 

From  analysis  of  ARCES  data  it  became  very  soon  clear  that  two  different  seismic  events 
occurred  about  2  minutes  and  16  seconds  apart  in  the  same  area,  which  later  was  confirmed  as 
the  Kursk  submarine  accident  area.  The  first  of  these  two  events  (Kursk-1)  was  about  two  mag¬ 
nitude  units  smaller  than  the  second  one  (Kursk-2),  which  had  a  local  magnitude  of  about  3.5 
(Ringdal  et  al ,  2000);  for  details  see  Table  6.6. 1 .  To  get  a  better  understanding  of  the  accident, 
the  relative  location  between  these  two  events  is  investigated. 

We  assume  that  the  sources  of  both  events  were  in  the  submarine  or  in  the  surrounding  water. 
In  this  case,  S-waves  can  only  be  generated  from  P  energy  converted  at  the  bottom  of  the  Bar¬ 
ents  Sea.  Events  occuiring  at  different  depths  in  the  water  will  have  approximately  the  same  S- 
P  time  difference  at  recording  stations.  Our  event  locations  can  therefore  only  provide  informa¬ 
tion  on  the  horizontal  positions  as  projected  down  to  the  sea  bottom.  By  measuring  with  high 
accuracy  the  time  difference  between  both  events  at  many  stations  and  for  different  phases  we 
should  in  principle  be  able  to  provide  information  on  the  relative  horizontal  position  between 
Kursk- 1  and  Kursk-  2.  Therefore,  we  will  calculate  relative  coordinates  of  Kursk- 1  with  respect 
to  the  well  located  Kursk-2  event  by  applying  the  master-event  location  technique. 

The  Kursk- 1  was  best  observed  at  ARCES  and  correlation  analysis  between  different  onsets  for 
the  two  events  show  correlation  coefficients  of  up  to  0.78  (see  Table  6.6.2).  This  indicates  sim¬ 
ilar  but  not  identical  propagation  effects  and  source  characteristics  for  Kursk- 1  and  Kursk-2. 
The  correlation  method  was  used  to  measure  the  travel-time  differences  between  the  two  events 
accurately.  The  first  event  was  less  visible  at  other  stations  than  at  ARCES,  but  using  the  corre¬ 
lation  method,  signals  from  the  first  event  could  also  be  identified  at  the  Apatity  array,  the 
FINES  array,  and  at  the  3C-broadband  stations  APZ9,  KEV,  and  LVZ.  For  this  analysis  all  data 
were  resampled  to  a  common  digitalization  rate  of  400  samples  per  second  and  the  time  differ¬ 
ence  of  all  interpretable  signals  was  measured  with  an  assumed  accuracy  of  0.005  s  (two  sam¬ 
ples).  To  take  into  account  the  inversion  the  different  signal-to-noise  ratios  of  the  Kursk- 1 
signals,  the  measured  time  differences  were  weighted  with  the  observed  correlation  coeffi¬ 
cients.  Both  P-type  and  S-type  onsets  were  used  to  measure  the  source-time  difference.  These 
time  differences  scatter  around  a  mean  value  of  135.76  s,  and  details  with  respect  to  each  sta¬ 
tion  and  phase  are  listed  in  Table  6.6.2. 

Then  the  master-event  location  technique  was  appUed  and  the  deviations  in  the  time  differences 
were  inverted  for  a  relative  horizontal  location  between  the  two  events.  For  this  the  12  mea¬ 
sured  time  differences  could  be  inverted  by  applying  the  generalized-matrix  inversion.  The  data 
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were  weighted  with  the  standard  deviations  calculated  from  the  corresponding  correlation  coef¬ 
ficient.  The  inversion  reduced  the  variance  of  the  residuals  by  42.7%.  Table  6.6.3  shows  the 
results  of  the  inversion,  where  the  horizontal  distance  between  Kursk- 1  and  Kursk-2  was  esti¬ 
mated  at  about  145  m. 

Observations  by  analyzing  the  seismic  signals  of  the  Kursk  events 

As  mentioned  before,  many  authors  have  tried  to  estimate  the  yield  of  the  explosions,  which  hit 
the  Kursk  submarine.  Koper  et  al  (2001a,  b)  and  Savage  and  Helmberger  (2001)  proposed  a 
simple  explosion-like  source  for  Kursk-2.  However,  it  was  also  observed  that  all  readable  first 
movements  of  the  first  P  signals  have  a  negative  motion  for  this  event  (Koper  et  al.  (2001a,b); 
Northrop  (2(K)1)).  As  an  example.  Fig.  6.6.1  shows  the  clear  negative  first  motions  on  the  array 
beams  of  ARCES  and  FINES.  This  does  not  fit  with  the  idea  of  an  explosion  source.  Koper  et 
al.  (2001a,  b)  rejected  Northrop’s  (2001)  argument  against  an  explosion  source  mostly  by 
pointing  to  the  presence  of  a  bubble  signal  as  observed  in  the  spectra  of  Kursk-2  in  the  fre¬ 
quency  range  between  1  and  6  Hz  (Fig.  6.6.2,  red  curve).  The  wide  amplitude  maximum 
around  9  Hz  is  interpreted  as  the  signal  of  the  surface  reflection  and  its  reverberations  in  the 
water  layer.  For  a  water  depth  of  about  1 15  m  as  in  the  Kursk  accident  area  we  can  expect 
destructive  interference  for  the  frequencies  6.3  and  12.6  Hz  and  constructive  interference  for 
the  frequencies  9.5  and  15.8  Hz.  This  is  exactly  what  is  seen  in  Fig.  6.6.2  (red  curve),  support¬ 
ing  our  presumption  that  Kursk-2  happened  close  to  the  sea  bottom.  Fig.  6.6.2  also  shows  the 
spectra  for  Kursk- 1  (blue  curve)  and  an  equally  long  noise  sample  observed  just  before  the 
Kursk  accident.  It  can  easily  be  seen  that  the  data  of  Kurek-1  are  very  close  to  the  noise  level. 
Any  modulation  due  to  a  bubble  pulse  is  not  visible;  maybe  a  slight  amplitude  increase  for  the 
frequency  range  of  the  surface  reflection  /  water-layer  reverberations  signal  can  be  seen. 

As  reported  by  Ringdal  et  al.  (2000)  and  in  several  press  reports,  the  Russian  navy  conducted  a 
series  of  underwater  explosions  in  the  Kursk  accident  area  during  autrnnn  and  winter  2000/ 
2001.  For  the  largest  of  these  explosions  with  a  magnitude  of  about  2.5  (i.e.  about  one  magni¬ 
tude  unit  smaller  than  Kursk-2,  see  Table  6.6.1)  we  can  clearly  observe  a  positive  first  motion 
at  ARCES  and  at  FINES  (Fig.  6.6.1).  In  addition,  we  cannot  identify  a  bubble-pulse  related 
modulation  of  the  lower  part  of  the  spectrum  (compare  the  green  with  the  red  curve  in  (Fig. 
6.6.2).  The  signals  from  the  surface  reflection  /  water-layer  reverberations  at  about  9  Hz  are 
now  more  clearly  visible  than  for  Kursk-2  but  we  cannot  see  the  outstanding  amplitude  mini¬ 
mum  at  about  6  Hz.  We  can  also  identify  further  amplitude  maxima  at  about  1 1  and  15  Hz.  By 
interpreting  these  frequency  modulations  as  signals  from  the  bubble  pulse  and  assuming  an 
identical  explosion  depth  as  for  Kursk-2,  we  are  led  to  a  relatively  small  yield  estimate  of  about 
15  kg  TNT  equivalent  for  this  explosion  (after  relations  published  in  Gittermann  et  al.  (1998)). 
This  is  far  too  small  to  produce  a  magnitude  2.5  event. 

We  observe  these  apparent  discrepancies: 

•  clear  negative  first-motion  onsets  for  Kursk-2  and  spectral  evidence  of  a  bubble  pulse  for  a 
proposed  simple  explosion  inside  of  a  submarine. 

•  clear  positive  first-motion  onsets  for  an  underwater  explosion  with  no  bubble  pulse. 

Because  both  signals  came  from  the  same  source  area,  propagation  effects  can  be  excluded  as 
cause  for  the  observed  discrepancies.  This  leads  to  the  conclusion  that  the  source  histoiy  of 
Kursk-2  is  quite  complex.  It  also  indicates  that  an  explosion  inside  a  closed  steel  container, 
possibly  still  in  a  gas  volume,  cannot  be  described  with  the  standard  model  of  an  explosion 
source  in  water  as  done,  e.g.,  for  the  Dead  Sea  explosions  in  1999.  Therefore,  a  source  function 
with  an  initially  implosive  signal  may  be  considered,  as  also  proposed  by  Northrop  (2001).  In 
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the  following,  the  size  of  an  imploding  volume  will  be  estimated,  which  would  explain  all 
observed  seismic  energy  of  this  event  (as  an  extreme  case). 

Following  Muller  (1973,  2(X)1)  the  seismic  moment  Mq  of  a  volume  change  AV  due  to  an 
explosion  (or  implosion)  can  be  modeled  byMQ=  (k  +  2p)  ■  AV,  with  the  Lame’s  parameters 
X  and  p.  From  the  observed  local  magnitude  3.5  of  Kursk-2  a  seismic  moment  can  be  deduced 
of  about  Mq  =  9  •  10^^  Nm.  The  modeled  implosion  source  can  be  described  as  the  collapse 

of  a  gas  volume  at  normal  atmospheric  pressure  inside  the  submarine  due  to  the  sudden  pres¬ 
sure  change  (about  100  m  water  column)  after  a  leakage.  For  water  and  gas  the  shear  modulus 
p  =  0 .  Then  X  becomes  identical  to  the  bulk  modulus  k  and  the  seismic  moment  can  be  written 

as  Mq=  X  -  AV . 

The  change  of  the  gas  volume  AV  due  to  the  pressure  change  from  one  atmosphere  to  the  pres¬ 
sure  at  about  100  m  water  depth  can  be  written  after  applying  Boyle-Mariotte’s  law  for  ideal 

gasses  asAV=  V  ■  with  the  relative  pressure  change  ^ ,  and  the  original  volume  V.  The 

equation  for  the  volume  to  be  collapsed  to  radiate  a  specific  seismic  moment  is  then: 

K  •  Ap 

For  100  m  water  depth,  the  relative  pressure  change  is  about  a  factor  of  10  and  the  bulk  modu- 
lus  for  sea  water  is  about  k  =  2. 14  •  10  N/m"^.  Putting  all  results  together,  we  get  for  the  col¬ 
lapsed  volume  a  value  of  approximately  V  =  4200  m^.  Following  published  specifications  of 
the  Kursk  (e.g.  Federation  of  American  Scientists;  http:www.fas.org/nuke/guide/russia/theater/ 
949.htm),  the  150  m  long  Kursk  submarine  had  a  submerged  water  displacement  in  the  range 

16  400  to  24  000  tons,  which  corresponds  with  a  total  volume  of  about  16  700  to  24  400  m^.  In 
the  case  that  the  energy  from  Kursk-2  was  radiated  only  by  one  imploding  volume,  about  20  to 
26%  of  the  whole  submarine  must  have  been  cataclysmic^ly  flooded  during  this  event.  This 
volume  corresponds  well  with  news  reports  that  about  the  first  third  (bow)  of  the  submarine 
was  heavily  damaged. 

Discussion 

The  application  of  the  master-event  analysis  between  the  two  seismic  events  connected  with 
the  Kursk  accident  suggests  that  the  submarine  moved  about  145  m  to  the  north-west  during 
the  135.8  s  between  the  two  events.  The  azimuth  of  this  movement  is  about  302°.  After  the 
accident  not  only  the  exact  position  of  the  Kursk  submarine  became  known  but  also  the  direc¬ 
tion  in  which  the  submarine  was  lying  on  the  sea  bottom  (Lind,  2002).  This  direction  was 
reported  as  288°.  This  is  in  good  agreement  with  our  results  about  the  relative  movement  of  the 
submarine  during  the  time  interval  between  Kursk- 1  and  Kursk-2.  How  much  the  correspond¬ 
ing  change  in  depth  was,  cannot  be  resolved.  Assuming  a  pure  horizontal  movement,  the  mini¬ 
mum  average  velocity  of  the  submarine  was  about  1.1  no/s  (or  2.2  knots).  In  case  the  first  event 
occurred  when  the  submarine  was  close  to  the  surface  and  the  second  event  occurred  when  the 
submarine  was  close  to  the  bottom  of  the  Barents  Sea,  the  depth  difference  would  be  about  1(X) 
m.  Then,  the  total  change  in  position  was  about  180  m  and  we  get  a  maximum  average  velocity 
of  about  1.3  m/s  (or  2.6  knots).  Both  extreme  average  values  for  the  velocity  of  the  submarine 
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during  the  accident  suggest  that  the  submarine  was  already  in  a  more  or  less  motionless  state 
when  it  was  hit  by  the  first  event. 

The  observed  differences  between  the  explosion  events  during  autumn  and  winter  2000/2001, 
and  the  proposed  simple  explosion-like  source  derived  for  Kursk-2  are  obvious.  The  spectral 
analysis  of  the  signals  firom  Kursk-2  and  the  observed  negative  first  motion  of  the  onsets  clearly 
indicate  a  very  complex  source  function. 

Johannes  Schweitzer 
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Table  6.6.1.  Source  parameters  for  the  analyzed  events  (three  first  lines).  The  depth  of 
all  events  are  unknown.  The  epicenter  of  Kursk-2  is  assumed  to  coincide 
with  the  known  location  of  the  Kursk  submarine  on  the  sea  bottom  after 
the  accident  (provided  in  the  fourth  line). 


Event 

Date 

Time 

Latitude 

n 

Longitude 

n 

Depth 

[km] 

Magnitude 

(NORSAR) 

Reference 

Kursk-1 

12.08.2000 

07.28.26.6 

69.6160 

37.5740 

? 

1.50 

This  study 

Kursk-2 

12.08.2000 

07.30.42.4 

69.6166 

37.5708 

B 

3.50 

This  study 

Explosion 

15.11.2000 

06.23.16.8 

69.703 

37.001 

? 

2.49 

NORSAR 

Submarine 

12.08.2000 

— 

69.6166 

37.5708 

0.115 

“ 

Lind 

(2002) 

Table  6.6.2.  The  table  shows  the  data  used  for  the  master-event  inversion: 

App.vel.  is  the  applied  apparent  velocity  for  this  observation, 

Cor.coeff  is  the  measured  cross-correlation  coefficient  between  Kursk-1 
and  Kursk-2  for  the  phases  considered, 

5t  is  the  measured  time  difference, 

6t-M  is  the  measured  time  difference  after  removing  the  mean  value. 

Residuum  is  the  time  difference  after  the  inversion,  and 

a-5t  is  the  assumed  uncertainty  of  the  measured  time  difference. 


Station 

Backazimuth 

Phase 

App.  vel. 

Cor. 

coeff. 

5t 

6t-M 

Residuum 

a-6t 

APAO 

219.938 

Pn 

7.91 

0.455 

135.770 

0.0088 

0.0057 

0.0110 

APAO 

219.938 

Lg 

3.56 

0.407 

135.755 

-0.0063 

0.0128 

0.0123 

APZ9 

216.610 

Lg 

3.56 

0.311 

135.805 

0.0438 

0.0150 

0.0161 

ARCES 

273.839 

Pn 

8.61 

0.656 

135.760 

0.0188 

0.0071 

0.0119 

ARCES 

273.839 

PnPn 

8.14 

0.755 

135.757 

-0.0013 

-0.0098 

0.0076 

ARCES 

273.839 

Sn 

5.72 

0.780 

135.743 

-0.0043 

-0.0103 

0.0066 

ARCES 

273.839 

Lg 

4.12 

0.774 

135.740 

-0.0183 

-0.0147 

0.0064 

FIAO 

214.939 

Pn 

8.07 

0.421 

135.780 

-0.0213 

-0.0204 

0.0065 

KEV 

276.374 

Pn 

8.62 

0.496 

135.747 

-0.0143 

-0.0104 

0.0101 

KEV 

276.374 

Lg 

3.07 

0.734 

135.731 

-0.0303 

-0.0291 

0.0068 

LVZ 

210.569 

Pn 

7.74 

0.675 

135.769 

0.0078 

0.0087 

0.0074 

LVZ 

210.569 

Lg 

3.71 

0.770 

135.778 

0.0168 

0.0181 

0.0065 
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Table  6.63.  Results  of  the  master-event  location  between  Kursk-1  and  Kursk-2  and 
associated  standard  deviations  a.  The  distances  are  given  relative  to 
Kurek-1. 


Kursk-1  >  Kursk-2  [km] 

a  [km] 

East-West 

-0.123 

0.021 

North-South 

0.076 

0.013 

Fig.  6.6.1.  Butterworth  band-pass  (1.5  -  8  Hz)  filtered  beams  of  Kursk-2  (Jcurs)  and  the  discussed 
explosion  in  the  Barents  Sea  (_exp)  as  observed  at  the  regional  arrays  ARCES  and  FINES. 
The  source  details  are  listed  in  Table  6.6.1. 
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Fig.  6.6.2.  Power  density  spectra  from  the  Kursk  main  event  (Kursk-2,  red),  the  Kursk  precursor 

event  (Kursk- 1,  blue),  the  magnitude  2.5  explosion  in  the  Barents  Sea  close  to  the  Kursk  site 
(green ),  and  a  noise  sample  measured  just  before  the  Kursk  accident  (magenta).  All  spectra 
are  mean  spectra  of  all  three  components  of  the  central  site  ARAO  of  the  ARCES  array.  The 
time  series  were  all  90  s  long,  starting  with  the  P  onset,  and  including  most  of  the  Lg  wave 
train.  The  data  were  not  filtered  or  processed  before  calculating  the  power  density  spectra 
using  the  Welch  method. 
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